{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "from collections import defaultdict\n",
    "\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "from Bio.PopGen.GenePop import Controller as gpc\n",
    "\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "ctrl = gpc.GenePopController()\n",
    "my_pops = [l.rstrip() for l in open('hapmap10_auto_noofs_2.pops')]\n",
    "num_pops = len(my_pops)\n",
    "pop_iter, loci_iter = ctrl.calc_allele_genotype_freqs('hapmap10_auto_noofs_2.gp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "num_pops = len(my_pops)\n",
    "fix_pops = [0 for i in range(num_pops)]\n",
    "num_reads = [defaultdict(int) for i in range(num_pops)]\n",
    "num_buckets = 20\n",
    "MAFs = []\n",
    "for i in range(num_pops):\n",
    "    MAFs.append([0] * num_buckets)\n",
    "for locus_data in loci_iter:\n",
    "    locus_name = locus_data[0]\n",
    "    allele_list = locus_data[1]\n",
    "    pop_of_loci = locus_data[2]\n",
    "    for i in range(num_pops):\n",
    "        locus_num_reads = pop_of_loci[i][2]\n",
    "        num_reads[i][locus_num_reads] += 1\n",
    "        maf = min(pop_of_loci[i][1])\n",
    "        if maf == 0:\n",
    "            fix_pops[i] += 1\n",
    "        else:\n",
    "            bucket = min([num_buckets - 1, int(maf * 2 * num_buckets)])\n",
    "            MAFs[i][bucket] += 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA7kAAAIYCAYAAACokCs/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XuYJVV97//3R0AholxkQC6jQxAvaBB1BERzQsQLl2PAxAskAhoMJoGoz4kkaPI7EKOR5EQxRsWQQASiAklUUPCCKDEoKoMicpE46iDDIAygiIoo8P39Uatl09PT0zO7u3d39fv1PP3sXatWVX336n377rVqVaoKSZIkSZL64CGjDkCSJEmSpOlikitJkiRJ6g2TXEmSJElSb5jkSpIkSZJ6wyRXkiRJktQbJrmSJEmSpN4wyZUkaUSSrEjyvLWs+/Uk1892TBNZVyxJHpPkx0k2ms24JEmaiEmuJGlkWpL38yTbjCu/MkklWTKayEavqv67qp4w6jhgzVjGJ+dV9b2q2ryq7htNhJIkPcAkV5I0at8FDhtbSPJrwGajC2d6JNl41DFIkrQQmeRKkkbtLOCIgeUjgTMHKyTZIsmZSVYnuSHJXyZ5SFv3yiSXJvn7JD9I8t0kBwxsu0OS85PckWR5kj8YWHdikn9P8m9J7kryjSSPT/LGJLcmuTHJC9ZjX//R9vUj4JUDZee0/X81yVPHPf49klyV5M5Wb9O2v32TrBzY/+IkH25tcHuSd0/UmOs6ZpInJbkkyQ+TXJPktwbWHZjk2rbdTUneMD6WJGcBjwE+1oYo/1mSJa3nfeMpttO57f95V4th6cD6P2/HvivJ9Un2m+hxSpK0Nia5kqRR+xLwyJZ8bQS8HPi3cXX+EdgC+FXgN+iS4lcNrN8LuB7YBvg74LQkaes+BKwEdgBeAvzNuMTpRXSJ9lbA14BP0X0+7gi8Gfingbrr2tfBwH8AWwIfGCj7d2Br4IPAR5NsMrDNy4D9gZ2B3YFXjm+g1i4fB24AlrTYzh5fb1wcaxyzHfdjwKeBbYE/AT6QZGwo8mnAa6rqEcBTgM+O33FVHQ58D3hRG6L8dxMcf13t9Fst/i2B84F3t8f5BOBY4JkthhcCKyZ5nJIkrcEkV5I0F4z15j4f+CZw09iKgcT3jVV1V1WtAN4OHD6w/Q1V9c/tnNAzgO2B7ZIsBp4D/HlV/ayqrgT+Zdy2/11Vn6qqe+kSw0XASVX1C7pEbEmSLae4r8uq6qNVdX9V3d3Krqiq/2j7ewewKbD3wDbvqqpVVXUHXQK6xwTtsyddwnhcVf2kHf/SSdpzbcfcG9i8Pb6fV9Vn6ZLnseHivwB2S/LIqvpBVX11kmNMaIrtdGlVXdj+X2cBYz3N9wEPazFsUlUrqurb6xuDJGlhM8mVJM0FZwG/S9eLeea4ddsAD6XrxRxzA11v5pjvj92pqp+2u5vTJYZ3VNVdk2x7y8D9u4HbBiZQGktUp7qvGyd4bL8sq6r7eaCHc43YgZ+2Y423mC6Rv3eCdRNZ2zF3AG5sZRM9ht8BDgRuSPJfSZ41xeMNmko7jX/MmybZuKqWA68HTgRuTXJ2ksG2kiRpnUxyJUkjV1U30E1AdSDw4XGrb6PrYXzsQNljGOjtncQqYOskj9iAbTdkXzXBdovH7rTziHdq+1ofNwKPWY/JrNZ2zFXA4rHzmZtfPoaquryqDqYbyvxR4Ny17H+ixzlmqDavqg9W1XPo/t8F/O1UtpMkaYxJriRprjgKeG5V/WSwsPWqngu8NckjkjwW+D+sed7uGqrqRuCLwNuSbJpk93acD0y+5bTu6xlJfrslqK8H7qE7D3l9fAW4GTgpycPb8Z+9Acf8MvAT4M/aObr70p2TfHaShyb5vSRbtGHOP6IbPjyRW+jOj17DMG2e5AlJnpvkYcDP6HrSvSyRJGm9mORKkuaEqvp2VS1by+o/oUvOvgNcSjeZ0ulT3PVhdJM1rQI+ApxQVRdtYJgbsq/z6M4p/gHdeam/3ZLIKWuJ/ouAx9FN+rSy7XO9jllVP6eb9OkAuh7y9wJHVNU323aHAyva7NB/CLxiLft/G/CXbYbmN0ywfkPb/GHASS2279P1KL9pCttJkvRLqZpsxJEkSdpQSU4EHldVa0sWe3FMSZLmEntyJUmSJEm9YZIrSZIkSeoNhytLkiRJknrDnlxJkiRJUm+Y5EqSJEmSemOqF5UfiW222aaWLFky6jAkSZIkSTPgiiuuuK2qFk3nPud0krtkyRKWLVvbJRMlSZIkSfNZkhume58OV5YkSZIk9YZJriRJkiSpN0xyJUmSJEm9YZIrSZIkSeqNOT3xlCSpX5Ycf8GoQxi5FScdNOoQJEnqNXtyJUmSJEm9YZIrSZIkSeoNk1xJkiRJUm+Y5EqSJEmSesMkV5IkSZLUGya5kiRJkqTeMMmVJEmSJPWGSa4kSZIkqTdMciVJkiRJvbHxqAOQNPOWHH/BqEMYuRUnHTTqECRJkjQL7MmVJEmSJPWGPbma8+yFtBdSkiRJmip7ciVJkiRJvWGSK0mSJEnqDZNcSZIkSVJvmORKkiRJknrDJFeSJEmS1BvrTHKTLE7yuSTXJbkmyeta+dZJLkryrXa7VStPknclWZ7kqiRPH9jXka3+t5IcOXMPS5IkSZK0EE2lJ/de4E+r6knA3sAxSXYDjgcurqpdgYvbMsABwK7t72jgFOiSYuAEYC9gT+CEscRYkiRJkqTpsM4kt6purqqvtvt3AdcBOwIHA2e0amcAh7T7BwNnVudLwJZJtgdeCFxUVXdU1Q+Ai4D9p/XRSJIkSZIWtPU6JzfJEuBpwJeB7arqZugSYWDbVm1H4MaBzVa2srWVjz/G0UmWJVm2evXq9QlPkiRJkrTATTnJTbI58J/A66vqR5NVnaCsJil/cEHVqVW1tKqWLlq0aKrhSZIkSZI0tSQ3ySZ0Ce4HqurDrfiWNgyZdntrK18JLB7YfCdg1STlkiRJkiRNi6nMrhzgNOC6qnrHwKrzgbEZko8EzhsoP6LNsrw3cGcbzvwp4AVJtmoTTr2glUmSJEmSNC02nkKdZwOHA99IcmUrexNwEnBukqOA7wEvbesuBA4ElgM/BV4FUFV3JPlr4PJW781Vdce0PApJkiRJkphCkltVlzLx+bQA+01Qv4Bj1rKv04HT1ydASZIkSZKmar1mV5YkSZIkaS4zyZUkSZIk9YZJriRJkiSpN0xyJUmSJEm9YZIrSZIkSeoNk1xJkiRJUm+Y5EqSJEmSesMkV5IkSZLUGya5kiRJkqTeMMmVJEmSJPWGSa4kSZIkqTdMciVJkiRJvWGSK0mSJEnqDZNcSZIkSVJvmORKkiRJknrDJFeSJEmS1BsmuZIkSZKk3jDJlSRJkiT1hkmuJEmSJKk3THIlSZIkSb2xziQ3yelJbk1y9UDZiUluSnJl+ztwYN0bkyxPcn2SFw6U79/Klic5fvofiiRJkiRpoZtKT+77gf0nKD+5qvZofxcCJNkNOBR4ctvmvUk2SrIR8B7gAGA34LBWV5IkSZKkabPxuipU1eeTLJni/g4Gzq6qe4DvJlkO7NnWLa+q7wAkObvVvXa9I5YkSZIkaS2GOSf32CRXteHMW7WyHYEbB+qsbGVrK19DkqOTLEuybPXq1UOEJ0mSJElaaNbZk7sWpwB/DVS7fTvw+0AmqFtMnEzXRDuuqlOBUwGWLl06YR1JGoUlx18w6hBGbsVJB406BEmSpEltUJJbVbeM3U/yz8DH2+JKYPFA1Z2AVe3+2solSZIkSZoWGzRcOcn2A4svBsZmXj4fODTJw5LsDOwKfAW4HNg1yc5JHko3OdX5Gx62JEmSJElrWmdPbpIPAfsC2yRZCZwA7JtkD7ohxyuA1wBU1TVJzqWbUOpe4Jiquq/t51jgU8BGwOlVdc20P5o5yOGNDm+UJEkaz++IfkfUzJnK7MqHTVB82iT13wq8dYLyC4EL1ys6SZIkSZLWwzCzK0uSJEmSNKeY5EqSJEmSesMkV5IkSZLUGya5kiRJkqTeMMmVJEmSJPWGSa4kSZIkqTdMciVJkiRJvWGSK0mSJEnqDZNcSZIkSVJvmORKkiRJknrDJFeSJEmS1BsmuZIkSZKk3jDJlSRJkiT1hkmuJEmSJKk3THIlSZIkSb1hkitJkiRJ6g2TXEmSJElSb5jkSpIkSZJ6wyRXkiRJktQbJrmSJEmSpN5YZ5Kb5PQktya5eqBs6yQXJflWu92qlSfJu5IsT3JVkqcPbHNkq/+tJEfOzMORJEmSJC1kU+nJfT+w/7iy44GLq2pX4OK2DHAAsGv7Oxo4BbqkGDgB2AvYEzhhLDGWJEmSJGm6rDPJrarPA3eMKz4YOKPdPwM4ZKD8zOp8CdgyyfbAC4GLquqOqvoBcBFrJs6SJEmSJA1lQ8/J3a6qbgZot9u28h2BGwfqrWxlaytfQ5KjkyxLsmz16tUbGJ4kSZIkaSGa7omnMkFZTVK+ZmHVqVW1tKqWLlq0aFqDkyRJkiT124Ymube0Yci021tb+Upg8UC9nYBVk5RLkiRJkjRtNjTJPR8YmyH5SOC8gfIj2izLewN3tuHMnwJekGSrNuHUC1qZJEmSJEnTZuN1VUjyIWBfYJskK+lmST4JODfJUcD3gJe26hcCBwLLgZ8CrwKoqjuS/DVweav35qoaP5mVJEmSJElDWWeSW1WHrWXVfhPULeCYtezndOD09YpOkiRpmi05/oJRhzByK046aNQhSNKMme6JpyRJkiRJGhmTXEmSJElSb5jkSpIkSZJ6wyRXkiRJktQbJrmSJEmSpN4wyZUkSZIk9YZJriRJkiSpN0xyJUmSJEm9YZIrSZIkSeoNk1xJkiRJUm+Y5EqSJEmSesMkV5IkSZLUGya5kiRJkqTe2HjUAUiSJEnShlhy/AWjDmHkVpx00KhDmHPsyZUkSZIk9YZJriRJkiSpN0xyJUmSJEm9YZIrSZIkSeoNk1xJkiRJUm+Y5EqSJEmSemOoJDfJiiTfSHJlkmWtbOskFyX5VrvdqpUnybuSLE9yVZKnT8cDkCRJkiRpzHT05P5mVe1RVUvb8vHAxVW1K3BxWwY4ANi1/R0NnDINx5YkSZIk6ZdmYrjywcAZ7f4ZwCED5WdW50vAlkm2n4HjS5IkSZIWqGGT3AI+neSKJEe3su2q6maAdrttK98RuHFg25Wt7EGSHJ1kWZJlq1evHjI8SZIkSdJCsvGQ2z+7qlYl2Ra4KMk3J6mbCcpqjYKqU4FTAZYuXbrGekmSJEmS1maoJLeqVrXbW5N8BNgTuCXJ9lV1cxuOfGurvhJYPLD5TsCqYY4vSdJCs+T4C0YdwsitOOmgUYcgSZrDNni4cpKHJ3nE2H3gBcDVwPnAka3akcB57f75wBFtluW9gTvHhjVLkiRJkjQdhunJ3Q74SJKx/Xywqj6Z5HLg3CRHAd8DXtrqXwgcCCwHfgq8aohjS5IkSZK0hg1OcqvqO8BTJyi/HdhvgvICjtnQ40mSJEmStC4zcQkhSZIkSZJGwiRXkiRJktQbw15CSJIkSQuMs3w7y7c0l9mTK0mSJEnqDZNcSZIkSVJvzOnhyt+46c6RD4dxKIokSZIkzR/25EqSJEmSesMkV5IkSZLUGya5kiRJkqTeMMmVJEmSJPWGSa4kSZIkqTdMciVJkiRJvWGSK0mSJEnqDZNcSZIkSVJvmORKkiRJknrDJFeSJEmS1BsmuZIkSZKk3jDJlSRJkiT1hkmuJEmSJKk3THIlSZIkSb0x60lukv2TXJ9keZLjZ/v4kiRJkqT+mtUkN8lGwHuAA4DdgMOS7DabMUiSJEmS+mu2e3L3BJZX1Xeq6ufA2cDBsxyDJEmSJKmnUlWzd7DkJcD+VfXqtnw4sFdVHTtQ52jg6Lb4BOD6WQtww2wD3DbqIOY523B4tuH0sB2HZxsOzzYcnm04PNtweLbh9LAdhzfX2/CxVbVoOne48XTubAoyQdmDsuyqOhU4dXbCGV6SZVW1dNRxzGe24fBsw+lhOw7PNhyebTg823B4tuHwbMPpYTsObyG24WwPV14JLB5Y3glYNcsxSJIkSZJ6araT3MuBXZPsnOShwKHA+bMcgyRJkiSpp2Z1uHJV3ZvkWOBTwEbA6VV1zWzGMAPmzdDqOcw2HJ5tOD1sx+HZhsOzDYdnGw7PNhyebTg9bMfhLbg2nNWJpyRJkiRJmkmzPVxZkiRJkqQZY5IrSZIkSeoNk9xJJHlxkkryxLb8kCTvSnJ1km8kubxNovW6JO8c2O6fknxmYPlPkrxrFI9h1JI8OsnZSb6d5NokFyZ5fJK7k1w58HdEq//jcdu/Msm7RxP93DBJG149rt6JSd7Q7r8/yXdb234zyQmjiX5uGLINv57kf5KcmWTH0TyCuSPJdkk+mOQ7Sa5Icll7r9w3ycdbnTVet0kuSbKgLl+wLkl+nGTJwPvhtUnel+SpA++Ndwy8lj+z7r322wSfEVsmuT1J2vKz2uf2Tm15i9aGD2mv6Ze08q2TfC3Jq2b/UYxWa5+zBpY3TrJ63Ot39bjP6N2SLG3ffx7a6u3S3gceOarHMkpTbMd3t/sPSXJGktPTWZFkm7buGe01/rTRPJLRS/Kogefa95PcNLB8QpJrklzVlvdq2/iZMqA9ry5NcsBA2cuSfDLJfa3trk7ysSRbtvVLxn8P6hOT3MkdBlxKNws0wMuBHYDdq+rXgBcDPwS+COwzsN0ewBZJNmrL+wBfmJWI55D2peMjwCVVtUtV7Qa8CdgO+HZV7THwd+ZIg52j1tGG63JcVe1B93w8MsnOMxjqnDUNbfhU4AnA14DPjX3BW4haW34U+HxV/WpVPYPu/XGn0UY27327vVZ3B3YDdhl7b6S7AsFxbfl5I41yDqqqHwLfB57Uivahe62OfSbvDXy5qu4f2ybJFnQTYJ5aVf86i+HOFT8BnpJks7b8fOCmcXXOGfcZfW1VLQM+D7yh1XkP8BdV9aPZCXvOmUo7jr1vvg/YBHh1DUyGk2R34D+Al1fV12Y+5Lmpqm4feM97H3Byu/9HwP7A06tqd+B5wI0jDHXOas+rPwTekWTTJA8H3gocA9zd2vcpwB2trPdMctciyebAs4GjeCDJ3R64eezDsqpWVtUP6D5QH59ks/bh+VPgSuDX2nb70CXCC81vAr+oqveNFVTVlfgGtT6mow03bbc/mc7A5pGh27A6J9N9mT5gXfV77LnAz8e15Q1V9Y8jjKk3qupeus+Kx406lnnmCzyQ1O4DnDxuefDzd3PgE8AHq+qUWYtw7vkEcFC7fxjwoSlu9ybg1Un+DNikqqa6XV9NpR3/AXgUcMTgjy10P8x8FDi8qr4yo1HOX9sDt1XVPQBVdVtVrRpxTHNWVV0NfAz4c+AE4Myq+va4apcBC2JUmknu2h0CfLKq/ge4I8nTgXOBF7Uu/7ePDS1pX0yuBJ5J+9UY+BKwT5Id6GaxXoiJ3VOAK9aybpdxQ6F+fTYDm0em3IZ0v+AN+n+tfCVwdlXdOpOBzmHDtOF4XwWeOK3RzS9PpmuDqXj5uLZ1WNk6JPkVYD/gG6OOZZ4ZHE31q8C/88DzbfxIqncAl7YfrRays4FDk2xKN4Lgy+PWv3zcZ/Rm8Mue878F3gb88axGPDetqx1/F3gGcGj7rjjoPODYqrp05sOctz4NLE53ytB7k/zGqAOaB/6K7nl3APB3gyvaCNP96EYI9Z5J7todRvfmRbs9rKpW0g1bfCNwP3Bxkv1anbFfkveh+5Xksnb/2SzMXtx1GT9c+b8nqet1rib2oDakG+IzaGy48qOB/ZLss+YuFrx1teF4mY2g5osk70l3zvLlE6w+Z1zbLpvt+OaRXdoPAV8ALqiqT4w6oHnmC3Q/Ku8MrKiqn9GNEt2cLsEY7CX7LHBwkm1HEOecUVVXAUvovutcOEGV8cOV7x5YdwBwC93Q+gVtCu34VeCxwJ4TrPsMXa/4RhOsE1BVP6Z7DR8NrAbOSfLKkQY1x1XVT4BzgLPGesCBzdpnzO3A1sBFo4pvNpnkTiDJo+iG5f1LkhXAcXS/aqaq7qmqT1TVccDf0PX4wgO/JD+LLsG9ju4DYEGej9tcQ/fmtD7uHnfO49bAbdMX0ryzIW34IO1D4hLgOdMR0Dw0dBsOeBrda3uhugZ4+thCVR1D96vwopFF1A9jP7Y8rapOHHUw801VfQvYCngR3ecvdKM3XgV8t70HjjkbOAW4MMkjZjXQued84O+Z+lBlkvxvYAvghXSjhX5lhmKbTyZrx28CL6NLzp48bt2x7fa9MxjbvFdV91XVJVV1Al2b/c6oY5oH7m9/Y+5uPzY/FngonpO7oL2Ebhz7Y6tqSVUtBr4L/K82/JgkD6EbmnJD2+aLdEOVF1XVre0E8NXAwSzcntzPAg9L8gdjBUmeSfciW5v/Al7R6m5G9+HwuZkMco7bkDZ8kCQbA3sB48/LWCimow2T5LV05wd9cvpDnBvSzfb7q5NU+SywaZI/Gijr1ZfcdDPwvmVEx06Sf03ygyTrfY5euhlbn9fun5jk36a43SVJXr2+x5sN6Wawnco5ypcBr+OBJPcy4PVM8PlbVe8ELgY+siETya1P206yj99L8ulh9jENTgfeXFVTGh7fPpPfDhzTtjkP+IsZjG++mLQdq+qLdKfCXJDkMQOr7qfrAX5CkjfPfJjzT5InJNl1oGgPHvjerfVUVXcCrwXekGSTUccz00xyJ3YY3Wysg/4TeD/wsTbd9lXAvcC7AdoEVKvpejrGXAZsC3x9huOdk1qi/2Lg+eku3XINcCKwijXPyX1t2+x1wG+3YRVfAv69qj4/ivjngnW04bqMnZN7Fd05fh+esUDnsGlow68D/0N3zv1vVtXPZyzYSbQE5udpl50YKL+yJQJLhj1GVW1eVd+ZZH3RjV75jXSXvPgKcAbdJBdaD+3Hp3vGFT+HbobWnapqouGNgl9JsnLg7/+08i8Ai3lgWPxldOfnTvgjc1X9Od3kc2e1H61nVVV9oKpeMNvHHRfDyqr6h7WsHn9O7j7A/wd8tKqubXVOpDsfdde17GNBWEc7jtX5ON25kp9sowXHyu+h6wz5rSQLondtPW0OnJHu8mpX0Y2QPHG0Ic1vbRbvr/PApLq9lYGZzCVJc1Q7deIe4N1jsxkn+TW6y088Hti5qlaMLMCeSPJ+YGVV/eU07GujqrpvLeueCvzzYDKb5BXAH1bVBp1a0J4jr66qzyQ5EXhcVb1iCttdAvxbVf3Lhhx3YD8bTzC5zlCSFLBrVS2fzv2ux/HXeEzr07aSpNGwJ1eS5o+zgCMGlo8EHnSN6SQHJflakh8lubF9IR9b9/Ik30nyyLZ8QJLvJ1nUln85NLQN231vkk+0YcxfSPLoJO9sw2m/mTbD/PhtB7Z/S7u/b+t1+7Mktya5OckhSQ5MN2vmHUnetLYH3fb1niQXJLkryZeT7NLWLWnH3nig/i+H3yZ5ZYv95CQ/bI9/n1Z+Y4vnyHGH3CbJRe1Y/5XksQP7fmJbd0eS65O8bFycpyS5MMlP6C5fNf6x7NBGB3wV2DFtGH2So4B/AZ7V2vuvJth2lySfTXJ7ktuSfCDJlmtrt3Hb7p3ki60Nvp5k30nq/n6S69r/+VODj39cvbG2PyrJ9+iGs096rCSvavu+q/0vXjNun8e158eqJL8/bt2B6Xp07kpyU5I3MIGB//k/JrmzPVf3G1i/RZLT2nFuSvKWtMl/xj1f7mAKvUZJfivJNe3xXpLkSQPrFif5cJLV7f/27oHjOKuuJM0Qk1xJmj++BDwyyZPal/KXA+PPDfwJXSK8Jd31G/8oySEAVXUO3TDOd6UbMncaXc/f6rUc72XAXwLb0PUiX0aXnG1D14P8jvWI/dF012zeEfi/wD/TnX//DODXgf+byc8HPoxuuN9WwHK6i9xP1V50w/YfBXyQbvKhZ9Jdj/YVwLvTzcQ75veAv6Z7nFcCHwBI8nC6WSk/SHcqymHAe/PgCWV+t8X2CGCiJOZDdENrH073//mbJPtV1Wl05+1d1oaNnzDBtqG7dMsOdNfYXMzUkrAdgQuAt9BN5vcG4D/TftwYV/cQumuh/jbdhGL/zbonJvqNFs8Lp3CsW4H/DTySbmKok9Ndoo8k+7f6zwd2BZ437jinAa+pqkfQXRrss5PEtBfwHbr/4QnAh5Ns3dadQXe60ePoJpN7AfDqCbbdlnU8z5I8nq59Xk/XXhfSndb00PYa/TjdOYRL6J77Z69lV5KkaWSSK0nzy1hv7vPpZu68aXBlm4XyG1V1f7u8xYfokpAxx9DNHn8J8LF2rtjafKSqrmiXZPkI8LOqOrMNwT2HLkGYql8Ab62qX9B90d8G+IeququqrqGbz2D3Sbb/cFV9pQ0d/QDdBCRT9d2q+teBuBfTTRRzT1V9Gvg5XcIz5oKq+nw7X+4v6HpXF9MlZyvavu6tqq/SzdfwkoFtz6uqL7T2/9lgEG0fzwH+vKp+VlVX0vXeHj6VB1FVy6vqohb3arofGaZy3chXABdW1YUtrovozl09cIK6rwHeVlXXtbb+G2CPtfXmNidW1U/aZWYmPVZVXVBV367Of9FdB3PsOukvA/61qq6u7jIYJ447zi+A3ZI8sqp+0Np/bW4F3llVv2g/7lwPHJRkO7pL4Ly+xXwrcDIPPj9tVVX9Y/sf373mrh/k5XTPl4vac/vvgc3orqywJ90PEse1Y/2svCaqJM0Kk1xJml/OoustfCXjhioDJNkryefa8Mg76XoHfzlZVVX9EPh3up6wt6/jWLcM3L97guXNmbrbB85PHUsc1md/3x+4/9P1PPb441BVkx37xrE71V1+5g66ZOWxwF5tWOoPk/yQrtf30RNtO4EdgDuq6q6BshvoevjWKcm2Sc5uQ2x/RNeLv826tmtxv3Rc3M+hmy18orr/MFDvDroe5MliHHzMkx4r3RD5L7Xh3j+kS37HHsMO4/Y1fhbV32n1b0g3jPxZk8R0Uz140pEbeOB/uAlw80B8/0TXazvR41mXHQbjrKr72/Y70v2YcsN0n6csSVo3k1xJmkeq6ga6S5odyMQzZn+Q7rqNi6tqC+B9dEkKAEkkhSXGAAAgAElEQVT2AH6frof3XdMY2k958OWEHr22itPsJ+12Oo+9eOxOG8a8Nd1s3DcC/1VVWw78bV5Vg5dUmmw2x1XA1nnw9Vkfw7je+Em8re1/96p6JF2vaSbfBFrcZ42L++FVddJa6r5mXN3NqrsMytoMPua1HivJw+h6vv8e2K6qtqQb3jv2GG5moO3p2uaBg1RdXlUH0yWkHwXOnSSmHZMMts1jeOB/eA+wzUB8j6yqwSHn6zMj5yoGLkfWjrmY7n96I/CYDJwvLkmaHSa5kjT/HAU8tw3pHO8RdL2FP0uyJ12vLwBJNqXr/XsT3fmQOyb542mK6Urgd5Ns1M6tnMow2qG1Ybs3Aa9ox/59YJchd3tgkueku4bqXwNfrqob6c6vfHySw5Ns0v6eOTjR0DpivZHukjZvS7Jpkt3p/pcfmGJcjwB+DPywnft63BS3+zfgRUle2Npo03STge00Qd33AW8cO8+4TdL00ikeZ13HeijwMLrL7d2b5AC682HHnAu8MsluSX6F7lxaWhwPTXdt2S3asOAfARPOXN1sC7y2/Y9eSnfO8IVVdTPdEOm3J3lkkoekm9BrQ5+v59INg94v3XUn/5Quif4i8BW6xP2kJA9vbfHsDTyOJGk9mORK0jzTzmlctpbVfwy8OclddBM8DfZ2vY3u8jintPNNXwG8JdNzncvXAS8CxobwfnQa9jlVf0CX8N0OPJm1XBt1PXyQLsG6g25irN8DaMOMX0B3/uYquiHUf0uXuE3VYXSTEK2iO8/5hHbe6lT8FfB04E66yZ2mdO3rllwfTPfjxmq6HsbjmOA7QFV9hO4xnd2GRF9Ndw7rlEx2rNZ+r6V7Tv6A7geY8we2/QTwTroJpZaz5sRShwMrWlx/SPf8XZsv001edRvd5FEvqarb27oj6BLua1sc/8HEQ7en8nivb3H8YzvWi4AXVdXP2/D8F9Gd7/09YCXdObySpBnmdXIlSVJvJHkl3azhG3S9YUnS/GdPriRJkiSpNzY4yW3nlnwl3YXer0m7cH2SnZN8Ocm3kpzTzmkiycPa8vK2fsn0PARJkiRJkjrD9OTeQzfxyVPprle4f5K96c7lObmqdqU71+WoVv8o4AdV9Ti6a9L97RDHliRJWkNVvd+hypK0sG1wktsu5P7jtrhJ+yvguXSTOACcARzS7h/clmnr9xs3vb8kSZIkSUMZ6tptSTYCrqCbOfA9wLeBHw5c+HwlD1xAfkfaBdar6t4kdwKPopuNcHCfRwNHAzz84Q9/xhOf+MRhQpQkSZIkzVFXXHHFbVW1aDr3OVSS26bH3yPJlnSXQpjoWoFj0zdP1Gu7xtTOVXUqcCrA0qVLa9mytV0lQ5IkSZI0nyW5Ybr3OS2zK1fVD4FLgL2BLZOMJc870V0LELpe3cUAbf0WdNcglCRJkiRpWgwzu/Ki1oNLks2A5wHXAZ8DXtKqHQmc1+6f35Zp6z9bXqRXkiRJkjSNhhmuvD1wRjsv9yHAuVX18STXAmcneQvwNeC0Vv804Kwky+l6cA8d4tiSJEmSJK1hg5PcqroKeNoE5d8B9pyg/GfASzf0eJIkSZIkrcu0nJMrSZIkSdJcMNTsypIkSVp4lhx/wahDGLkVJx006hAkrYU9uZIkSZKk3jDJlSRJkiT1hkmuJEmSJKk3THIlSZIkSb1hkitJkiRJ6g2TXEmSJElSb5jkSpIkSZJ6wyRXkiRJktQbJrmSJEmSpN4wyZUkSZIk9YZJriRJkiSpN0xyJUmSJEm9YZIrSZIkSeqNDU5ykyxO8rkk1yW5JsnrWvmJSW5KcmX7O3BgmzcmWZ7k+iQvnI4HIEmSJEnSmI2H2PZe4E+r6qtJHgFckeSitu7kqvr7wcpJdgMOBZ4M7AB8Jsnjq+q+IWKQJEmSJOmXNrgnt6purqqvtvt3AdcBO06yycHA2VV1T1V9F1gO7Lmhx5ckSZIkabxpOSc3yRLgacCXW9GxSa5KcnqSrVrZjsCNA5utZIKkOMnRSZYlWbZ69erpCE+SJEmStEAMneQm2Rz4T+D1VfUj4BRgF2AP4Gbg7WNVJ9i81iioOrWqllbV0kWLFg0bniRJkiRpARkqyU2yCV2C+4Gq+jBAVd1SVfdV1f3AP/PAkOSVwOKBzXcCVg1zfEmSJEmSBg0zu3KA04DrquodA+XbD1R7MXB1u38+cGiShyXZGdgV+MqGHl+SJEmSpPGGmV352cDhwDeSXNnK3gQclmQPuqHIK4DXAFTVNUnOBa6lm5n5GGdWliRJkiRNpw1OcqvqUiY+z/bCSbZ5K/DWDT2mJEmSJEmTmZbZlSVJkiRJmgtMciVJkiRJvWGSK0mSJEnqDZNcSZIkSVJvmORKkiRJknrDJFeSJEmS1BvDXCdXkiRJkkZmyfEXjDqEkVtx0kGjDmHOMcmdYb7wfOFJkiRJmj0OV5YkSZIk9YZJriRJkiSpN0xyJUmSJEm9YZIrSZIkSeoNk1xJkiRJUm+Y5EqSJEmSesMkV5IkSZLUGxuc5CZZnORzSa5Lck2S17XyrZNclORb7XarVp4k70qyPMlVSZ4+XQ9CkiRJkiQYrif3XuBPq+pJwN7AMUl2A44HLq6qXYGL2zLAAcCu7e9o4JQhji1JkiRJ0ho2OMmtqpur6qvt/l3AdcCOwMHAGa3aGcAh7f7BwJnV+RKwZZLtNzhySZIkSZLGmZZzcpMsAZ4GfBnYrqpuhi4RBrZt1XYEbhzYbGUrG7+vo5MsS7Js9erV0xGeJEmSJGmBGDrJTbI58J/A66vqR5NVnaCs1iioOrWqllbV0kWLFg0bniRJkiRpARkqyU2yCV2C+4Gq+nArvmVsGHK7vbWVrwQWD2y+E7BqmONLkiRJkjRomNmVA5wGXFdV7xhYdT5wZLt/JHDeQPkRbZblvYE7x4Y1S5IkSZI0HTYeYttnA4cD30hyZSt7E3AScG6So4DvAS9t6y4EDgSWAz8FXjXEsSVJkiRJWsMGJ7lVdSkTn2cLsN8E9Qs4ZkOPJ0mSJEnSukzL7MqSJEmSJM0FJrmSJEmSpN4wyZUkSZIk9YZJriRJkiSpN4aZXVmSJM2yJcdfMOoQRm7FSQeNOgRJ0hxmT64kSZIkqTfsyZWkKbIHzR40SZoufqb4maKZY0+uJEmSJKk37MmVJM0aey7suZAkaaaZ5EoLgImFiYUkSdJCYZKrOc8EzQRNkiRJmirPyZUkSZIk9YZJriRJkiSpN0xyJUmSJEm9YZIrSZIkSeqNoZLcJKcnuTXJ1QNlJya5KcmV7e/AgXVvTLI8yfVJXjjMsSVJkiRJGm/Yntz3A/tPUH5yVe3R/i4ESLIbcCjw5LbNe5NsNOTxJUmSJEn6paGS3Kr6PHDHFKsfDJxdVfdU1XeB5cCewxxfkiRJkqRBM3VO7rFJrmrDmbdqZTsCNw7UWdnKHiTJ0UmWJVm2evXqGQpPkiRJktRHM5HkngLsAuwB3Ay8vZVngrq1RkHVqVW1tKqWLlq0aAbCkyRJkiT11bQnuVV1S1XdV1X3A//MA0OSVwKLB6ruBKya7uNLkiRJkhauaU9yk2w/sPhiYGzm5fOBQ5M8LMnOwK7AV6b7+JIkSZKkhWvjYTZO8iFgX2CbJCuBE4B9k+xBNxR5BfAagKq6Jsm5wLXAvcAxVXXfMMeXJEmSJGnQUEluVR02QfFpk9R/K/DWYY4pSZIkSdLazNTsypIkSZIkzTqTXEmSJElSb5jkSpIkSZJ6wyRXkiRJktQbJrmSJEmSpN4wyZUkSZIk9YZJriRJkiSpN0xyJUmSJEm9YZIrSZIkSeoNk1xJkiRJUm+Y5EqSJEmSesMkV5IkSZLUGya5kiRJkqTeMMmVJEmSJPWGSa4kSZIkqTeGSnKTnJ7k1iRXD5RtneSiJN9qt1u18iR5V5LlSa5K8vRhg5ckSZIkadCwPbnvB/YfV3Y8cHFV7Qpc3JYBDgB2bX9HA6cMeWxJkiRJkh5kqCS3qj4P3DGu+GDgjHb/DOCQgfIzq/MlYMsk2w9zfEmSJEmSBm08A/vcrqpuBqiqm5Ns28p3BG4cqLeyld08uHGSo+l6ennMYx4zA+FJkqSFbMnxF4w6hJFbcdJBow5BkmbMbE48lQnKao2CqlOramlVLV20aNEshCVJkiRJ6ouZ6Mm9Jcn2rRd3e+DWVr4SWDxQbydg1WQ7+sZNd47811Z/6ZQkSZKk+WMmenLPB45s948EzhsoP6LNsrw3cOfYsGZJkiRJkqbDUD25ST4E7Atsk2QlcAJwEnBukqOA7wEvbdUvBA4ElgM/BV41zLElSZIkSRpvqCS3qg5by6r9JqhbwDHDHE+SJEmSpMnM5sRTkiRJkiTNKJNcSZIkSVJvmORKkiRJknrDJFeSJEmS1BsmuZIkSZKk3jDJlSRJkiT1hkmuJEmSJKk3THIlSZIkSb1hkitJkiRJ6g2TXEmSJElSb5jkSpIkSZJ6wyRXkiRJktQbJrmSJEmSpN4wyZUkSZIk9YZJriRJkiSpNzaeqR0nWQHcBdwH3FtVS5NsDZwDLAFWAC+rqh/MVAySJEmSpIVlpntyf7Oq9qiqpW35eODiqtoVuLgtS5IkSZI0LWZ7uPLBwBnt/hnAIbN8fEmSJElSj81kklvAp5NckeToVrZdVd0M0G63Hb9RkqOTLEuy7L6f3jmD4UmSJEmS+mbGzskFnl1Vq5JsC1yU5JtT2aiqTgVOBXjY9rvWDMYnSZIkSeqZGevJrapV7fZW4CPAnsAtSbYHaLe3ztTxJUmSJEkLz4wkuUkenuQRY/eBFwBXA+cDR7ZqRwLnzcTxJUmSJEkL00wNV94O+EiSsWN8sKo+meRy4NwkRwHfA146Q8eXJEmSJC1AM5LkVtV3gKdOUH47sN9MHFOSJEmSpNm+hJAkSZIkSTPGJFeSJEmS1BsmuZIkSZKk3jDJlSRJkiT1hkmuJEmSJKk3THIlSZIkSb1hkitJkiRJ6g2TXEmSJElSb5jkSpIkSZJ6wyRXkiRJktQbJrmSJEmSpN4wyZUkSZIk9YZJriRJkiSpN0xyJUmSJEm9YZIrSZIkSeqNWU9yk+yf5Poky5McP9vHlyRJkiT116wmuUk2At4DHADsBhyWZLfZjEGSJEmS1F+z3ZO7J7C8qr5TVT8HzgYOnuUYJEmSJEk9laqavYMlLwH2r6pXt+XDgb2q6tiBOkcDR7fFJwDXz1qAG2Yb4LZRBzHP2YbDsw2nh+04PNtweLbh8GzD4dmGw7MNp4ftOLy53oaPrapF07nDjadzZ1OQCcoelGVX1anAqbMTzvCSLKuqpaOOYz6zDYdnG04P23F4tuHwbMPh2YbDsw2HZxtOD9txeAuxDWd7uPJKYPHA8k7AqlmOQZIkSZLUU7Od5F4O7Jpk5yQPBQ4Fzp/lGCRJkiRJPTWrw5Wr6t4kxwKfAjYCTq+qa2Yzhhkwb4ZWz2G24fBsw+lhOw7PNhyebTg823B4tuHwbMPpYTsOb8G14axOPCVJkiRJ0kya7eHKkiRJkiTNGJNcSZIkSVJvmOROIsmLk1SSJ7blhyR5V5Krk3wjyeVtEq3XJXnnwHb/lOQzA8t/kuRdo3gMo5bk0UnOTvLtJNcmuTDJ45PcneTKgb8jWv0fj9v+lUnePZro54ZJ2vDqcfVOTPKGdv/9Sb7b2vabSU4YTfRzw5Bt+PUk/5PkzCQ7juYRzB1JtkvywSTfSXJFksvae+W+ST7e6qzxuk1ySZIFdfmCdUny4yRLBt4Pr03yviRPHXhvvGPgtfyZde+13yb4jNgyye1J0paf1T63d2rLW7Q2fEh7Tb+klW+d5GtJXjX7j2K0WvucNbC8cZLV416/q8d9Ru+WZGn7/vPQVm+X9j7wyFE9llGaYju+u91/SJIzkpyezook27R1z2iv8aeN5pGMXpJHDTzXvp/kpoHlE5Jck+SqtrxX28bPlAHteXVpkgMGyl6W5JNJ7mttd3WSjyXZsq1fMv57UJ+Y5E7uMOBSulmgAV4O7ADsXlW/BrwY+CHwRWCfge32ALZIslFb3gf4wqxEPIe0Lx0fAS6pql2qajfgTcB2wLerao+BvzNHGuwctY42XJfjqmoPuufjkUl2nsFQ56xpaMOnAk8AvgZ8buwL3kLU2vKjwOer6ler6hl07487jTayee/b7bW6O7AbsMvYeyPdFQiOa8vPG2mUc1BV/RD4PvCkVrQP3Wt17DN5b+DLVXX/2DZJtqCbAPPUqvrXWQx3rvgJ8JQkm7Xl5wM3jatzzrjP6GurahnweeANrc57gL+oqh/NTthzzlTacex9833AJsCra2AynCS7A/8BvLyqvjbzIc9NVXX7wHve+4CT2/0/AvYHnl5VuwPPA24cYahzVnte/SHwjiSbJnk48FbgGODu1r5PAe5oZb1nkrsWSTYHng0cxQNJ7vbAzWMfllW1sqp+QPeB+vgkm7UPz58CVwK/1rbbhy4RXmh+E/hFVb1vrKCqrsQ3qPUxHW24abv9yXQGNo8M3YbVOZnuy/QB66rfY88Ffj6uLW+oqn8cYUy9UVX30n1WPG7UscwzX+CBpHYf4ORxy4Ofv5sDnwA+WFWnzFqEc88ngIPa/cOAD01xuzcBr07yZ8AmVTXV7fpqKu34D8CjgCMGf2yh+2Hmo8DhVfWVGY1y/toeuK2q7gGoqtuqatWIY5qzqupq4GPAnwMnAGdW1bfHVbsMWBCj0kxy1+4Q4JNV9T/AHUmeDpwLvKh1+b99bGhJ+2JyJfBM2q/GwJeAfZLsQDeL9UJM7J4CXLGWdbuMGwr167MZ2Dwy5Tak+wVv0P9r5SuBs6vq1pkMdA4bpg3H+yrwxGmNbn55Ml0bTMXLx7Wtw8rWIcmvAPsB3xh1LPPM4GiqXwX+nQeeb+NHUr0DuLT9aLWQnQ0cmmRTuhEEXx63/uXjPqM3g1/2nP8t8Dbgj2c14rlpXe34u8AzgEPbd8VB5wHHVtWlMx/mvPVpYHG6U4bem+Q3Rh3QPPBXdM+7A4C/G1zRRpjuRzdCqPdMctfuMLo3L9rtYVW1km7Y4huB+4GLk+zX6oz9krwP3a8kl7X7z2Zh9uKuy/jhyv89SV2vczWxB7Uh3RCfQWPDlR8N7JdknzV3seCtqw3Hy2wENV8keU+6c5Yvn2D1OePadtlsxzeP7NJ+CPgCcEFVfWLUAc0zX6D7UXlnYEVV/YxulOjmdAnGYC/ZZ4GDk2w7gjjnjKq6ClhC913nwgmqjB+ufPfAugOAW+iG1i9oU2jHrwKPBfacYN1n6HrFN5pgnYCq+jHda/hoYDVwTpJXjjSoOa6qfgKcA5w11gMObNY+Y24HtgYuGlV8s8kkdwJJHkU3LO9fkqwAjqP7VTNVdU9VfaKqjgP+hq7HFx74JflZdAnudXQfAAvyfNzmGro3p/Vx97hzHrcGbpu+kOadDWnDB2kfEpcAz5mOgOahodtwwNPoXtsL1TXA08cWquoYul+FF40son4Y+7HlaVV14qiDmW+q6lvAVsCL6D5/oRu98Srgu+09cMzZwCnAhUkeMauBzj3nA3/P1Icqk+R/A1sAL6QbLfQrMxTbfDJZO34TeBldcvbkceuObbfvncHY5r2quq+qLqmqE+ja7HdGHdM8cH/7G3N3+7H5scBD8ZzcBe0ldOPYH1tVS6pqMfBd4H+14cckeQjd0JQb2jZfpBuqvKiqbm0ngK8GDmbh9uR+FnhYkj8YK0jyTLoX2dr8F/CKVnczug+Hz81kkHPchrThgyTZGNgLGH9exkIxHW2YJK+lOz/ok9Mf4rzxWWDTJH80UDZrX3LbbJqvnqVjvT/JW9r9fZOs3MD9rGvba+gmpFnf/S5JN7vrxmtZ/6Yk/7K++x1Gkmcn+Va6WaMPWfcW0+oy4HU8kOReBryeCT5/q+qdwMXARxbyRHLA6cCbq2pKw+PbZ/LbgWPaNucBfzGD8c0Xk7ZjVX2R7lSYC5I8ZmDV/XQ9wE9I8uaZD3P+SfKEJLsOFO3BA9+7tZ6q6k7gtcAbkqz35858Y5I7scPoZmMd9J/A+4GPtem2rwLuBd4N0CagWk33hWXMZcC2wNdnON45qSX6Lwaen+7SLdcAJwKrWPOc3Ne2zV4H/HYbVvGl/5+9+w6rsnwDOP592EMZIk5QVNyouLdl2tBKrZyZaTZsmJrtvX4t26Xt0rJcWWplw1Jz5gBFcIsLEBBk73We3x/vAUFBkXUY9+e6uDjnPe+4XzgW93me576BH7XWmy0Rf3VwmZ/h5eSvyQ3GWOP3c6UFWo1VwM9wH3AUY839UK11dqUFW82Zf5ZjgKuUUjlKKROwBKPIBQDm99xC4MJRMh9gt1KqyLQ9ZbTZyDMnRvlftb5tmDk5zcJY55xT0efXWr+uta6SDwQKeQWYr7Wup7VefeGLymjbs0oplaaUOq2Uur2kE5k/WHpLGe2B4pRS85RSCnBSRouRDPP7Jl0p9RfGv1FvIMA8nfF7jPW595rfU1cXPr/W+kmM4nOLzR9a1znm4pkflvDyhWtyBwDPA6u11gfN+7yEsR61bQnnqBMu83PM3+c3jLWSf5pnC+Zvz8IYDBmllKoTo2tXqB7wrTLaqwVjzJB8ybIh1WzmKt77OF9Ut9ZShSqZCyGEEKViXsqRhZHUfGze1gWjHUY7oJXW+pR5uwJCATeMImgPFTrPNIy2GpedTq+U+hf4Xmtd6SOUSqlFQITW+jlzgvS91vqKWyWVdKxSqhvwpda6uLV6pTmvD8YMI9tiCtpYhFIqFLhfa11sP1+l1FKMD9fvxhiRWQsM0FofKGbfGcBcjOnwGmMN2Uda68/MH5R0wfgwOgV4ARintc7vaT+NUr6nhBBC1E518tNLIYQQFWIxcGeh51OB4npeD8boMT4bY+SnPFNE2yildimlkpRSa5RSDaD4acFKqVNKqeHmxy8ppVYopb5TSqUopQ4opXoV2re7UmqP+bXlnG+9dRGlVDOl1E9KqVil1MlCM1FQRiu5RUqpBKXUQYwZABcefz/G+r3nzFOOfc3bFymjmNdacxw7lVJtLvPzmK6UilRKRSmlHi10jZeUUt+bH+dPbZ6qlApTSp1TSj1baN8+SqkApVSyUuqsUuq9S9z7vUqpUKVUvFLqF3V+Cc9xjJHTX80jp/YXHOeMsZbuea11qrmi7C/AlBIuNRV41zxKdgZjmuw0AK31Lq3111rreK11DkbLoPaFR8iEEELUbZLkCiGEKKsdgItSqqMyKoROwJgmeqGpGL37lpuf31SOa94JTMdImnOBj67g2FEYhYfcMBKs+QDmpHs1RtLeAKMFTbHFTcxTW3/FmO7VHGOkcY5S6nrzLi8Cbcxf12PcexFa68+01p201uuKucQkjGmN7hij369d5p6GAm2B64Cn8pP6EgzC6BAwDHhBKdXRvP1D4EOttYs57hXFHayUugajdcx4jPXppzF3IdBatwHCgJvN05WzLji8HZBnbsuXbx/GdO3idKboUp9L7TsEiNZaxxXa1t2czB9VSj2vSli7LIQQonaSJFcIIUR55I/mXotRSfRM4ReVUX11HLDEPOq2kosTv35KqcRCX/0udT2t9X5zm4TngfGq9C04tmqtf9da55nj7pZ/fYziTx9orXO01iuB4toigTEy66m1fkVrna21PgF8yfn1TeOB18yjjOFcWRIO8LN5pDIX+AFjWu+lvKy1TjMXvVmIkSRfat8MrfU+jKQx//5zAF+lVEPzKOuOEo6fDHyjtd5jTmKfBvqbp05fTj0g6YJtSVy8druk/ZOAeuap7wWUUl7AAoypzfk2Y/THboTxYcUkjC4JQggh6ghJcoUQQpTHYozG89MofqryLRgjrvk9JH8ARiilCrcd2qG1div0VVKSBUaxoHynMZLThqWMNbrQ43SMStE2GKPCZ3TRIhUlVfBsCTQrnJQDzwCNza83KybGK3FhjPUus/+F12pWhnPfjTHSelgptVsZbWKK04xC92NuzROHMaJ9OamAywXbXDDW1JZmfxcgtfDvyPweWgd8orUuaN+itT6htT6ptTaZk/9XMLomCCGEqCMkyRVCCFFmWuvTGAWQRlJ8Be+pGMlUmFIqGmMqsC2XHnG8FO9Cj1tgjEKeA9Io1M7IPLpb2v69UUDzC0YJW5SwbzhG79XCSXl9rfXIQue6MMbKdOG1SlM1vAit9TGt9SSMkc+3gJXmNbQXiqRQ6y3zPh5cMHpfgqOAzQWVeLtRtCNBYQc4P9J80b5KKXeMBPcXrfXlpnRrQF1mHyGEELWIJLlCCCHK627gGvMU4gJKqfw1qzdhTLv1x0hW3qKYtaqldIdSqpN5GvQrwErz9OOjGCOzNyqj/99zgP2lTlTIfxijzbOUUjZKqVuBkqoe7wKSlVJPmotMWSul/JTRexmM9axPK6XczVNpHy7jfZbW80opJ6VUZ+Auzq97LjWl1B1KKU+ttQlING/OK2bXJcBdSil/c2Gp14Gd+VW0L8X83vgZeEUp5ayUGojROmVxCYd8B8xVSjU3F7d6FKONH0opF+AvYJvW+qli7meEUqqx+XEHjGntay4XoxBCiNpDklwhhBDlorU+rrUOKOalKUCQ1nqd1jo6/wtjnWpXpZRfGS63GCPZicaogDzLHEMS8CDwFcbIYhoQUfwpLoo/G7gVY8p1AkYBrWL7SpsT6psxEvaTGKPIXwGu5l1expjSexJjpLGkJK6ibMIoULUeeKeEYlaXcwNwQCmVilGEaqLWOvPCnbTW6zESxp8wRqzbcGW9Fh8EHIEYjOrSD+S3D1JKDTZfP9/nGAW+QoD9GO2GPje/dgvG2ui7VNH+yvmj5sOAYKVUGsY0+Z8xEnIhhBB1hPTJFUIIIYQQQghRa8hIrhBCCCGEEEKIWkOSXCGEEEIIIYQQtYYkuUIIIYQQQgghag1JcoUQQgghhBBC1Bo2lg7gUho2bKh9fHwsHYYQQm9E/QoAACAASURBVAghhBBCiEoQGBh4Tmtd2t72pVKtk1wfHx8CAorrSiGEEEIIIYQQoqZTSp2u6HPKdGUhhBBCCCGEELWGJLlCCCGEEEIIIWoNSXKFEEIIIYQQQtQa1XpNbrUXvR9yM8Grl6UjEUIIIYQQQtRAOTk5REREkJmZaelQKpWDgwNeXl7Y2tpW+rUkyS0rkwlWTAEHV7h3Iyhl6YiEEEIIIYQQNUxERAT169fHx8cHVUtzCq01cXFxRERE0KpVq0q/nkxXLisrKxj0CETuhdD1lo5GCCGEEEIIUQNlZmbi4eFRaxNcAKUUHh4eVTZaLUlueXSdCK7esHkeaG3paIQQQgghhBA1UG1OcPNV5T1KklseNnYwcDaE74STmy0djRBCCCGEEELUeZdNcpVS3yilYpRS+wtta6CU+lspdcz83d28XSmlPlJKhSqlgpVSPQodM9W8/zGl1NTKuR0L6D4F6jWBzW9bOhIhhBBCCCGEKJNVq1ahlOLw4cMAmEwmZs2ahZ+fH126dKF3796cPHmSDz/8kDlz5hQcN2PGDIYPH17w/OOPP2bWrFlVHn9hpRnJXQTccMG2p4D1Wuu2wHrzc4ARQFvz133Ap2AkxcCLQF+gD/BifmJc49k6wMBZcGoLhO2wdDRCCCGEEEIIccWWLl3KoEGDWLZsGQDLly8nMjKS4OBgQkJCWLVqFW5ubgwYMIDt27cXHBcUFERSUhJ5eXkAbN++nYEDB1rkHvJdtrqy1nqzUsrngs2jgavNj78F/gWeNG//TmutgR1KKTelVFPzvn9rreMBlFJ/YyTOS8t9B9VBz7tgy3uwaR5M+dnS0QghhBBCCCFqoJd/PcDByOQKPWenZi68eHPnS+6TmprKtm3b2LhxI6NGjeKll14iKiqKpk2bYmVljIt6eXkB0L17d44ePUpGRgbZ2dk4OTnh6+tLSEgI/v7+bN++nXnz5lXoPVypsq7Jbay1jgIwf29k3t4cCC+0X4R5W0nbL6KUuk8pFaCUCoiNjS1jeFXMzgkGzITj6yEi0NLRCCGEEEIIIUSprV69mhtuuIF27drRoEED9uzZw/jx4/n111/x9/fn0UcfZe/evQDY2Njg7+/P7t272bFjB3379qVfv35s376dyMhItNZ4e3tb9H4quk9ucSWz9CW2X7xR6y+ALwB69epVc0oW974Htn4AW96BSbVjgFoIIYQQQghRdS434lpZli5dWrDOduLEiSxdupS3336bI0eOsGHDBjZs2MCwYcP48ccfGTZsGAMHDmT79u1kZGTQv39/2rZty+uvv46npycDBgywyD0UVtYk96xSqqnWOso8HTnGvD0CKJy2ewGR5u1XX7D93zJeu3qyrw/9HoR/X4foEGjSxdIRCSGEEEIIIcQlxcXFsWHDBvbv349Siry8PJRSzJs3D3t7e0aMGMGIESNo3Lgxq1evZtiwYQwYMIDPP/+czMxMHnroITw9PTl48CCenp4WX48LZZ+u/AuQXyF5KrCm0PY7zVWW+wFJ5unMfwHXKaXczQWnrjNvq136zgB7F6m0LIQQQgghhKgRVq5cyZ133snp06c5deoU4eHhtGrVis2bNxMZGQkYlZaDg4Np2bIlAAMGDGDHjh3ExsbSqFEjlFJ4enqyZs2amjGSq5RaijEK21ApFYFRJflNYIVS6m4gDBhn3v13YCQQCqQDdwForeOVUq8Cu837vZJfhKpWcXSDPvfBlnch5jA06mDpiIQQQgghhBCiREuXLuWpp54qsu22225j2rRpNGjQgKysLAD69OnDzJkzAXB3d8fT05POnc9Pr+7fvz/btm2jW7duVRd8CZRRCLl66tWrlw4ICLB0GFcmLQ4+6AIdboTbvrR0NEIIIYQQQohq7NChQ3Ts2NHSYVSJ4u5VKRWote5Vkdcp63RlURJnD+g9HfavhLjjlo5GCCGEEEIIIeoUSXIrQ/+HwdrO6J0rhBBCCCGEEKLKSJJbGeo3hp7TIHgZJJy2dDRCCCGEEEIIUWdIkltZBswCZQXbPrB0JEIIIYQQQghRZ0iSW1lcm4P/ZNj7PSRHWjoaIYQQQgghhKgTJMmtTIMeAW2CbR9aOhIhhBBCCCGEqBMkya1M7i2h60QIXASpMZaORgghhBBCCCEuopRiypQpBc9zc3Px9PTkpptuAmDRokV4enri7+9f8HXw4EECAgLw8/MjOzsbgOPHj9O6dWuSk5Mtch/5JMmtbIPnQl42bP/Y0pEIIYQQQgghxEWcnZ3Zv38/GRkZAPz99980b968yD4TJkwgKCio4KtTp0706tWLIUOG8M477wDw0EMP8dprr+Hi4lLl91CYjUWvXhd4tAG/22D31zBwjtFHVwhxEa01z67eD8BrY/xQSlk4IiGEEEKIKvbHUxAdUrHnbNIFRrx52d1GjBjB2rVrGTt2LEuXLmXSpEls2bLlsse9/vrr9OjRAxsbG3Jycpg0aVJFRF0uMpJbFQY/BjlpsPNTS0ciRLX1zbZTLNkZxpKdYfx14KylwxFCCCGEqFMmTpzIsmXLyMzMJDg4mL59+xZ5ffny5UWmK+eP+rq5ufHkk0/y9NNP88knn1gi9IvISG5VaNQBOo6CnZ9D/5ng6GbpiISoVvaGJfDG74cY3rEREQkZvPTLAQa1bUg9e/lPlBBCCCHqkFKMuFaWrl27curUKZYuXcrIkSMven3ChAnMnz+/2GP/+OMPGjduzMGDB2nfvn1lh3pZMpJbVYY8DlnJsOsLS0ciRLWSmJ7NzCV7aeLqwLvj/Hntli6cTcnkvXVHLR2aEEIIIUSdMmrUKB577LErmnL822+/kZSUxF9//cXjjz9Oenp6JUZYOpLkVpWmXaHdCNjxCWSlWDoaIaoFk0nz6Ip9xKRksuD2Hrg62dKzpTuT+rRg0faT7D+TZOkQhRBCCCHqjOnTp/PCCy/QpUuXUu2fkZHBo48+yoIFC+jSpQujR4/mtddeq+QoL0+S3Ko05HHISDCKUAkh+HLLCdYfjuGZkR3p5n1+Gv+T13eggbMdz6wKIc+kLRihEEIIIUTd4eXlxezZs4t97cI1udu3b+fVV19lzJgxdOrUCYCXXnqJZcuWcezYsaoM+yJK6+r7B2SvXr10QECApcOoWItvMSqmzQ4GOydLRyOExQScimfCFzu4rlNjPpnc46JqymuCzjB7WRAvj+rM1AE+lglSCCGEEKKSHTp0iI4dO1o6jCpR3L0qpQK11r0q8joyklvVhjwBabEQuMjSkVRLOXkmEtKyLR2GqGTxacY6XC93R94a27XYdkGjujVjkG9D3v7rCGeTMy0QpRBCCCGEqIkkya1qLfuDz2DY9iHkyB/uhWmtuefbAIa8vZETsamWDkdUEpNJ88jyIOLTsllwew9cHGyL3U8pxf/G+JGdZ+KVXw9WcZRCCCGEEKKmkiTXEoY8BqnREPS9pSOpVhbvOM2mo7Fk5ZiYsTiQtKxcS4ckKsGnm46z6Wgsz9/cCb/mrpfc16ehMzOH+rI2JIqNR2KqKEIhhBBCCFGTSZJrCa2uAq8+sPUDyJWpuQChMam8/vshrm7vycK7enM8NpUnVgZTndeMiyu340Qc7647ws3dmnFH3xalOmbGVa1p7enM86v3k5GdV8kRCiGEEEKImk6SXEtQCq56ApLCIXiZpaOxuJw8E3NXBOFoa82827oy0LchT9zQgbUhUXy15aSlwxMVJDYli1lL9+Lj4cwbt3Ypdh1ucextrHltTBciEjL4aINlK/UJIYQQQojqT5JcS/EdDk39Yct7kFe3p+V+vCGU4Igk3ri1C41cHACYMaQ1I/ya8MYfh9h+/JyFIxTllWfSzFm+l6SMHBZM7kE9e5srOr5/Gw9u6+HFl5tPcCRa+kwLIYQQQoiSSZJrKUoZfXMTTsL+nywdjcXsCUtgwcZQbuvhxQ1+TQu2K6V4e1w3WnvW4+Ele4lMzLBglKK8Pt5wjG2hcbwyujMdm7qU6RzP3tiReg42PLsqBJP0zhVCCCGEqDBKKaZMmVLwPDc3F09PT2666SYAFi1axMyZMwEwmUxMnTqV6dOno7XGx8eHc+eMQanAwEBatWrF3r17q/4mCpEk15Laj4RGnWHLO2Cqe2sN07Jymbs8iCYuDrw4qtNFr9ezt+HzKT3JyjXxwA97yMqtez+j2mBb6Dk+XH+MW7s3Z3wv7zKfp4GzHc+M7EjA6QRWBIRXYIRCCCGEEHWbs7Mz+/fvJyPDGFj6+++/ad68+UX7aa25//77ycnJ4auvviqy/Cw4OJixY8eyfPlyunfvXmWxF+fK5gyKimVlZVRaXnkXHFwDfrdaOqIq9drvhzgdn87Se/uV2EamjWc93hnXjfu/D+SlXw7yxq1dqjhKUR4xyZnMXraXNp71+N8tfqVeh1uScT29WBkYwRt/HGZ4p8Y0rGdfQZEKIYQQQljeW7ve4nD84Qo9Z4cGHXiyz5OX3W/EiBGsXbuWsWPHsnTpUiZNmsSWLVuK7DN79mzi4uJYvnw5Vlbnx0sPHTrE1KlTWbx4MX369KnQ+MtCRnItrdNoaNgONr8DJpOlo6kyGw6fZcnOMO4b3Jp+rT0uue8Nfk148Oo2LN0VxvLdYVUUoSiv3DwTDy/dS1pWHp9M7oGTXfk/U1NK8fotfqRn5/L62kMVEKUQQgghhACYOHEiy5YtIzMzk+DgYPr27Vvk9SVLlhAYGMiyZcuwsSn6d93o0aOZP38+gwYNqsqQSyQjuZZmZQ2DH4VVM+DoH9DhRktHVOniUrN4YmUIHZrUZ+517Up1zKPXtSfkTBLPrzlAhyYudPN2q+QoRXl98M8xdp6M551x3WjXuH6Fnde3UX1mDGnD/I2hjO3pxQDfhhV2biGEEEIISyrNiGtl6dq1K6dOnWLp0qWMHDnyotd79OjB4cOH2bVrFwMHDizy2vDhw/nqq6+4/vrrsba2rqqQSyQjudWB31hwbwWb5kEt7wurtebpn0NIzsjh/Qn+2NuU7h+BtZXio4nd8axnzwPfBxKXmlXJkYry2HQ0lgX/hjK+lxdje3pV+PlnXuNLSw8nnlu9n8wcWasthBBCCFERRo0axWOPPcakSZMueq1Dhw6sWLGCCRMmcODAgSKvzZ8/H4AHH3ywSuK8HElyqwNrGxg8F6KCIPQfS0dTqVYGRrDu4Fkeu77dFVfZdXe24/MpPYlLy2bWsr3k5tWd6d01SVRSBo8sD6Jdo/q8PMqvUq7hYGvNq6P9OHEujU//PV4p1xBCCCGEqGumT5/OCy+8QJcuxdfBGTBgAJ999hk33ngjYWHnlxFaWVmxdOlSjhw5wgsvvFBV4ZZIktzqoutEcPWu1aO54fHpvPzrQfq2asDdg1qX6Rx+zV353xg/toXG8fa6IxUcoSivnDwTDy/ZS2ZOHgsm98DRrvKmqwxp58nN3Zrx6b/HOR6bWmnXEUIIIYSoK7y8vJg9e/Yl97npppt48cUXueGGG4iLiyvYbm9vz5o1a/jll19YsGBBZYd6SUqXI6FSSj0C3ANoIAS4C2gKLAMaAHuAKVrrbKWUPfAd0BOIAyZorU9d6vy9evXSAQEBZY6vxtn1Jfz+GNz5C7S+ytLRVKg8k2bSFzs4GJXMn3MG4+XuVK7zPbc6hO93hPHJ5B6M7NL08geIKvHGH4f4fNMJPpzoz2j/i8vOV7SYlEyGvbuJLs1d+eGevuWu3iyEEEIIUdUOHTpEx44dLR1GlSjuXpVSgVrrXhV5nTKP5CqlmgOzgF5aaz/AGpgIvAW8r7VuCyQAd5sPuRtI0Fr7Au+b9xOFdZ8C9ZrA5rctHUmF+3LLCXadiuflUZ3LneACvHBTZ7q3cOPxH/dx7GxKBUQoymv9obN8vukEt/dtUSUJLkCj+g48cUMHth+PY9XeM1VyTSGEEEIIUb2Vd7qyDeColLIBnIAo4Bpgpfn1b4Ex5sejzc8xvz5MybBLUbYOMHA2nNoCp/+zdDQV5mBkMu+uO8IIvybc2qNikh87Gys+ndwTRztrZiwOJCUzp0LOK8omIiGduSv20ampCy/c1KlKrz25Twv8vd14be0hEtOzq/TaQgghhBCi+ilzkqu1PgO8A4RhJLdJQCCQqLXONe8WAeRnNc2BcPOxueb9L90gtS7qOQ2cGtaa0dzMnDweWR6Eu5Mdr9/SpUKnkzZxdWD+7T04HZ/Ooyv2YTLVzrXM1V12romZS/aSZ9J8MrkHDrZVWzbeykrx+i1dSMzI4c0/KrZ5uhBCCCFEVSjPEtKaoirvsTzTld0xRmdbAc0AZ2BEMbvm301x2c1Fd6qUuk8pFaCUCoiNjS1reDWXnRMMmAnH10NEoKWjKbd31x3hyNkU5o3tiruzXYWfv19rD54Z2ZF1B8/y6SapsmsJb/15mKDwROaN7YpPQ2eLxNCpmQvTB/qwbHc4u0/FWyQGIYQQQoiycHBwIC4urlYnulpr4uLicHBwqJLr2ZTj2OHASa11LIBS6mdgAOCmlLIxj9Z6AZHm/SMAbyDCPL3ZFbjor1Gt9RfAF2AUnipHfDVX73tg6wfGaO7tyywdTZltP36Or7aeZEq/llzdvlGlXWf6QB+CwhN5d90RujR3ZUg7z0q7lijqz/3RfL31JFP7t7R4AbA5w9uxNjiKZ1eF8NvDg7GzkeLxQgghhKj+vLy8iIiIoLYP8Dk4OODl5VUl1ypPkhsG9FNKOQEZwDAgANgIjMWosDwVWGPe/xfz8//Mr2/QtfnjivKwrw/9HoR/X4eoYGja1dIRXbGkjBweW7GPVh7OPD2yQ6VeSynFW7d14Wh0CrOW7eXXmYPwblD+4lbi0sLi0nl85T66ernyzI2WrwjobG/Dy6P9uPe7AL7aeoIHr/a1dEhCCCGEEJdla2tLq1atLB1GrVKeNbk7MQpI7cFoH2SFMQL7JDBXKRWKseb2a/MhXwMe5u1zgafKEXft13cG2LvAlncsHUmZvPTLAc6mZPHeBH+c7MrzWUrpONnZ8PmUnuSZNA/8EEhmTl6lX7Muy8rN46ElewBYcHsP7G2qdh1uSa7t1JjrOzfmo/XHCI9Pt3Q4QgghhBDCAso1n09r/aLWuoPW2k9rPUVrnaW1PqG17qO19tVaj9NaZ5n3zTQ/9zW/fqJibqGWcnSDPvfBwV8gpmYV0/ktOJJVe8/w8DW++Hu7Vdl1fRo688EEf/afSea51ftr9boGS3t97SFCziTxzrhu1W7U/KVRnbFWiufXyHtACCGEEKIukkVr1Vm/B8HWqUaN5kYnZfLsqv1083bjoaFVP110WMfGzBrWlpWBEfywM6zKr18XrA2O4tv/TnP3oFZc37mJpcO5SFNXR+Ze155/j8Tye0i0pcMRQgghhBBVTJLc6szZA3pPh/0/QVz1rxystebxlfvIys3j/fHdsLUuw9vLZILs8k0znTOsLUPbe/LyrwcIPJ1QrnOJok6eS+PJn4Lx93bjyRsqd611eUzt35LOzVx4+dcDJEsPZSGEEEKIOkWS3Oqu/8NgbQdb3rN0JJe1eMdpthw7x7M3dqK1Z70rP0FyJHx9LXzUHZKjyhyHlZXigwndaerqyIM/BBKbklXmc4nzMnPyeOiHPdhYKxZM7lGtqxfbWFvx+i1diE3N4t2/jlg6HCGEEEIIUYWq71+pwlC/MfScBsHLIOG0paMpUWhMKq//foir2nlyR98WV36CsB3wxdUQcwiykuHne8FU9uJRrk62fHZHT5IycnhoyR5y8kxlPpcwvPLbQQ5GJfPe+G40d3O0dDiX1c3bjTv7teS7HafZF55o6XCEEEIIIUQVkSS3JhgwC5QVbH3f0pEUKyfPxNwVQTjYWvP22K4opa7sBAELYdFNYOcM966HG9+DU1tg01vliqtTMxfevLUru07G88bvNat4V3WzJugMS3aGMeOq1lzTobGlwym1R69vj2c9e55ZFUKufNAhhBBCCFEnSJJbE7g2B//JEPQDJJ2xdDQX+XhDKMERSbxxSxcauTiU/sDcbPh1Dvw2B1pfBfdugEYdwX+Scb+b5sHxjeWKbUz35kwb4MM3206yJqj6/exqgtCYVJ7+OYTePu48dl17S4dzRVwcbHnx5s4ciEzm2/+q70wIIYQQQghRcSTJrSkGPQLaBNs/snQkRewJS2DBxlBu7dGcEV2alv7AlLPw7c0QuBAGzoHbV4Cj+/nXR74Nnu3h5/uMfcvh2Rs70tvHnad+CuFwdHK5zlXXZGQb63AdbK35eFKPshUTs7CRXZpwdXtP3lt3hKikDEuHI4QQQgghKlnN+4u1rnJvCV0nQuCicid9FSU9O5e5y4No4uLAS6M6l/7AiEBj/W3UPhj7DVz7MlhZF93HzhnGLYKsFPj5nnKtz7W1tmLB7T2o72DDjMWBJGVItd3SevGX/RyNSeH9Cf40cb2CUfpqRCnFq6P9yNOal345YOlwhBBCCCFEJZMktyYZPBfysuG/+ZfcLTsvm5j0GI4lHGN39G7Wn17P6tDVrD+9nsCzgZxIOkFiZiImXb41iq+tPcTp+HTeHd8NFwfb0h209wdYOAKsbeCev8HvtpL3bdQRbnwHTm6GzW+XK9ZGLg58ekcPziRkMHd5ECaTLtf56oKVgRGsCIjgoat9uaqdp6XDKRfvBk7MGtaWvw6c5e+D1eNDIiGEEEIIUTlsLB2AKFmOKYekrCSSs5JJzEokKSuJxPZDSD70A4kuziTpbBKzEs+/np1EUlYSGbmlm5JpraxxtXelgUMD3OzdcHdwp4FDA9wd3HG3Nx67ObgVeWxrZSSzGw/H8MPOMGYMaU2/1h6Xv1heDqx7DnZ+Bq2GwNhFRh/gy/GfDCe3wL9vQssBxrFl1LNlA164uRMvrDnAxxtCmT28bZnPVdsdPZvCc6tD6Ne6AXNqyc/p3sGtWb33DC+u2c+ANh4428t//oQQQgghaiOldfUd0erVq5cOCAiwdBjllmfKIyU7pUgimpSVdD5xLSZRTcpKIjUntcRzWqNwdXDH1d4VN3s3XO1ccbV3Pf/cvuhzZ1tnUrNTSchMICErgYTMBOIz4wse5z/Pj0lT/Puivl19XO3cORNnhb2qz/UdfWno2OCiBNndwfhytHGEtHPw4zSjYnL/mTD8ZWMkt7SyUuHLoZCZBPdvhXqNrvA3cJ7WmkdWBLIm+CTvjO9A95bOpOemk5GbQXqO+XtuepHHJb2WmZtJG7c2DG0xlMHNB+Nq71rmuKqTtKxcRi/YRmJ6Dr/PGnRlxcSquYBT8Yz97D/uHdyKZ2/sZOlwhBBCCCHqPKVUoNa6V4WeU5LcstsdvZuY9JiLktQiz7OSSMlOKTFpVChc7F0umagWvObgiqudK26/P4nzic2oOSHg6Fbh95VryiUpK6kgIY7PjCcxM5H4rHgSMhJYd+Q4senx+DTSpOclkZiZSK7OLfZcjlZ2uOdk4Z6Xi5tHBxp4dipIgBs4NCiSELs7uFPftj5KKbTWZJuyycgxJ5pnQ0j/aToZjTuTfvVTpOdlXpSEFpeQXvhaRk5GibEWx1pZ42TjhKONI062xvf8x7ZWtgTHBhOXGYe1sqZn454M9R7K0BZDaV6veUX9OqqU1ppHV+xjVdAZvr+7LwN9G1o6pAr31E/B/BgYwa8zB9GpmYulwxFCCCGEqNMkya1mbvvlNo4mHC14Xt+2fpHEtCB5zX9uV/S5q70r9e3qY6WucGl0VDB8PhiGPgtXPVHBd3VpPwaE8/jKYJ4Z2YH7hrQBjMQoJSelyIhwQmYCCac2k3D0NxLsHIlv1o0EcknITCAxK7HEKdU2VjY4WDuQkZtBni59sSkbZVOQhBZ8t7n4ef7jrGwbvtwUgZtDPZ4b6Y+bQ72CfQsfZ2tle8m+vyZtIuRcCP+G/8vGsI0cTzoOQFv3tgz1Hso13tfQyaPTlfcOtpBlu8J46ucQ5gxvy5zh7SwdTqVITM9m2Lub8G7gxE8PDMDaqmb8bsoqKDyRA5FJTOrdAqtafq9CCCGEqHkkya1mQhNCsbYy1rW62LlgY1WFa/yWTITwHTAnBOzrV8klw+PTGfHhFjo3c2HJvf1KTg7ycuGfF40CWS0HwrhvoV7RwkUZuRlFk+JC06az8rIuTlZtHXGydsTxv09wOr4Rx9Gf4NRyUMHrttalLHxVyKajsUxbuItR3ZrxwQT/CklEw5LD2Bi+kY3hG9kbsxeTNtHIsRFXe1/N0BZD6dOkD3bWduW+TmU4FJXMmAXb6O3TgG+n96nVyd+qvRE8snwfr47xY0q/lpYOp1Ikpmfz1p9HWLY7DK1htH8z3hnXrUa2gRJCCCFE7SVJrjgvIhC+ugaGv2T00K1keSbNpC92cDAqmT/nDMbL3an4HdPjYeVdcOJf6H0v3PAGlCEBLVFWitF+KCvVvD63fFV/5284xjvrjvLizZ24a2CrionRLCEzgS1ntrAxbCPbIreRkZuBk40TA5sPZKj3UIZ4Dak263hTMnMYNX8baVm5/D57MA3r2Vs6pEqlteaOr3cSHJHE+kevolH92rPuWGvNz3vO8Prvh0jMyOGuAT7Ud7Dl/X+OMrS9J59M7omjnfXlTySEEEIIUQUkyRVFLb7FmLo8J9joK1uJPtt0nDf/OMy747pxW0+v4neK3g/LboeUKLjxXehxZ+UEEx0CXw4Dn4Ew+SewKvvIlMmkuW9xIP8eiWHJvf3o06pBBQZ6XlZeFjujdrIxfCP/hv/LuYxzWCtrejTuYazj9R6KV/0Sfq6VTGvNw0v38ntIFEvu7Ve6atm1wInYVG74YAvX+zXh40ndLR1OhQiNSeG51fvZcSKe7i3ceG1Ml4J1x0t2hvHs6hB6tnDn66m9cXWqwA+fhBBCCCHKSJJcUdTp/2DhDXD9G9D/wUq7zMHIZEYv2Mrwjo35ZHKP4qf1HlgNqx8AB1cYvxi8e1daPAAEfAO/PQLDXoDBsRcOJQAAIABJREFUj5brVMmZOYyev42UzFzWzhpE40quJmzSJvaf21+Q8IYmhgLg6+ZrrONtYazjveK12mW0eMdpnl+9n8evb89DQ32r5JrVxQf/HOWDf47x3fQ+DKnBvYAzsvP4eMMxvtxyAic7G54a0YEJvbwvWoO7NjiKOcv30sazHt9N71OrKmcLIYQQomaSJFdcbNFNcO4YzN4HthX/B2tmTh6j528jPj2bv+YMoYHzBetJTXmw8TXY8i549YEJi6F+kwqP4yJaw093G8n1tN+MHrrlcPRsCmMWbKNDk/osu68/djZVt24xPDmcDeEb+Df8X/bE7ClYx3uV91UM9R5Kn6Z9sLeunOnD+88kcesn2xng68E3U3vXucJEWbl5jPhgC7kmzbpHhuBgW/Om8W44fJYX1hwgIiGD23p48fTIDpecbr7lWCwzFgfSsJ49i+/uQ0uPyp0FIoQQQghxKZLkioud2ATfjYKR70Cfeyv89K+tPciXW06y8K7eDG1/QX/ajET46R4I/Rt6TIWRb4NNFa7lzEyGL66CnExjfa5z+abZ/hYcycwle7mzf0teGe1XQUFemcTMRDaf2VziOt7BzQfj5nCFbaP2fg/bPoJbv4Bm/gWbkzNzuOmjreTkmVg7a/DFH2DUEduPn+P2L3cyc6gvj13f3tLhlFpkYgav/HqQPw9E49uoHv8b41fqqeZB4YlMW7gLW2srvpveh45NpZWSEEIIISxDklxxMa3h6+uMdbAP7wGbiktUth8/x+SvdjK5bwv+N6ZL0RdjDhvrbxNPG8ltr+kVdt0rErUPvhoOra6C21eUa30unE/qL7n2uIoUXse7KXwTsRmxWCtrujfqXrCO19vF+9In2fk5/PEEKCtw8oDpf4FHG7TWPPD9Hv4+dJbl9/Wjl0/lrEWuKeYuD+LX4Ej+mD0Y30ZVU628rHLyTCzador3/zmKSWtmDWvLPYNaX/Hsg2NnU5jy9S7SsnP5Zlpvetfx90BtprWuMW3MhBBC1D2S5IriHfsbfhgLoz6usGJPyZk53PD+ZuxtrVk7axBOdoXaIx36DVbNAFsnGP8dtOxfIdcss91fwdpHYfjLMGhOuU6Vm2fijq93sjcskZ8eGIBf8+pR/dikTRw4d6CgPdGF63iHeg+lc8PORdfxbn0f/nkJOtwEQ5+Bb2822k1NX8fC4HRe/vVgkX7Hddm51CyGvbuJ9k3qs/y+fqVKCHLycsrUuqo8Ak8n8OyqEA5HpzCsQyNeGtUZ7wYlVDovhYiEdO78eheRSRl8MrkH13RoXIHRCks7HZfG4yuDSUrP4aupvcr1XhFCCCEqiyS5onhaG211MhNhZiBYl79f79zlQazZF8nK+/vTvYW7sdFkgk1vwaY3oVl3mPADuDYv97XKTWv4cRoc+hXu+h1a9CvX6c6lZnHzx1uxtlL8OnMQ7tVwGm94cnhBwpu/jtfT0dNYx+s1lL7HNmG/+R3wuw1u+dxo43QmEBbdTEb9lgw4+yg927Xkyzt7yQiP2bJdYTz1cwjzxnZlfK/zI+SZuZmcTDpJaGIoxxKPcTzxOKEJoUSmRdLEuQn+nv74N/LH39Ofdg3aYWtV8Ymv0fP2MEt3hdPU1YGXRnXmuk6NK+R3dy41i2kLd3EoKoV3x3VjTPdq8G9alIvWmh92hvH674ewtlJYKYWtteKbab3p6nWFyx2EEEKISiZJrijZod9g+WQjoek2sVynWhscxUNL9jB7WFseubadsTEz2Ri9PfI7dLsdbnq/UgpdlVlmEnw+BPJyjPW5TuWberk3LIEJn++gXxsPFk7rjXU1LsiUmJlo9OMN38jWM1vJyM3A0WRioH0jhvaezRDvqwvW8aYdWIfdj5MIUe1p/cifuLnIWsx82bk53PbVL4SlnmDKEHvOpBmJbVhKGCZtAsDGyobWrq3xdfPFu743p5NPExQbRHRaNAAO1g74NfQrSHq7enbF3cG9zDFprfnJ3PM2KSOHuwe1Yvawtjjbl/+DrMJSMnO497sAdpyIr5Se0aLqnE3O5ImVwWw6Gsvgtg2ZN7YraVm5TFu4m7jUbObf3p1hHWXEXgghRPUhSa4omckEnw0CUw48uAOsylYl9mxyJtd/sJmWDZxY+cAAbK2t4Fyosf42LhRueAP63AfVcfQvcq+xPrn1UJi0rNzrc5fsDOOZVSE1pyCRyUTW2kfYdXAZG3168q/KJDYjFitlRZeGXWjp0pLA4+By9jRz9N80aTGYJrcuwsm+nqUjr1ImbSIyNZLQxFBjdDbhGKGJoZxMOkmOKce8l8LHpSW+br74uvvi6+ZLW7e2eLt4FztSG50WTVBsEPti9hEUE8Th+MPk6lwAfFx86ObZrSDxbe3WulTtoY6eNXre7joZT8+W7vxvjF+lFojKzMlj1tK9rDt4llnD2vLI8LYyyl/D/LIvkudX7ycrN49nR3bkjn4tC36HMSmZ3L0ogAORSbw6xo/JfVtaOFohhBDCIEmuuLT9P8PKu2DsQvC79YoP11ozdeFudp2MY+2swbTxrAdH1xkVlK1tYNy30GpwJQRegXZ+AX88Dte+CgNnletUWmue+imE5QHhfDGlJ9d1roLWSGVlyoNfHoagH2DgbBj+MiY0B+MOsiFsAwFnAzgWF0ZKbjxKFf0372bvRlPnpsZXPeN7E+cmBds8HD2qrGdvRdJacy7jHMcSjxGaEFqQ1IYmhpKRm1GwX1PnpgXJbFu3tmw9aMPy7Vksu3dwqasVXygjN4ODcQcJigkqSH4TshIAqG9bn66NuhqJr3m019n2fBufjOw8PtpwjC83n6Cegw1Pj+jAuJ4X97ytDLl5Jp7+OYQfAyOY0q8lL4/qXOfaStVEienZPLd6P78FR9G9hRvvjutGa8+LP7xKy8rl4aV72XA4hgeubsPj17WX368QQgiLkyRXXJopDz7pB1a2xpTdKxzJ/O6/U7yw5gCvjvFjSt8WsPU9WP8qNOkCE38AtxaVE3dF0hpWTIEjf8Bdf4J373KdLjMnj/Gf/8fJ2DTWzBxY7B+OFpeXAz/fCwdWwdXPwFVPXDTSHng6ngmf7+CaDh68cEtzotOiidr1KdHH1xHl3ZMot6ZEp0UTmRpJem56kWNtrGxo4tSkIAEu/NWkXhOaODXBydayBW2SspKMBDbBWDebn8wmZSUV7NPAoQFt3doWjMz6uvnSxq0N9e2KVlPOyM7jug82YWdtxe+zB2NvU/7euVprwlLCCpLeoJggjiceR6OxUla0dWuLfyN/bLJbsWanLVFxTozr6c1TIzrgcYmet5VBa82bfxzm880nuLlbM94d161K+0aLK7PxSAxPrgwmPi2bR65tx4whrbGxLvn3lZtn4oVfDrBkZxij/Zsxb2zXCnmPCyGEEGUlSa64vH3LjLWzE36AjjeV+rDQmFRu+ngLfVt5sGhyR9Sah+DgGvAzV222q0FVOTMSjfW52gQzNpd7fe6ZxAxu/ngrHs52rH5oYIWvhyyXnEyj6NbRP+C6/8GAhwHIM2nOJmcSHp9OREIG76w7go214reHB+PqaJ5uqzX8Ogv2fAcj3oa+96G1JiUnhajUKCPpTYskKi2K6NRootKiiEqLIjYjtmCNar780eDCI8BN6jWhmXOzCh0NTs9JNwo/mYtA5Y/QxmbEFuxTz7beRdOM27i1wcOx9KOyG4/EcNfC3Tx6bTseHtb24h1iDsHax6DvfdBpdJnuJSU7hZDYEIJig9gZuYd9sfvIIxMAF9sG9GpiLmjVyJ9OHp2wt67aZPezTcd584/DXNXOk0/v6FG0wrqwuLSsXF77/RBLdobRvnF93h3frdTV4LXWfLrpOPP+PELfVg34YkovXJ2qtlK4EEIIkU+SXHF5ebkwvxc4uMJ9/5Zq7WxOnonbPt1OWHw666e1wOO3uyD2MFz7CvSfWT3X317OmT3G+ty218LEJeW+h22h55jy9U5G+DVl/u3dLb5W0WTSxCYk4PjTFFwit/Kv75P84XATEYlGUhuZmEFO3vl/2/UdbFhyTz+6eF3wR3BeLvw4FQ6vhbFfG9WYLyPHlENsemxB0hudFk1UalTB86i0KNJy0oocc+FocOFkOP954dHg7LzsgorGhUdoz6SeKdjH3tqe1q6taevetmBktq17Wxo7VUzV4Yd+MPoIr5szBJ+G56cTE7YTlow3ip0pBTd/WObWXTl5JhZuO8n7fx9Dk8fkwQ608YrlQHwIQTFBhKWEAcbPr5NHJ/w9/QvW9zZyalTue7ycZbuMden+3m58M603bk7Vr9J4XRRwKp65K/YRnpDOfYNb88i17XCwvfLR2DVBZ3jsx334eDiz8K7eeLnXoA8zhRBC1BqS5IrS2fOdsT5z8kojybuM9/4+ykfrj7FieAZ9Ah41No79BnyHVXKglWzHp/DnU3D969D/oXKfLn9kqyp6y2qtOZeaTURCOuEJGUQkGMlreHw6ZxIySEiM4zOrefRSR3gy9z5W5l1Fw3r2eLk74t3ACS93R+Oxu/G4ubtjyVMSczJg8a0QsRsmr4A215Q79sKjwYWT3/znMekxJY4GZ+VlcTr5NHk6DwBrZY2Pi0+RkVlfd1+86nlhXcYCa6VxNjmT4e9uwr+FG99N72Mkzkf+NEbOXZoZH56sexZC/4HrXoMBM6/o/AGn4nl21X6OnE1heMfGvDSq00VJRlxGHPti9xWs691/bj/ZpmwAmjk3o1ujbgUtjNq5t8PGquJHW//cH8WspUG0aujMd3f3obFLNaqqXsdk5ebx3t9H+WLzCbzcHXl3nD99WpVvpsp/x+O4b3EADrbWLJzWu9r0BhdCCFF3SJIrSic3Gz7uAfWbwt3rLjmKuScsgXGfbed9ry2Miv0cPDvCxO+hQesqDLiSaA3L74Cjf8H0v8CrZzlPp3nwhz38dSCa7+/uywDfhuU6V2J6DuGFkteIhIyC5xEJ6WTmFE0CGzjb4eXuSDuXHOaefYYm6Uc51P8d7P3H0tzNCUe7ciR8GYmw6EaIPwnTfoPmPcp+rlLINeUWGQ2OSosiKjWKyLRIbKxsjETWPOXYx8UHO2vLjCDmr1P/cKI/o/W/xodHTbsaHyA5NzT+rf18LxxcDUMeh6HPXnbWQEJaNm/+cZjlAeE0y+95W8qiZjl5ORyKP0RQTJCR/MYEEZMRA4CjjSNdGnYpGOnt5tkNV/uKSVi2h57j3u8CcHe24/u7+xYd2RZV4mBkMnNXBHE4OoVJfVrw7I0dqVdBSyeOnk1h2je7SMzIYcHkHgxtX/mzBIQQQoh81S7JVUq5AV8BfoAGpgNHgOWAD3AKGK+1TlDG/MEPgZFAOjBNa73nUueXJLccdn0Jvz8Gd/4Cra8qdpf07Fxu+eAfHs1cwHWmzcbawtGfQG1qKZORAJ8NMR7fvxkcy96zFCA1K5fR87eSkJ7Dbw8PopmbY4n7JmXkFCSvEYWS1/B443tadl6R/V0cbAqNwjrhnf+9gRPN3R2NP2hTY2HxLXDuiFHtusPIct1PESnR8PW1kJ0G09dBQ9+KO3cNlWfS3LpgK8MSljHL9L3RnmrCYrAvVKzKlAe/zoa9i432Wje8VWzRN5NJszIwgjf+OERKZi53D2rFrHL2vNVaF7Qvyi9qdST+SMEoeCvXVvh7+tO9UXf6N+tPE+eyVwgPjkhk2sLdWCn4dnofOjeTEb+qkJtn4vPNJ/jgn6O4Odkx77auDO1Q8Uno2eRM7lq4myNnU3htjB8T+9SAQoNCCCFqheqY5H4LbNFaf6WUsgOcgGeAeK31m0qppwB3rfWTSqmRwMMYSW5f4EOtdd9LnV+S3HLIyYQPu0HDtsbIXDHeXv43Iw48Rmer06hhz8OguTVz/e3lRATAN9dDuxtgwvflvsfQmFTGLNhGG09n3ri1K2cSiyav+SOyKZm5RY5ztrM2J7FOF00r9nJ3Ol8QqiTJUfDdKEgMN6pdV8Z08rjjxlpmWydjFoBL04q/Rk1iMnFu1RM0DPmSfW7D6TZzKdgUM6qsNax7Dv6bD10nwugFRtstsyPRKTy3OoTdpxLo7ePO/8Z0oX2T+hefpwKk56RzIO5AwUhvUGxQQZVpHxcfBjQbQP9m/endpHeR1kWlERqTyp1f7yQlM5evpvaibxlbLInSOXUujbkrgtgTlsiNXZryvzF+uDtX3qyG1KxcHvxhD5uPxvLwNb7MvbadxesPCCGEqP2qVZKrlHIB9gGtdaGTKKWOAFdrraOUUk2Bf7XW7ZVSn5sfL71wv5KuIUluOf33Cfz1tNFKp2X/Ii/t3fQrLTY8gLONxmHCQmh3nYWCrCLb5xvrJ294C/rdX+7T/bk/mvu/DyyyzdHWuoQ1sU54N3DE1dG27H8wJpw2Ety0c3D7CvAZWO57KFHkXlh0E7i1hLt+B0e3yrtWdZabDWsegpAV7Gw0nolho1j5wCB6tixhNoDWsOUd2PA/aH8jjP2GdG3Dh+uP8fWWk9R3sOHpER0Z29OrSnuTaq05lniM/yL/47+o/wiMDiQzLxMbZUNXz64FSW9nj86lWuMcmZjBlK93EpGQwYLbezC8U+MquIu6RWvN9zvDeH3tIWytFa+O8WNUt2ZVknDm5Jl4btV+lgeEc2v35rx5W1dpISWEEKJSVbck1x/4AjgIdAMCgdnAGa21W6H9ErTW7kqp34A3tdZbzdvXA09qrQMuOO99wH0ALVq06Hn69OkyxSeA7HT4oIuxhnDKKmOb1qRt/QT79c9xxqo5TWf8hF3j9paNsypoDUsnGUWC7v4LmpdvfS4YBVvi0rIKRmY9nO0q54/QuOPw7SjIToE7fgavCv1vQPGOb4QfxhnXmrIKbEuell0rZaXCijvh+HoY9gJpvWcx/P3NuDra8uvDg7C9RB/S/KUCcY36MTFpFseSYEIvb54c0YEGlTgKV1pZeVkExQTxX+R/bI/czqH4QwDUt6tPv6b96Ne0HwOaDcCrvleJ54hPy2bawl0ciExm3m1dua1nyfuKKxOdlMkTPwWz+Wgsg9s25O2x3WjiWrXFvrTWzN8Qyrt/H2Wgrwef3tETFwdpMSSEEKJyVLcktxewAxiotd6plPoQSAYeLiHJXQu8cUGS+4TWOrC484OM5FaIre/DPy/BPRugcWf0b3NQ+5byj6kn3ncvpn3L5paOsOqkxxv9c5WV0T+3JoxQxhyC70aDKRemrDY+sKgq+3+GldOh/QgYv7jI9NtaLS0OlowzRrQLtQf660A0MxYH8vSIDsy4quTq2hEJ6fy55EOmxcwj1MaXjPHL6d6++hZyi8+MZ2fUzoKk92z6WQC863sbo7xN+9O7aW9c7FyKHJealct93wWw/Xgcz9/UibsHtbJE+LWG1ppf9kXy/Or95ORpnrmxI3f0bWHR6cI/BUbw5E/BtPGsx8K7el+yBoEQQghRVtUtyW0C7NBa+5ifDwaeAnyR6crVR1aKMZrb2M8oKBS5h/dzbsPp2qeZcXVbS0dX9cJ3wcIR0H4kjP+ueq9BjgwyikxZ28HUX8DTAiPu+QXMut8Bo+ZX759XRUgMM9opJYUbbbQ63Fjk5Xu+DWBb6Dn+njvkonY/OXkmvtpyko/WHwPg/W5nuP7Q0yiPNsZoeP2yF32qKlprTiafNKY2R/7H7ujdpOemY6Ws6NKwC/2b9WdAswH4NfTD1sqWrNw8Zi8N4s8D0cwc6suj19XuNZwmbSIqLYpTSaeITovGzd6NRk6NaOTUCA9HjzK3cEpIy+a5NftZGxxFjxZuvDven1bVpIL1ttBz3L84EGd7G76Z1ptOzVwuf1Atdjw2lff+PkrgqQTem9CNAW3KXmVfCCGEoVoluQBKqS3APVrrI0qpl4D8/yvHFSo81UBr/YRS6kZgJucLT32kte5zqfNLkltBNs2Dja9hsnVmdvaDnG02jKX39sO6CtcFVivbPoK/n4cRb0Pf+ywdTfHCd8H3Y8HBBe5cAx6V25f3kja8BpvnGYXJhr9ouTgq29mD8P1txodBty+DlgMu2uVMYgbXvreJ/q09+Gpqr4KEbtfJeJ5bHcLRs6lc16kxL47qTHM3RzixyZgmX8/T+D26+1TxTZVPTl4O+2L38V+UkfQeiDuASZtwtnWmd5PeDGg2gL5N+vH5P8ksD4jg9r4teHW0X43/b0taThqnkk5xMvmk8T3pJKeSTxGWHEZmXmaxx1gpKzwcPPB08qSRUyMaOzXG09GzIAnO/3KxcynyQcDGwzE88VMwienZzBnejhlDWmNzqenwFnA4Oplp3+wmNSuXT+/oweC2npYOqcqdSczgo3+OsXJPBPY2VjRytiYiKYf/SSVqIcotz6Q5eS4N30a1qLuHuCL/b++8w6Mqvgb83vSEhIQOgUDovXdBRAFBBelFARFRVPSzK/afvWMvWJAiUsUuHRERCB0SOiEkEEgnvW2b74/ZhA2kbJJtCfPm2Wfv3jL3ZHb23nvmNFdUcrshSwh5AVHATMANWAU0Bc4BE4UQl8wlhD4HRiBLCM28Mh73SpSSayPyMjBte48nT3diU1It1j16PSG1/co+rrpiMsHyKRC1FWZtguBuzpaoKGe3w7LJENBAloAKCnGuPELAn4/B/kUw4h3o96Bz5bEHMbtg+WTw8IXpP0ODjiXu+u2/Uby59jjzp/WgT/M6vL32OKv3x9I4yJdXb+94dSKm2P3w43jw8JEu5/Xb2fmfsR/p+enStdms9F7IugBAcI1g/EztiTjdgJuaDeCzKQPw9ihn3ebkSFj3NDTsDEP+B1YkwaoMBVbZs+lSkY3OiC5cLqg9DFJ5beLfhNDAUEJrhtI8sDmhNUMJ9g8mLT+NpJwkEnISSMpNIjEnUS7nyOW0/LSrzuvj7kM9v3rU8alLYqoP0Qke1PWtx919u9KjcbNCZdjb3duu/395iUvPZebCvUQmZvH2uM5M7OXk65KDSM7K54utkfwYdg6A+3oF8UjeV3id3cK8mnP5PLYlswe1YO6IdlV+ckehcAZ5eiOPLD/IxmMJPD28LQ/dqMoXXou4nJJrb5SSazvmbzvDO+tO8MHErkxQSWJkfO78gdIV+P5t4OMiNT9Pb4aVU6XF767fXMfF1WSE1TPg+B8w7jvoMtHZEtmOk+tg9d0Q2EQm9qrVrNTdDUYToz7fQVJmHgaTICvPwL3Xt+CRIa3w8yrBXTXhmHQ9N+pg2k82SXzmbIQQnM88XxjLuyd+D1n6LITQ8KMZkzoO4YaQAXSr1w1P91KSFplMsOdrmTsADQy50P52GPeNTRKeFVhlo9KjiM6ILrTQnss4R74xv3C/AK+AQgW2eWBzmtdsTmhgKCEBIXi5VyxhWL4xv1DhTcxNJDE7kaTcJI4nnWd/bDQ60vD0ysSI7qpjA70DpcLrW/8qa3A9v3o08GtAbZ/auGmlWH2FgL3fQWo03DBXeoZUgsw8PQ8uPcB/kck8NrQ1jw5pXW3d09Nz9Xy3PYoF/50l32BiYs8mPNXyPHU3PwE5KRAUgkiN4bcmT/HY6a4M69CAjyd3q1Tda4XiWiM738DsH/axIzKFriFBHD6fxjMj2jJnsFJ0rzWUkqvAaBJk5ulJz5WvtJzLy4WvHD1puTrzZwPpOTriM/K4uUNDvprWo9o+lJSbc2Gw8FZoPwomLnJ+vOnxP6WyVb+dtPjVcLFYL30e/DgBzu2CO1dCq6HOlqjyHFwKvz8CjbrC1NVW9/nBc6lMnL+LHk1r8cbYTrRpYEXN20tRsGSMfEC+YwU0v76SwrsWBpOBI8lH+HbferbG/Ie773nQTPh6+NKrQa/CUkUtAltcvgalRsOvD0HMf7KO9ahP4Mga2PAChPSFO5aDX+0yz200GWWsrIU1tmA5KTepcL8Cq2yBMhsaeNkyW9untt2vjfkGIx9uPMU326MIqeXHvEld6dWsFhm6jKLKcM7Vr5S8FEzCVKQ9D82DOr51pGu03xWu0d61qb9vEfUP/4S/EFCzMYz6FFpX7nerM5h47ucI1hyIZWLPJrw1rnPp2carGLk6I4t2RjN/2xnSc/WM7NKIJwY3ocWBd2DfAqjXXk7A1G4Oq2bAmS0cbPEA449fT7uGgSy4uxeNAlWCLoWiLNJydNy9cC8RF9J5b3wXRncL5olVh/n98EXmjmjHg4OdGKalcDhKya0mCCHI1hlJy9EVUUwLFdcrFFa5Tkd6jp7MfAOlfWXeHm4E+XkS6Fvw8iLQ15OGgd7Mvr4lgX6qDEQRCrJP3zYPet/rPDkifoKfZ0Nwd2np8y2hFquzyUuHhbdJhW3GH9CkilokhZDf/ZZXoeVNMnu0d/ligdJydOWvfZwRBz+MgUtnYdJimbm6GrLhaDz/t3InDepd4PqulwhP2UtMhiwHV9+vPv0b9ad/no5++5ZRxwTc8g50m3p5ounoL/Dz/RDUVP4ezLHMWbqsQuW1IE42OiOamPQYdKbL1tDirLLNA5sTEhBSulXZjhy9mM4TKw9zMiGTO/s25YVb25fL6mcwGUjJTSEpV7pHJ+YkXnaVLlCQcxLJ1GdedWyAhx/N8/MIzcmged1OhPacRWjdTjSt2bRCVmohBB9tPs2nW05zfeu6fDm1BwFVvMSQzmBixd5zfPZ3JEmZ+dzYth5PDW9LR+Mp+OV+ec3r/xDc9BJ4mks6GfXw+//B4eVcbDmJWyLH4u3lxXczetGlSRXI3q+oEFtPJnImMYuZA5orF/UKkpiRx/QFezibnM1nd3ZneEfptWYwmnh81WH+OHyRZ29pxwOlVDJQVC+Ukuti5OmNRaypRZTW3KutrRlmBTYjV4/BVHK/e7hpUkG1UFaDfC0UVz+vy+v9Lm+v6euJj6d949iqHSYTLJsEZ/+FezdJi56jObBEWhObDZAJj7ytsAo6k8wE+P5myMuAezZAvTbOlqh8mEyw8QUI+xI6TYAxX4GHA+vX5lySCa7iDsPY+dBlkuPFP1zUAAAgAElEQVTO7UB2nUnhviX7CPT15IdZffD2TZeuzTFb2H1xJxlIq2S7wBb0D7mB/o3606NBDzw0Dy5mXyT69Dqi/3ufs14eRAd3IjovpYhV1l1zp0lAE2mRLYiVNcfNOsIqay0Go4mv/43i482nqOXnxbvju3Bju/p2O19O0kmSVk8jMTOWxN4zSazfmgtZF4hOi+Js8hESjbmF+7ppbjT2b1yk/wrctOv41CmzD1ftPc9zv0TQpkEAC+/u7fB6vrbAaBL8evACH20+RWxqLn1Ca/P0iLb0buIvE+5tnyet4GO+Kt77Qgj4+3XYPo+sZkMZHT+LCzkaH03qxi2dGzn+H1LYjex8A2/8dZzle2R89pB29fn0ju7KRb2cnL+Uw7QFu0nKzOfbu3oxoFVRDypLRbeskn2K6oNScl2MYR9u43RiVrHbNA1q+ngWUUZrXqGsFiiocr1XoVJbw8vdZR7Qrgmyk2V8rqcvzN5W6bi1crH7a1j3jLQmTv4RvKpIQrBLUbDgZplM6Z4NEFhF6i0bdPDbHIhYDX0fhOFvgZsTXC3zM2XW5ej/4Nb3oc99jpfBARy5kM6M7/cAsHhmbzolr4N1czGa9Bwf+BA7g+qxKy6MQ0mHMJgMeLt7I4QoYpWtaYLmej2hTa4jNOS6QsusM62y1nI2OZsnVh3i4Lk0buvSiDdGd6JWDTtOqJzfCyvukLHfk5dC80FX7ZJ9bifRax/nbOY5oht3JbpBG85mXyAmI6ZojLJnQJGJgwLL+JXW322nkpizdD81fT1ZOLM37RpWjRJDQgg2HI1n3sZTnE7MolPjmjw9vB2DWtdFSzoJv8yWE1Fd75SeBmXlbdj7Hax9Gn2DbtxreJptsYKnh7dlzuCW6n5eDTh4LpXHVx4i5lIOs69vQcNAH17/8xjtGtZULurl4HRCJtMW7CZPb2LhzN70aFq815rBaOKxlYf4MzyO529tx+xBStGt7igl18X4+UAsOTrjFRZVaWUN8PHATbmxVB1idsKi26DjWBi/wDHxuQWu0m1vg4kLwcO1sqmWSdxh6boc2ARmrrUqdtKp5GfBqrvgzBaZvXfg486Nw9bnyRjsU+tgyMuyRFM1fBiOSsrike828njeFwzR9kHT/jD6iyJlsXL0OexL2EdYXBjumnuReNlaOh3a8kkQfwRGfgg973beP2MlQgiWhsXw1toTeHm48fqYTtzeNdi+Jz36C/zygExWd+fq0j0sDDppodz+AfjWhpEfYmp3G/HZ8YVu4AWxzWczzpKYUzTbdIH1t+A7QleP9/5IJS/Pl6+n9+K6Vi6WT8ACIQT/RSbz/oaThMem07JeDZ68uS23dGqIJgTsni+vy97+Mka8/SjrGz/+J6yZhSkgmDeCXuf74xoTejbhrbGd8fKoPnHL1xIGo4nPt0by2d+RNKzpw7xJXenXog4A/5xM5OFlB/HzcmfBjN50buIiCSxdlPDYNGZ8vwcPdzd+mNWnzAkxg9HEoysP8Vd4HC/c2p77BrVwkKQKZ6CUXIXCnvz7gXQ7G/kx9Jppv/MIAf+8DdvehU7jYezX4OIWqRI5+690vQ3uLpNluaolOjsZfpwIcYdk8p0e050tkcSoh1/nQMQqGPAoDH21+im6R3/B9McTGPIymWecQs/Jz3Fzp3Ja/vOzZHbvyM0w6Bm48XmX7ae49Fye+Smc7aeTGdSmHu+N72JfN14h4L8PYctrMlnXlGXWJ62Lj5DjLz5cTvDd+kGxx2brsy9npi6IhTYn97K0/momX4z5dekR3IZBoR0L3Z8rGvtra/bHpPL+hhOERV2icZAvjw1tzdjujWVd4rTz0svj7L/mJGifyjJu5eXcblg+GeHmwfJWH/D8bk/6NK/N/Gk9qW1PK77C5pxNzubxlYc4dD6Nsd0b8+rojtS8Ivb8ZHwm9yzay6VsHR9P6VYYW+oyxOyEOq3A334hEtYQFpXCvYv3EeTnydJZfQmtW8Oq4ywV3Rdva8+91ytFt7qilFyFwp6YTLKeafQOuG+LrNdpa4SAjS/Crs+h+zT5IGXneqB259hvMstom+HSRdLVFPa0c7J8T3osTFgI7W51tkRFMZmky/reb6WV8rYPq/6YABl7vPYpmS05uAfpIz5jxh/phMem8e74LuWvs2rUw5+Pw8EfpAvp7Z+61FgTQvDboYu8/NsR9EbBC7e1Z2rfpvZ1VTXoZJ8cWirjy0d/cTkpkrUY9bDjY9j2nswHcOv70HGcVZMIJmEiPju+0OJ7MuUMG08fIdN4ATfPjML93DQ3gmsEX477tUgKZk3sb2U5HpfBvI0n2Xw8kbr+Xjx8Yyvu6NtU1nIWAsJXwtqnQZhkCEOPuyo3iZJ0Sk7+5aSwq9eHzNgeSKNAHxbM6E2r+uVLcKdwPEIIVuw9z2t/HMPTXePNsZ0ZVYonRlJmPvct2cfh2DSeHdGO2YNaON9FPSNOXn9P/Ak16sH476DFYKeIsuV4AnN+PEBIbT+Wzupb7kk/g9HEoysO8VeEUnSrM0rJVSjsTVaSjM/19ofZ/9g2CZTJJG86+xZAn9kw4l3nxIPag33fy4ftblPlg7azb/AFJByVD5v6HLhjJTTr72yJikcI+PsN6T7acZy07jsyGZatObke/nhEKrqD58KAx8Hdg+x8Aw8s3c/208kVcz8TQipj/7wFLW6ESUscG0NfApeydbz4awRrI+Lp2awW8yZ2tdpSUWFyU2HldIjeLmvgDn6ucr+7xOPSqnvxALQbKSdbKmDJ1BlMzF0Tzi+Hori1hzu3dHfnXFZMoRU4JiOGPGNe4f4BngFF4n671utKt/rdbGL5PZuczUebTvFH+EUCvD24/4aWzBwQermedXYK/PkYHP8dQvrJRHC1m1f6vABkxsuSawnHiB7wDhPCmqMzmPhqWs+rEu0ozJiMTp/gS87K59k1EWw+nsB1Leswb1JXq+Jt8/RGnlx9mL/C45jSO4TXx3RyTmktkwkOLoGNL4MxH657RE5Ep5yGwc/D9U869Lnjt0MXeHLVYdo3qsnie/pU2JtBbzTx6IqDrI2IV4puNUUpuQqFI4j+DxaPkq7E4761jcJmMspSE4d+rL5uqf+8K5WPAY/CsNecLQ3E7ILlk8HTD6atgQYdnS1R2ez4FDa9BK2GSQXOVd2/SyIvHdY/Ly2LDTrJjLSNuhTZJd9g5ImVh/krIo4HB7fkmeFty2/1OPijVKLrtZf1jWs6L4vt3ycSmLsmgrQcHU8Ma8vsQS3sX1bkUhT8OEnWGR79OXSdYpt2jQYI+wL+flOOvRHvyuzf5fx+hBDM23iKz7dGMrhtPb64s0dhBtorrb+W7s8JOQkA+Lj70LNBT/o16ke/4H60qdUGN836B/O49Fw+3XKaVfti8XJ3Y+aAUO4fdEUJvVMb4feH5UTMTS9IZcDWClZeBqyaDlH/kNZ/LpOPDiAyOZvXR3fizr5NbXuuqowhX4YL7fwUOk+QOROc4F7794kEnvkpnIw8A88Mb8s9A5qXK7eKyST4aPMpPvs7kuta1uGrqT0dW7Yx5Yys1BDzH4ReL2PK67SU4R5/PiYTLrYaJus8OyCHxtKwGF767Qh9Qmvz3YxelS4zpjeaeGT5QdYdieelkR2YNdBGE1IKl0ApuQqFo9j2Hmx9U7oT95xRubaMevj5PpkYZvDzcMMz1U/BBWllW/u0dLu9+U247mHnyXJiLfw0EwJDYPrPst5qVWH/IvjjMZmg6c4VZWd1dRXObIXfHobMizKJ1g1zS7RGG02Cl347wrLd55jSO4Q3x3Yuv2IYuUUmEvMJkrV067e3wT9hPVn5Bt748xgr9p6nXcMAPpzUjQ7BDrAqnwuTmbkRMiN76ADbnyPpFPz2EMTukfGpIz+CmuVPnLVs9zle+u0I7RsF8P2M3tSvWbqbYqYuk/0J+wmLC2PXxV1EpUcBUNunNn0b9qVfcD/6N+pPI//iJzVSsvL56p8zLAmLQQjB1L7NmHNjS+oHWJw3P0uGjOxfCPU7wriv7ROaUoBBJ/syYhW67ndzf8oUtp66xL0Dm/Pcre1VndXY/bJ/ko5LxexcmKx0MPg5mXXeASEJOTpZGmjZ7nO0axjAJ1O607Zhxb241uyP5dmfwwmp7cfCu3vTrI6dvTqMetj5Gfzzjqx4cPPrV7vcCyE9rtY/CzXqyzrtTWyqTxThy38ieW/9SYa0q88XU3vYrLylpaL78sgO3KMU3WqDUnIVCkdhMsLScfKGe9/fFbcCWmbQHfY6DHjEpmK6HCYj/HQPHPtVutzaysJUHg4sgT8elcmw7lwNNeo4XobKcuRn+Hm2VNym/2J9IiFnkJ8Fm/8nS6jUbQNj5kOTnmUeZmntu6VTQz6e0k3GSJaHuHCZUEyfC1N+LL6OqRUIIcjTm8jM05ORZyAzT09mnsH80he+Z1isO3IhnfiMPO6/oSWPDW1dftkrQvhqmRwpMERasOvYsayGySgzDW95Hdy9YMRbMhyhnBN0W08k8tCyA9Ty82LRzN60bmC98pCQncDu+N3suriLsLgwknOTAWhWs5m08jbqR++GvXETfny7/SwLtkeRqzcyrkcTHh3SmpDaV3hCnNsNv9wvLeDX/R/c9KJjstqbTLDlVdjxMaY2t/K231N8GxbP0Pb1+WTKNVpnVZ8LW9+S+SkCGsmEj21uhuTTUhGL3Az12sEt79o1lvTw+TQeW3mI6JRs7ru+BU/e3MYmv+XdUSncv3Q/GvD19F70aW4ny+nFg9JLLD4C2t8uY+oDSkl+dfGgnBzMiIPhb8rQKRtOugsheHf9SeZvO8PobsF8MLGrzd229UYT/7fsIOuPxvO/UR2YOUAputUBpeQqFI4kK1HG5/oEwn1bZZxuedDlwIo7IWqrzFpaTWuhXoUhXyoeMTvgjhXQephjzmuZZbblEOnuW97vzJU4vUnGXAaFSEU3sImzJbqamJ3w64OQGgP9H5JKg2f56kV+tz2KN/46zoBWdfh6ei/8y/nAL1JjEEsnoKVFk3DTx8SF3FqmgpqZZyDDYltmngGDqfR7oaaBv7cHNX1kibi6/t48NrQ1vUIdUDrLMha52QCZ4M1RJbtSzsiH6Jgd8nc16hM5JstBRGw6MxftRWcw8s1dvQpLsJQHIQRn0s6wK04qvHvj95JryEXDDZHfhPzMlvSq34eXht5C+0ZXtG/QyYz2Oz6Gmk1k7K09LOBlsftrWDcXmvRmZev3eX7DRdo0CGDBjF4EB11DdVZjdkqvj0tnZLK9Ya8V9VgRAk6tl8puarQs43Tzm1Crmc1EMBhNfPnPGT7ZcpoGAd58MKkr17W07WRidHI29yzaS2xqLu9O6MzY7ja8huty5PVg1xfSMnvbB9aXu8pNhV8elJPvHcbA7Z/ZJLeBpYfO1L5NeX10J7uV0tQbTTy87AAbjibwyqgO3K0U3SqPUnIVCkcTtQ2WjIYuk+WDkbUznnkZsGwynA+D2z+H7lPtK6erkZ8p6w4nn4a7foeQ3vY9n8kEG56H3V9B54kw+suqnbipgJidchz5BMJdv9nXclce9LkyUdauL+SD55ivoNl1FW5uzf5YnlkTTqfgmjw+rA3Z+UazEnrZonqlUmq57GfK5FuvD+nrdoK39HfwjXEkcPm3eqWCKl/FLXtS84p1BcfU8HJS7XNDvlQyw1dC1zukkunomtomk7TUb34FNDfpDtnz7nJZgM5fyuHuhXs4fymX9yd2YXS3cpaRskBvNLFsTxSf7fibDI5Ru24MuW5nMQnj1fG8OgNuvz4gLV3dpsGIt52brOzYb7DmPghqStiAb7nvt0R8vNz59q5edAsJcp5cjiA/S1q093wDQc1khvQWg0veX58nLb3b58nM1wMekzkfKpmrICYlm8dWHuLguTRGdwvmtdGdCPS1j1t0Wo6OB5ceYFdUCo/c1IrHh7WpfOblM1tljG1qtPwdDn0VfMs5dkwmGQO95TWZbG3SkkrlrdAbTTyx6jB/HL7InMEteboiuRbKic4gFd2Nx5SiWx1QSq5C4Qy2vg3b3pFZg7tPK3v/nEsyq2bcYZngodN4+8voimQlwfc3y1njezZAvbb2OY9BJ62JR36CfnPkjH91yVoNcPGQzBCtadKia8/4QWuI3Q+/PgDJp6DXLGmFsYHFfPOxBB5adoB8g6nIemsU1Jo+ngR6mRh09CWaXFhHQru7SBv0OgF+3s5VUCtLziVYMRXO7YQbX4RBTzk3nj81WircZ/+F5jdIC1A5rGvpOXru+2Efe85eYu6IdjxwQ/lKrRhNgj8OX+TDTac4dymHXs1q8fTwtvRtUYcsXRZ74/cSFhdGWFzY5Xheo4m+ehP92k+mf4/7S4zndSgxO2H5FPDwIWbEYqb+lUtSZj4fTurGbV1cQD57cGarTBaXdh76PgBDXgIvK2NV0y/AppflNT4wRLrZtr+9QgnRVu49z2t/HsPdTeONMZ0qNdliLTqDiRd/jWDVvlhGdQ3m/QldKhajmnMJNr4kE/vVbiknCUIHVk646B0yf0VeBoz8ELrdWe4mcnVG5vy4n60nk5g7oh0PDnbcZKzOYOKhZQfYdCyBV2/vyIzrQu1/0qht8v7n30C62gc0kMuOnnysZiglV6FwBiajtObG7oPZW0tPcJOVJGuyJp+EiYtdryaro0mNhgU3g5sHzNpoe5fb/Ezp0hu1FYa+Imf6q2NSr+TTcgzqsmSccdO+jpfBoINt78J/H8mYr9GfQ8ubbHqK2NQc4tPziiix5VJQTSbY/LJMwtJupMyOXtUyVBeQHAnLJsoH/DFfyqyzroAQMjnaxpekdW3Yq3Kyw8qJpXyDkadWh/PH4YtM69eUV0Z1xKOMmD0hBJuOJTBv4ylOJmTSvlFNnhnelsFt6xWvJKedI+HX2YQlhxPWoAVhXu4k510Cro7nDfR2UmK3xBNy8iovjYzbv2fm9gD2x6Ty1M1teOjGVs6vs2orctNkoq+DP0CdVnKyuGm/irUVvUPWFE84IidZbnnX6oRzKVn5PPtzBJuOJdC/hSwN5EgXcSEE87dF8e76E/RoGsQ3d/Wirr+VSpEQMs/F2mcgJ0Vas294ptyhISWSmQBrZslyZN2ny7heK9vOzNMza/E+9kZf4o0xnZja13Yu5dZiqei+Nrojd/UPtc+JkiOlx9jpDcVv960t740BDcG/4eVly8/+Dcpfy/waQSm5CoWzyEyA+QPkRWz21uJnoDPiYMntcqZ6yo/Qaojj5XRF4iNg4a1yxvOe9baLJcxONlvMw+WMtjVW9qpM2jlYMgYy4+T4srGCWSrxEfDLA/LhsttUGP5W+d3jHElh7GMvWR+5qiUfi/4PVk6TrsFTljtnUqMs0s7LBG9ntkCzgfI3aKU7vckkeG+DTE4zpF19Pruz++XatVewMzKZ9zac5ND5NJrXrcETw9pwW+dGxU98CAGHl8vvXphgxDvQfRoCio3nddPc6FinI/0a9aN/cH+61utqk/q8VpNxEZZOgOST6EZ+xtzTHfjl4AXGdW/M2+M7OyaZmT05uV661WYlyBJNg5+tvGJmNMjM2H+/ISc5+8yW7ZZyPdp6IpGnfwonI1fP08PbMmtg+UoD2ZJ1EXE8vuoQdf29+f7u3rQpKxFb+gVY+xScXAuNuknviSvKstkEo0HG+G6fBw06y+zLZfyeU7LymbFwDyfiMvlwcjdu71r+DOy2wlLRfX10R6bbUtHNS5c5EXbPBw9fOcHQeSJkJ8l62Fnx8r3gVfA5KwFMhqvb861lVnrNlmBLi3Dh54a2m8SoIiglV6FwJme2SitttzulZcWS1Bip4GYnw52rnJPUxJWJ3iH7rlEXGVtqrZtaSaTGyOzX6bEwcRG0vcUmYro8WYnwwzhIOgETFkCH0fY9n9EAOz6SNZB9a0lFpqr09fE/YM29ULOxLDFUu4WzJbKOQ8ulS3Dt5nDnSteWWwhZ+3v982DUwZCXoe/9Vteb/WFXNP/7/SidGgeyYEZv6gVctmwdPJfKBxtPsiMyhUaBPjw2tDXjezQp2eqbnSyV7hN/QtPrYOxXUCu02F31Rj3hyeHStfliGBHJERiviOftH9yf1rVal6s+b4XIS5cu6dHbEUNe4bP82/hw82l6h9bi6+m9qF2jCuYWyE6RSaMiVskyTaM/h8Y9bHuOnEvw9+uwbyH41YGh/5Mx1xYeBbk6I2+uPcbSMFka6KPJ3WjfyInx2GbCY9OYtXgfeTojn0/twQ1t6l29k8kE+7+HTa9IRemmF6Dvg+Bu50zcpzbIzP7CJK3uHW4vdre49Fymfbeb2NRcvprWg5vaNbCvXFagM5iY8+MBNh9P4PUxnZjer5JWZZMRDi6Vccs5KXIifcjL1tdwNpnkcVnx0lCSGVeMQpwg3036q4/3CbxaCS7OQlxVvZWuQCm5CoWz+ftN+Pc9mWinIHYl5Qwsvh10mTDtZ7vWnqvSHP8TVk2XGVrvWF7x+ofxR6SbnyFXTihU1PWtqpKbBssmQexeOatvLwt20klpvb14ADqOg9vmOS6jr604t1vGPmpucqxYUdrIaQgha3P/+z40HyQTwfjWcrZU1pFxUdZ2Pr0BQvrKh+O6ra06dPOxBP5v+UHqBnixaGYfDEbBBxtPsulYAnVqeDHnxlZM7du09BjGk+vlxEBemszw3f9hqxVtkPV598Xvuzqe11yft39wf/o16me/eF5DvjmvwBroM5s/gh/hqZ+O0KCmD9/f3YtW9Stes9XhHP0F/npKfheDnpY1s+2ZBDDusLTcn9sly8bd8j6E9CY8No3HVhwiKjmbewc256nhbW1Wq9UWXEzL5Z5FezmdmMUrt3csqpAlnZLxy+d2SbfsUR87drIr7RysmiGv/f0flqFAFvfrs8nZTPtuN+m5ehbM6EXfCmRLtxdS0d3P5uOJvDGmE9MqqujG7IL1c+X4CukHt7wjx5c9MJlk7pIrleAss2KcmXDZQmzUXX28d6DZClzUTVpXoy5/JaXSp9ONNG5gp5woNkQpuQqFszEZpUJ78YAsK4SQsZImA0z/1T5uRNWJ/YvlzbvLFDlRUN4EUTE7YdkUOXM57Wdo0ME+cro6umzpznrmbxj+NvSfY7u2TUYI+1LWR/WqIZXbTuNs176jSY6UVv+sRJi40DUt0fo8Wf/2yBo5aXHbR1UvO7gQEL5Kxkzqc+HG5+UDshWWp8Pn05i1eC85OiO5eiP+Xh7MHtSCmQObl15SKj9TxsgdWAINOsna3A07Vfpfic+OZ3fc7kKlt7j6vF3rdcXDzQOBQAhR7DuASZiKrhNg3qtwGwKEyYjY+Sni8HJMLW7gRMfHeXXdGfRGI8/e0pZuTYNAWLRnbhMoXC5sDy7LgsDTzZManjWo4VEDX09fanjWwMfdx7Zxv5kJsPZJ6UHRqJuc6LDBd2EVQkDET7DpJciM43iDkcw8fxtaQAPmTezKda1cs854Vr6BR5Yf5O8TidwzoDkvjGiJ+85P5US6p59MsFWB2tQ2wZAvY6n3fCMnriYshMDGHI/LYPqCPRhNJpbc05fOTZwU014K+QYjc5YeYMuJRN4cW8444bTzsu77kTXSC2jYazJ5qCvEyAthVobjzQrxZSVYZFwkNvsiEXlJhIscIjzdOe7thV7TuNPjOp6b+rWzpS8TpeQqFK5ARpysn+tbS7qiuHtJF9z67ZwtWdXg3/dlPFX/h+VN3FpO/AWrZ0JQU5j+s3y/ljHkS3fc47/DDXNh8HOVvxFfioJf50gLQttbYeTHcoa4qpOVKK3fcYdlzeres5wt0WWyk2H5HRC7p3okT8tMgL+ekG7DwT1kaIcVyYHOpeQwd004XUICefCGlgT5laHkx+yCX+6XVqcBj0ql2g7ZTYUQRKZFFiq8BfG8VR03za2I0lvDowY1PIt+9vP0w8/Tr3Cbn6effPfwu/zZww+/42vx3fQyWjknN2zN+bgk9v7wPLdl/4zJ3QsGzcV34ByXnjAymgRv/HWMAzs387n/QkL0Z6HjWBjxrmtce4+sgd8fAQ9vTg38mAkbvfHz8mDpvX1c2sMg32DkwaUH+PtEIm+N7cydfct4XtDlyJJK/30MCHO5qkcqH1plJ7J0WRxJOUJ4UjjhSeFEJEdwyZxczw0vTDkNqKWvz22NmjJ9wASCg13fAKOUXIXCVYjcIl1mA5u4Vv3SqoAQ0r1sz9dylnTAo2Ufs3+xTGAS3F1mF65qiYTshdEAfz4q44b6PiCtuhUpn2Qywb4FskyHm6fMWtp1StVWtq5Elw0/3QOn1ksXyiEvO///SzoJP06UM/Jjv4aOY5wrj60QAo7+DGufltbWG56RD40VDVGwxJAPW9+CHZ/Iia6xX0Oz/pVv10oK4nlPXDoBgIaGpmkU/mnFfDbXbHbT3Ip8Llh209zkOs3c3vndaLu/wc2/AfkDnuLb/ZmEx6YzrH1DJvduirubZtV5C+TTm/Tk6HPINmTLd3022fpscgw5lz+btxV8zjHI93xjvlX9ognw8/SlhmfAZWXYrCAXqxwX89nXwxchBHqTHqMwYjAZMJgMGIURvUkvl01GDMJQuK3gtTc6mb+OxOKmmbittQehaf9hTDmDwa82huaDMAQ2lvsKcxsFx1q0ZXnOgvWF+15xTgBPd0+83LzwcvfC080TL3e57OXmVWRbke1XrhPgdWojHlH/kGvyY4v3LYwcMZqGNf2tasu9HG75FSb5NNk/3IFvWiSLvSYz9P4PCKnrugpuAVYpugXXqo0vQ0asnGAY9ppLTaIbTUYi0yKJSI4oVGjPpJ0p9NoIrRlKm6COxCc2YM8JfzwMwcwa2IrZN7Sgpo996j/bA6XkKhSuRMwuGSfjCrOtVQ2TCX6+V84SW8Y3X4kQMtvj36/LWN7JP7jszKrTEAI2vABhX0DXO2WcbnmsKGnn4feHIeofmbH59s8h0P61I52C0SBdKvk/MuIAACAASURBVPcvgs6TpEuls6w8Uf/Ayrvk+e9YUT1j+bOTZWbYo79Awy7SqluZOs8JR2VSnIQjstTJiLfB2/UftivE2e0yIZWXH8Y7V/PWfncW/HeWm9rV59M7upfuxm1DDCZDEWVYvmeRfeJPcg4vJ0cTZLceRnajTuQYcguV44J9LT/nGnIdbgX3EOAhTLi7eeDhWQMPd6kYerp54q654+HmgYebB+6aXOfh5oG7mzsemvndzaPIvgXvAHqTHr1Rj86kQ2fUoTPp5Gfzss6oK3YfQ3EZdyuIu+ZerBJtqRgHeAUQ5B1EkHcQtXxqEegdSC3vWoXrgnyCCPQOxNOteIVo/ZF45i7fxTy/xQzVb5X3iXHfQg3XdAO3JN9g5IEfZA3ft8d15o4+Fspr3GFY96ysQ96ws7Seu0DS0OTc5CIW2iPJR8gx5AAQ6B1I57qd6VKvC13qdqGpfzuW70ri+x1nMRgFd/Rpyv8NaUX9gKpXpkgpuQqFovpg0EkX0rP/wpRl0HZE0e0mk8zQuedr5yskro4Q0g1865uyPuyE78t23SzMjPucjMMd/gb0nOl866a9sZw4aT4IJi+VWSwdyYEl8OfjUKe1zKBcy/G1JR3Ksd+lC3NuKlz/FFz/ZPl+yyYj7PpCfmc+gTDq02ujBnnCUVliSJcFk5eyNDGU//1+lNb1/fluRi+a1HJCVtVLZ2VehbP/yt/PqE9lJnArMZqMhYpvEUuyWQl219xLVTItFdK90Wm8tz6SzFwTcwa3Zka/lnh5FFVINaNO5hjY9r7MnTHgURj4uFMz0pqyk9Gvfw7dkdXo6rREP+xVdI26oDPqOJ2Uxit/HCYjP5cHBofSNcS/RAXaKqXaqCNTn0laXhqp+amlTjIEeAYQ5BNELW+zIuxTi8Q0d/49kUPjmnV59IauBCfso9bOLwnyCiJw3Hd4uoBSWBZ5eiMPLN3PPyeTeGdcZ6Z08JXXkgNLZDLFIS/LSTNHWMWvIN+Yz/GU41KpTQ4nIimCi9kXAfDQPGhbu+1lpbZeF5oGNEXTNPL0RpbsiuaLrWdIz9Vze9dgnry5Dc3qVF0jgFJyFQpF9SI/ExaPgsQT0u27oB6oQQe/PiAtvf0egpvfqJgb7rXG7q9l4p8Wg2Hyj+DtX/x+mfGy3Mqp9dBsgJxAKMeDarXg8Ar47SGo2xamrnaM9dpkgi2vwo6PpTVk4iLHK9jOIueSDFOIWCWTRI3+AoK7lX1caozMPByzQ07gjPqkSliQbEZ6rFR0UyJh7Hy2+9zAnB8P4O3hzrd39aR7Uwdl4DYZZRKiLa+B5i4nxXrMcMqkWK7OyNvrjrNkVwxtGvjz8eTudAguozRQRpwMx4hYBTWbSPk7jHGs/ELIe9q6uTL79MDH5aSPZ1GrW1JmPvct2cfh2DSeHdGO2YNa2CxJWL4xn7S8NNLypdKblp9WqACn56eTmmdel59GbHoy6bo0NLdiMvqaCXDzItCvPrV8ilqGLRXlK63IJVmM7Ume3sicJWE0j/qRub6/4WXKgz73y1AKB9V8F0JwPvM84clmK21SBCdSTxRa9hvVaFREoW1fuz0+HkXHhsFoYs2BWD7efJq49DxuaFOPZ0a0pWNw1b+PKCVXoVBUP7KT4fvh8v2e9TLOeeU06c459FU5817drYu25NByqbw17iHL5liW/Sl4yFr7lMyAO+R/Mpb3Wp1AOLMVVk6XLq/TfoIGHe13Ll2OTJR0/HdpMb/1A6ck53E6J9fJckPZSTDwMZk0rTivgwJPg3XPys+3vCvDGq7Fa0FuqnRdjtkBN79BZKu7uWfRPuIz8pg3sSujugbb9/xJp2RIw/nd0PpmGPmRvE47gSMX0nl0xUHOJGUza2Bzni5vaaCYXbDuaYiPgNDr5biy5+++gLTz8NeTssxW454yrKSU8+bpjTy5+jB/hccxuVcIr4/phJeHY67TQgg+2XKajzefZnjHBrw/qQO5hsxC5Tc1P5W0zDjSDv9I2qVIUmuFkF6/Dan6rEIFusC9tjj8Pf2vdp32kYqwv6c/AV4B+Hv64+91eTnAK4AanjUKXcXLzamNmNY/h9ulSLYau5J702vceuPgirVlJRm6DI4kHbms1CZHkJafBoCvhy+d6nYq4npcz6+YeslmhBBsOBrP+xtOciYpm24hQcwd0Y7+LatPfhKl5CoUiupJaoxUdNFkofX4CHMN2KnOlqxqcvxP+GmmdIed/ouMG89OkS6jx36Fxr1g7Hyra5lWa+KPyORPZpdQWtxg+3NkJsCKO+DCAemV0P+ha1NZKyA3VcaRH/oR6rWD0V8WrWGclSQ9DU7+Bc0Gylje6u7SXRb6PDlJcuxX6DeHSwP/x/1LD7A3OpXHh7bhkSGtbFsSCGQM+67PYOvb4OkrFcIuk50ydo0mwfxtZ/ho0ynq+nvzwcSuDGxdQYu+yQgHFkurdF4G9L4XbnzOPnWpTSbY+5304BAmuOkl6Hu/Va6xJpPgo82n+OzvSPq3qMP8aT0J9LOvFdRkErz+1zEW7ohmfI8mvDu+Mx7uJSjXQsCuz2HT/+Tvc+LiwjKKlhbjAsU4PS/9svW4GAtyaYpxAb4evgR4BuDvZVaCPQOkImxe9vfyL6oo56QRsHcB/jE7CQgMxXPI6zwYVo9/Tyfx7rguTOodYpN+M5gMRKZFFomlLai3raHRIrAFXep1oXO9znSp24WWQS2tVth3nknm3fUnOXw+jVb1/Xl6eFtu7tDA9r93J6OUXIVCUX1JOAoLb5GuyhMXXR2jqygfUf/A8jvlpMHAx2UMUm6aLPFx3SPXphWxJCxdQkd/AV0n267thGOwbDLkJMtkLe1H2q7tqs7pTVKZzYyTZWdufF5a1/94BPLSZaxcv4euXU+DKzGZZF3g3V9Bx7Hkj/qC534/zc8HLjCmWzDvjO9SPqtmacQfkR4hcYeg/Si4dZ7Tkiyev5TD4ysPsS8mldu6NOLNMZ3KLjFlDTmXZKbufQukgmvr2MzEE/D7/8nyYC1vkhbwWqHlbmbN/lie/TmckNp+fD+jN6F17RN3aTCaePbnCH7aH8vd14Xy8sgOuLlZoUjF7JKTqrmp0kOlx/QKnV9v1JOpzyRLl3X5XZdJpi6TLH1WkfVZermtyLI+y6ps4L4evhj03uTrvAgJqk1orTqlK8rFWJQv5V0iIimCw8mHiUiK4GjK0cJ451retaRCa7bSdqrbiQCv8ifIO3IhnXfXn2D76WQaBfrw+NA2jOvRuORJhyqOUnIVCkX1JuWMnGWv18bZklQPzu+FHyfI+K+GnWHMfGjYydlSuSa5adJNPnq7fNgd+ETlLVaRm2VtZ08/uHOFLIGlKEpeOmx8SVrWatSTbswNOsO4b6BBB2dL53oIATs/g00vQbOBiClL+TIshfc3nKRns1p8Pb0ndf0rUS/YoJOJ2bZ/AD5BcNs8p5W2EkKw5sAFXvn9KBrw2piOjOnW2PYWrPgIGScbswMadYNb34eQPhVvz5AP/30E/34g8yKMeKfSFvDdUSncv3Q/GvD19F70aV67zGPKQ77ByKPLD7H+aDyPDmnNY0Nbl6+fs5JgzSw4uw26TZXKrhOSe+n1eWQe+J6s/z4kMz+dzLbDyeo0lkw3tyKKcnpeJtsiz5OSm07Tum54eerKpSgX4OHmQfva7YsotU38m1RqjEYnZ/PBxpP8GR5HkJ8nDw1uxfT+zWw3geWiKCVXoVAoFOUj6RTE/Afdpqns1GVhyJfWq4jV0OseuOX9ilu89y6QdWLrt5cZlJ0Uw1hlOLNVJgZqNRQGP6fGalmEr5YJueq0gmlrWHvOjSdWHaKuvzff392bNg0qUFrpwgH47WFIPCoz2o94x2k1yVOzdTz/SwTrjsTTt3lt5k3qat9s0gX5Cja+BJkXocsUGPYqBDQsXzvn90jrbdIJ6DRB9qF/ybGW5SE6OZt7Fu0lNjWXdyd0Zmx321xTcnQG7v9hP9tPJ/PSyA7MGljBJIQmI/zzjsz0X78DTFoCdVvZREariN4B6+fKSYum/WXfl5LcLk9v5L4l+/gvMpn3J3RlQk/Znzqjroil2NKqXKAo+3v506VeF9rVboe3eyUmlSxIzMjjky2nWbn3PJ7ubswa2LzK1bqtDC6p5Gqa5g7sAy4IIUZqmtYcWAHUBg4A04UQOk3TvIElQE8gBZgshIgurW2l5CoUCoXCoZhM8Pdr0hLT5haYsKB8tZlNRqms7fpcJumZ8H31reWqcC5R22RCKp+aMPUnDuuCuXfJPvJ0Rj6f2oMb2lipXOlzpXKy81PwbwAjP3ZquMi2U0k8vfowqTk6nry5Lfdd3wJ3a9xmbYEuG7Z/KPvC3Utm3+37YNmTLvmZsOV1mYG6ZmMY+SG0GW5z8dJydDy49AC7olL4v5ta8fjQNta5FJdAeo6emYv2cOh8Gu+M78KkXjaIUY3cDGvuA6MeRn8GHcdWvs3SSDsnr7lHf5GZs29+DTqOs8pyXpKi60jSc/V88+8Zvv8vGr3RJGvd3tSK+jWrXq3byuCqSu4TQC+gplnJXQX8LIRYoWnafOCwEOIrTdPmAF2EEA9omjYFGCuEKDXwSSm5CoVCoXAKe76V5ZiCu8MdK62zxuiy5cPdyb+gz2wY/raKfVbYl/gIGU+uz4U7lnExqCezFu/jZHwGr9zekbv6h5Z+/Lkw6b2QEgk97oJhrzuspMqV5OmNvLPuBIt2RtOmgT8fTe7mvNIol6Jg/fNwap20lo94F1oPLX7f05tk3ev0WOhznwx3sOPEls5g4sVfI1i1L5aRXRrxwcSuFXJlTczM464Fe4hKyubTO7oxolMj2wmZHgur74bYvXKSYNhrtvfO0OXIcmw7PgE0ma39ukfK7Sadpzdy7+J97DiTzAcTujLeQYpuQa3bL/85Q1qOrHX7xLA2dou5dnVcTsnVNK0JsBh4E3gCGAUkAQ2FEAZN0/oDrwghhmuatsG8vEvTNA8gHqgnShFAKbkKhUKhcBon/oKfZsmEO9N+hjotS943Iw6WT4H4cOkm1/d+x8mpuLZJOycV3dSzMO4bsluN4tEVB9l8PJG7rwvlxdva4+6mYTQJDCaB3mjCmJeJ17a38D34HcaAxiTd+D7ZTa5HbxQYjAKDyXR5X5NcV7CsNwmMJlPhvoXL5mPk8QKD0fKz3MdoMsm2LNowms8TnZLN+Uu5zBwQytwR7VwjBvH0Jlj/rJwEaHsrDH8TareQ27KT5baI1bLe9u2fXa71bmeEEHz9bxTvrDtB96ZBfDO9F/UCrHebjU3NYdp3u0nIyOebu3pyfWvbuFQXwaCDzf+DsC+hSW+YsBCCbGApLnAt3/QyZFyQVtthr1Wq7VydkXuX7GXnmRTmTezKuB72U3SLq3X79PC2dGpc9WvdVgZXVHJ/At4GAoCngLuBMCFEK/P2EGCdEKKTpmlHgBFCiFjztjNAXyFEckntKyVXoVAoFE4ldh8smyQfrO5cBSG9r94nPkJmUM5Ng4kL7eKmqFCUSs4lWH6HrGU74m2MfR7gnXXH+Xb72UIFt4Dr3I7wrse3hLglschwM+8ZppCD7V0jPdw0PNw1PN3ccHfX8HBzw9Ndw91Nw9PdDXc3DQ+L5Rre7jxwQ0v7KFyVwaCTGa23vQdGHVz3f1C7JWx8UbopX/8kXP9E8fWe7cz6I3E8tvIQdWp4s3CmdbHYkYmZTPtuDzk6Awtn9qZnM9smsbqKo7/KWG93T5lhviSLuDVcPChrZ58Pg4ZdZFmrZtfZRMxcnZFZi/eyKyqFDyd1tVnMcwHXQq3byuBSSq6maSOBW4UQczRNG4xUcmcCu65QctcKITprmnYUGH6FkttHCJFyRbuzgdkATZs27RkTE1Ox/0yhUCgUCluQcgaWjpelbsYvKFoG6NRGWT7Du6ZMMGWuE6lQOBx9Lqy5F078KRWxoa+x/lgCERfScXdzo4Ypm+ujP6FD3C+k+zUlrNOrXKrbq1AZvayEupk/X62YWq6/ap27m3mb3L+61fEkM17WhA1fIT836S2tt/XbO1Ws8Ng0Zi22Lhb7yIV07vp+D24aLLmnLx2CazpGyORIWD1Dlgoc9DQMfrZ8pZqyEmVd44NLwa+OudzTNNuVezJToOiGRaXw4aRujOne2CbtWta6bVmvBk8Pb8fwjtWv1m1lcDUl921gOmAAfICawC/AcJS7skKhUCiqE9nJ0lp7Yb8sMdLnPtj9jczm2bCzjNutacOYNoWiIpiMshzO3m9ldt8xX0oL46kN8MdjkBV/uSaxp6+zpa2anN8Ll85A54k2V7IqysW0XO5ZtJfTiVm8cntHpvdrdtU+u6NSuHfxPmr6erL03r40d3Tspy5HZpw/tBSa3yAnDMvKdWDQwe75MmOzPgf6PiCTgfnYz7XXUtH9aHI3RneruKJ75EI67204yb+nkq6JWreVwaWU3CKNmC255sRTq4E1FomnwoUQX2qa9hDQ2SLx1DghxKTS2lVKrkKhUChcBl2OrAV5ci2E9JWuoW1vhfHflS8Ds0JhT4SQ2cG3vArNB0FAsLQ+1msPo7+AJj2dLaHCDmTlG3hk+UH+PpHIzAGhvHhbh8Ks1FtPJPLA0v00ruXL0ll9CQ5y4gTHgR9g7VPgW0vG6Tbrf/U+QsDpjbD+OTmh0Hq4jIeu29ohIubqjNyzaC+7z1ZM0Y1OzmbeplP8cfjiNVXrtjJUFSW3BZdLCB0Epgkh8jVN8wF+ALoDl4ApQoio0tpVSq5CoVAoXAqTUVoj9i2QFrFhr7mMNUehKMLhFTJzMpjjRp90StyownEYTYI3/zrO9zvOMqRdfT65ozt/n0jkiZWHaNcogMUz+1DH3wXGQHwErLoLUmNg6CvSvb7AdTfpFGx4TpYiqtMaRrwNrYc5XMQcnYF7Fu1lz9lLViu6iRl5fPr3aVbsuTZr3VYGl1Vy7YVSchUKhULhcggBWQkQ0NDZkigUpXPxkHRLrtfW2ZIoHMgPu6J55Y9jNA7y5XxqDr2b1ea7u3u5lrKVly4TUh3/HdqNhOFvSdfkPd+AZw0YPBd632f70kPlwFpFNyNPz9fbLte6ndInhEduan3N1bqtDErJVSgUCoVCoVAoFKWy7VQSD/94gF6htfhyak98vVzQ40QICPsKNr0EJgOgQc8ZcOOL1tUmdwA5OgMzF+5lb/QlPp7Sndu7Bhduu7LW7aiuwTx5Dde6rQxKyVUoFAqFQqFQKBRlkqsz4uPp5vpZfM/vgf2Loe9saNTV2dJcRY7OwN0L97Iv+hKfTOnOLZ0a8vOBC3y0+RRx6XkMalOPZ1St20qhlFyFQqFQKBQKhUKhcCDZ+dKiu/9cKiG1fIlOyaFrSBBzR7TlupZ1nS1elcceSq6HLRtTKBQKhUKhUCgUiupEDW8PFs7szewf9pGUmc/8aT0Y3rGh61vJr2GUkqtQKBQKhUKhUCgUpVDD24Mf7+3nbDEUVqKqESsUCoVCoVAoFAqFotqglFyFQqFQKBQKhUKhUFQblJKrUCgUCoVCoVAoFIpqg1JyFQqFQqFQKBQKhUJRbVBKrkKhUCgUCoVCoVAoqg0uXSdX07QkIMbZcpRBXSDZ2UJYSVWRtarICVVH1qoiJyhZ7UFVkROqjqxVRU5QstqDqiInKFntQVWRE6qOrFVFTlCy2oO2QogAWzbo0iWEhBD1nC1DWWiats/WxYvtRVWRtarICVVH1qoiJyhZ7UFVkROqjqxVRU5QstqDqiInKFntQVWRE6qOrFVFTlCy2gNN0/bZuk3lrqxQKBQKhUKhUCgUimqDUnIVCoVCoVAoFAqFQlFtUEpu5fnG2QKUg6oia1WRE6qOrFVFTlCy2oOqIidUHVmripygZLUHVUVOULLag6oiJ1QdWauKnKBktQc2l9OlE08pFAqFQqFQKBQKhUJRHpQlV6FQKBQKhUKhUCgU1Qal5FqgadoITdNOapoWqWnas8Vs99Y0baV5+25N00Ittj1nXn9S07Th1rbpaFk1TRumadp+TdMizO83WRzzj7nNQ+ZXfSfLGqppWq6FPPMtjulp/h8iNU37VNM0zYlyTrWQ8ZCmaSZN07qZtzmrTwdpmnZA0zSDpmkTrtg2Q9O00+bXDIv1zujTYuXUNK2bpmm7NE07qmlauKZpky22LdI07axFn3arrJyVkdW8zWghz+8W65ubx8pp89jxcqasmqbdeMVYzdM0bYx5m8371Qo5n9A07Zj5O96iaVozi20OG6eVkdXRY7WSfepq47SkPnXoOLVS1gfM4+6Qpmn/aZrWwWKbw+7/FZVTc817f0myOvTeX0lZHXr/t3ZMaZo2QdM0oWlaL4t1LvWcWpKsjh6rlZDT5cZpKbK64nPq3ZqmJVmc916Lbba5/wsh1Eu6bLsDZ4AWgBdwGOhwxT5zgPnm5SnASvNyB/P+3kBzczvu1rTpBFm7A8Hm5U7ABYtj/gF6uVC/hgJHSmh3D9Af0IB1wC3OkvOKfToDUS7Qp6FAF2AJMMFifW0gyvxey7xcy4l9WpKcbYDW5uVgIA4IMn9eZLmvs/vUvC2rhHZXAVPMy/OBB50t6xVj4RLgZ49+tVLOGy3O/yCXf/sOG6c2kNVhY7UycrroOC1RVkeN03LIWtNi+XZgvXnZYff/Ssrpivf+kmQNxUH3/srKesU+dr3/WzumgADgXyCs4PyOHKc2kNVhY7WScrrcOC1JVkeOU2tlBe4GPi/mWJvd/5Ul9zJ9gEghRJQQQgesAEZfsc9oYLF5+SdgiHkWYTSwQgiRL4Q4C0Sa27OmTYfKKoQ4KIS4aF5/FPDRNM3bBjLZXNaSGtQ0rRHyhrNLyFG/BBjjInLeASyvpCxlUaasQohoIUQ4YLri2OHAJiHEJSFEKrAJGOGsPi1JTiHEKSHEafPyRSARsGfd7Mr0abGYx8ZNyLECcuxUtk9tKesEYJ0QIscGMlVUzq0W5w8DmpiXHTlOKyWrg8dqZfq0WJw8Tq2R1d7j1FpZMyw+1gAKEpo48v5fYTld8d5fSp8Wi5N//9bIau/7v7Vj6nXgPSDPYp0rPqcWK6uDx2pl+rRYnDlOrZTVJZ5TS8Fm93+l5F6mMXDe4nOseV2x+wghDEA6UKeUY61p09GyWjIeOCiEyLdYt9DsNvCSjdwrKitrc03TDmqatk3TtOst9o8to01Hy1nAZK6+eDijT8t7rLP6tEw0TeuDnAk8Y7H6TU26OH5ko5tfZWX10TRtn6ZpYZrZrRI5NtLMY6UibZaEra4rU7h6rNqyX8sr5yzkzGxpx9pjnJZ2vpKwlLUQB4zVysrpyuO02D7F/uMUrJRV07SHNE07g3yAfKSMY+1x/6+MnJa4yr2/NFkdde+3hawF2Pv+X6acmqZ1B0KEEH9aeazTnlNLkdUSe4/VysrpUuPUyj51pefU8eZr+U+apoWUcWy5+1UpuZcp7ku9cqaupH3Ku76yVEZWuVHTOgLvAvdbbJ8qhOgMXG9+Ta+knGXKUcY+cUBTIUR34AlgmaZpNa1ss7zYok/7AjlCiCMW253Vp+U91ll9WnoDcubuB2CmEKLAKvkc0A7ojXRnmVsZIQtOVcy68sjaVAjRC7gT+FjTtJY2aLMkbNWvnYENFqtt3a9Wy6lp2jSgF/B+Gcc6vU+LkbVgvSPGamXldMlxWkaf2nucWi2rEOILIURL8zlfLONYp11TS5BTNuBa9/6SZHXkvb+yssoGHHP/L+sZxA34CHiyHMc6pU/LkLVgH0eM1crI6VLj1Mo+daXn1D+AUCFEF2Azlz0lbTZWlZJ7mVggxOJzE+BiSftomuYBBCJjhUo61po2HS0rmqY1AX4B7hJCFFochBAXzO+ZwDKku4HTZDW71aSYZdqPtI60Me9v6dZmi36tVJ+aucri4MQ+Le+xzurTEjHfLP4CXhRChBWsF0LECUk+sBDn92mBmypCiChkfEt3IBkIMo+VcrdpL1nNTAJ+EULoC1bYoV+tklPTtKHAC8DtFrP1jhynlZXVkWO1UnK64jgtSVYzjhinVstqwQouu8k58v5fGTld8d5frKwOvvdXSlYLHHH/L0vOAGQM6z+apkUD/YDfNZl8yNWeU0uT1ZFjtcJyuuA4LbVPzbjMc6oQIsXimv8t0LOMY8vfr8KGgcZV+QV4IIObm3M5SLrjFfs8RNHEQ6vMyx0pGtAfhQy6LrNNJ8gaZN5/fDFt1jUveyLjsx5wsqz1AHfzcgvgAlDb/Hkv8gdcEHx+q7PkNH92Q/4AW7hCn1rsu4irE0+dRQbz1zIvO61PS5HTC9gCPFbMvo3M7xrwMfCOk/u0FuBtXq4LnMacYAFYTdGEPnOcKavF+jDgRnv2q5W/qe7IB4PWV6x32Di1gawOG6uVlNPlxmlJsjpynJZD1tYWy6OAfeZlh93/KymnK977S5LVYff+yspq/uyQ+395xxQWCYUcOU5tIKvDxmol5XS5cVqSrI4cp+X4TTWyWB4LhJmXbXb/r9Q/Ud1ewK3AKeRN9wXzuteQM8wAPsiHgUhkhi/LgfKC+biTWGT7Kq5NZ8qKdLHJBg5ZvOojEynsB8KRgf6fFPx4nSjreLMsh4EDwCiLNnsBR8xtfg5oTv7+Bxf8QC3WObNPeyMvZtlACnDU4th7zP9DJNK10pl9WqycwDRAf8U47Wbe9jcQYZZ1KeDvzD4FrjPLc9j8PsuizRbmsRJpHjveLvD9hyJvxm5XtGnzfrVCzs1AgsV3/LszxmllZHX0WK2EnK44Tkv7/h02Tq2U9RPkdfwQsBWLBzYceP+vqJy45r2/JFkdeu+3wfc/GAfd/8uS84p9/6GokuNSz6klyerosVoJOV1unJbx/TtsnFr5m3rbov+2Au0sjrXJ/V8zH6RQKBQKhUKhUCgUCkWVR8XkKhQKhUKhUCgUCoWi2qCUXIVCoVAoFAqFQqFQVBuUkqtQKBQKhUKhUCgUdYPrSQAAAFdJREFUimqDUnIVCoVCoVAoFAqFQlFtUEquQqFQKBQKhUKhUCiqDUrJVSgUCoVCoVAoFApFtUEpuQqFQqFQKBQKhUKhqDYoJVehUCgUCoVCoVAoFNWG/wfpkBb7eadMvgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1152x648 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, axs = plt.subplots(3, figsize=(16, 9), squeeze=False)\n",
    "axs[0, 0].bar(range(num_pops), fix_pops)\n",
    "axs[0, 0].set_xlim(0, 11)\n",
    "axs[0, 0].set_xticks(0.5 + np.arange(num_pops))\n",
    "axs[0, 0].set_xticklabels(my_pops)\n",
    "axs[0, 0].set_title('Monomorphic positions')\n",
    "\n",
    "axs[1, 0].bar(range(num_pops), [np.max(list(vals.keys())) for vals in num_reads])\n",
    "axs[1, 0].set_xlim(0, 11)\n",
    "axs[1, 0].set_xticks(0.5 + np.arange(num_pops))\n",
    "axs[1, 0].set_xticklabels(my_pops)\n",
    "axs[1, 0].set_title('Maximum number of allele reads per loci')\n",
    "\n",
    "for pop in [0, 7, 8]:\n",
    "    axs[2, 0].plot(MAFs[pop], label=my_pops[pop])\n",
    "axs[2, 0].legend()\n",
    "axs[2, 0].set_xticks(range(num_buckets + 1))\n",
    "axs[2, 0].set_xticklabels(['%.3f' % (x / (num_buckets * 2.)) for x in range(num_buckets + 1)])\n",
    "axs[2, 0].set_title('MAF bundled in bins of 0.025')\n",
    "None"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "#for pop_data in pop_iter:\n",
    "#    pop_name, allele = pop_data\n",
    "#    print(pop_name)\n",
    "#    for locus_name, vals in allele.items():\n",
    "#        geno_list, heterozygosity, allele_cnts, summary = vals\n",
    "#        print(locus_name, vals)\n",
    "#        break\n",
    "#    break"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2446/NA20127\n",
      "1375/NA12264\n",
      "NA18594/NA18594\n",
      "NA18140/NA18140\n",
      "NA20881/NA20881\n",
      "NA19007/NA19007\n",
      "NA19372/NA19372\n",
      "M005/NA19652\n",
      "2581/NA21371\n",
      "NA20757/NA20757\n",
      "Y105/NA19099\n"
     ]
    }
   ],
   "source": [
    "exp_hes = []\n",
    "for pop_data in pop_iter:\n",
    "    pop_name, allele = pop_data\n",
    "    print(pop_name)\n",
    "    exp_vals = []\n",
    "    for locus_name, vals in allele.items():\n",
    "        geno_list, heterozygosity, allele_cnts, summary = vals\n",
    "        cexp_ho, cobs_ho, cexp_he, cobs_he = heterozygosity\n",
    "        exp_vals.append(cexp_he / (cexp_he + cexp_ho))\n",
    "    exp_hes.append(exp_vals)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6UAAAIYCAYAAACVAkpZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xd4VFXXxuFnpSckdBCpoRcVLCAiSFGRomKnvChNREVs2LA3bCgoFmx07F1UsGBBFFGwgYoIAkpRitQkpO/vjxl48+YLEuUkOwm/+7py6cw5s/aakyHJM3ufM+acEwAAAAAAPkT4bgAAAAAAcOAilAIAAAAAvCGUAgAAAAC8IZQCAAAAALwhlAIAAAAAvCGUAgAAAAC8IZQCQClkZk+Y2c0B1aprZilmFhm+/YmZDQ2idrjebDMbGFS9fzDuaDPbbGZ/FvfYPpjZajM70XcfpZ2ZHWdmy3z3AQAHEkIpAJQw4XCxy8x2mtk2M5tvZheZ2Z6f2c65i5xzdxay1t8GFefc7865ROdcTgC932Zmz+Sr38M5N21/a//DPupIukpSC+dcjeIc+98ws85mtrYI6081s9H57ks2M2dmUb77K0mcc/Occ0133ybsA0DRI5QCQMl0qnMuSVI9SfdKuk7SpKAHKUwgKaXqSfrLObfRdyMo068zAEAACKUAUII557Y752ZK6iNpoJkdKv3vzJeZVTWzt8OzqlvMbJ6ZRZjZDEl1Jb0VXp57bZ7ZsfPN7HdJH+1lxqyhmX1lZtvN7E0zqxwe6//NmO2eSTKz7pJukNQnPN734e17lgOH+7rJzH4zs41mNt3MKoS37e5joJn9Hl56e+Pejo2ZVQg/flO43k3h+idK+kBSzXAfU/fy+FPM7Ls8s9Etw/c3DB/HI8O3a4Z76Zzn+dxT0PEJbz8mXG+bmX2/+3HhbZXNbIqZrTezrWb2hpmVkzQ7T78p4TEjzGyUmf1qZn+Z2Uv5xjkv/Lz/+rvjVFhmFmtmD4SP/QYLLRGP/zf9FfQ6C9/fy8x+DB+bT8ysefj+Pnlqp5hZRnh7m3AvUXn6PMvMvgv/f7yZTQsfy6Xh1/jaPPs2D9fZFh63V55tPc3sJwutSFhnZleH79/zGreC/w29Y2aX5jt2i83s9P39HgDAgYpQCgClgHPuK0lrJR1XwOarwtuqSTpIoWDonHPnSfpdoVnXROfcmDyP6SSpuaRuexlygKQhkmpKypb0cCF6fFfS3ZJeDI/XqoDdBoW/ukhqIClR0qP59ukgqamkEyTdsju4FOARSRXCdTqFex7snJsjqYek9eE+BuV/YDhwTpZ0oaQqkp6UNNPMYp1zvyo0M/2smSVImiJpqnPukzwlCjw+ZlZL0juSRkuqLOlqSa+aWbXw42ZISpB0iKTqkh50zqXm6zfRObde0mWSTg8/t5qStkp6LDxOC0mPSzovvK2KpNp7OU6FdZ+kJpIOl9RIUi1Jt/yb/vLY8zozsyaSnpd0hUKv1VkKhb0Y59zu10xiuNZKSc875xZK+ktS1zw1zw0fR0m6VVKyQq+BruFtkiQzi5b0lqT3FTrWlyr0Pd29NHeSpAvDKxIOVTg457WXf0PT8o3TKnysZv3NsQUA/A1CKQCUHusVCjr5ZUk6WFI951xW+Jw4t49atznnUp1zu/ayfYZz7odwILlZUm8LXwhpP/WXNM45t9I5lyLpekl97X9naW93zu1yzn0v6XtJ/y/chnvpI+l659xO59xqSWMVCmmFcYGkJ51zXzrncsLnvGZIOkaSnHNPS1ou6UuFjm3+mci9HZ9zJc1yzs1yzuU65z6QtEhSTzM7WKFwd5Fzbmv4ezX3b3q8UNKNzrm1zrkMSbdJOjt8rM6W9LZz7tPwtpsl5e7jOV8dnjHcZmbbJC3evcHMLHxMrnTObXHO7VToDYa+/7K/3fK+zvpIesc594FzLkvSA5LiJR2bp48ISc9J+sQ592T47j0hMDwT2y28jyT1lnR3+Hiu1f++eXKMQm963Oucy3TOfSTpbUn9wtuzJLUws/Lhx3+zj+O325uSGptZ4/Dt8xR6IyazkI8HAORDKAWA0qOWpC0F3H+/pBWS3jezlWY2qhC11vyD7b9JipZUtVBd/r2a4Xp5a0cpNMO7W96r5aYpFCzyqyoppoBatQrZRz1JV+ULaXXC/e32tEIzaI+EQ1deezs+9SSdk69uB4WCbR1JW5xzW/9Bj6/nqbNUUo5Cx6pm3h7C4fivfdR7wDlXcfeXpJZ5tlVTaAb36zzjvRu+/9/0t1ve4/Q/33vnXG54e97v2V2SkhSahd3tGUmnmlmiQiF0nnPujzw1846Rf7w14XF2y/saOUtST0m/mdlcM2v3N891j/Br4SVJ54ZDdD/9d+YWAPAvEEoBoBQwszYK/TH9Wf5t4ZnCq5xzDSSdKmmkmZ2we/NeSu5rJrVOnv+vq9Cs0mZJqQqFl919Rep/g8u+6q5XKMzkrZ0tacM+Hpff5nBP+WutK+Tj10i6K29Ic84lOOeel6RwAHpIoSWet+U9lzNsb8dnjUKzqHnrlnPO3RveVtnMKhbQT0HHbY2kHvlqxTnn1kn6I28P4WXGVQr53AuyWdIuSYfkGatCeDntv+mvoOf1P9/78OxsHYW/Z2bWV6GAd3Z4JjVUIFTvC0lnKDQrmTcA/qH/Xbac9/uyXlIdy3PVauV5jTjnFjrnTlNoae8bCgXNghT03KcpNOt/gqQ059wXe3ksAKAQCKUAUIKZWXkzO0XSC5Kecc4tKWCfU8ysUfiP/B0KzVbt/niXDQqdb/dPnWtmLcJh5w5Jr7jQR8b8IinOzE4On7N3k6TYPI/bICk5XxDI63lJV5pZ/XDw230OavY/aS7cy0uS7jKzJDOrJ2mkQrNqhfG0pIvMrK2FlAs/p6Tw9vGSvnbODVXoHNEn8j1+b8dn96xeNzOLNLO48IVzaodn92ZLmmBmlcws2sw6huttkFTFwhd9Cnsi/PzqSZKZVTOz08LbXpF0ipl1MLOYcA//+nd6eDbxaUkPmln18Hi1zGz3Ocf/tL+CvCTpZDM7IfzauUqhJdPzzewIhc4RPt05t6mAx06XdK2kwyS9nq/m9eHjWUvSiDzbvlToTZRrw8e6s0Jv2rxgZjFm1t/MKoQD8O5/NwX5f/+GwiE0V6El48ySAsB+IpQCQMn0lpntVGg26kZJ4yQN3su+jSXNkZSi0IzShDwX5blH0k3hJZZX/4PxZ0iaqtBS2jiFl1M657ZLGi5pokIzTqkKXWRpt5fD//3LzAo6R29yuPanklZJSlfoAjT/xqXh8VcqNIP8XLj+PjnnFil0DuWjCl2gZ4VCF2BSOFh1l3RRePeRko40s/55Suzt+KyRdJpCF5vapND37xr99/fteQrNqv4saaNCF/2Rc+5nhQL7yvD3qqZCwXimQsuyd0paIKlteP8fJV0Sfs5/hJ/D/n6O6HXh47DAzHYo9Jpq+m/6K4hzbplC54Y+otDM7KkKXUAoU6FjVknSZ/bfK/DOzvPw1xVeLhxeqrzbHeHnvSrc7ysKBV2F6/ZS6DzezZImSBoQfi5S6HuxOvxcL1Keixfls7d/Q9MVCsmFfSMEALAXtu9rYQAAgN3M7BOFZq0n+u7lQGJmvyp0tdw5f7PPxZL6Ouc6FUM/AyQNc851KOqxAKCsY6YUAACUaGZ2lkLndn6U7/6Dzay9hT4ztalCS4JfL6hGwP0kKLRi4KmiHgsADgSEUgAAUGKFZ6Yfl3RJvivpSqErMD8paadCgfVNhZbpFmU/3RRamr1B//1oGgDAfmD5LgAAAADAG2ZKAQAAAADeEEoBAAAAAN5E+Rq4atWqLjk52dfwAAAAAIAi9PXXX292zlXb137eQmlycrIWLVrka3gAAAAAQBEys98Ksx/LdwEAAAAA3hBKAQAAAADeEEoBAAAAAN4QSgEAAAAA3hBKAQAAAADeEEoBAAAAAN4QSgEAAAAA3hBKAQAAAADeEEoBAAAAAN4QSgEAAAAA3hBKAQAAAADeEEoBAAAAAN4QSgEAAAAA3hBKAQAAAADeEEoBAAAAAN4QSgEAAAAA3hBKAQAAAADeEEoBAAAAAN4QSgEAAAAA3hBKAQAAAADeRPluwAczC7ymcy7wmgAAAIXF3zcASqtCzZSaWXczW2ZmK8xsVAHbB5nZJjP7Lvw1NPhWg+OcK9TXP90XAADAF/6+AVBa7XOm1MwiJT0mqauktZIWmtlM59xP+XZ90Tk3ogh6BMok3tHefxzD/ccx3H9FcQylA+s4cgyBsoPfK/g3CrN892hJK5xzKyXJzF6QdJqk/KEUBxB+4Oy/wj5fMzvgjk1hcQz3H8dw/3EM998/OS4cRxQl/r7Zf/xM3H8H4uuwMKG0lqQ1eW6vldS2gP3OMrOOkn6RdKVzbk0B+6CM4AcOAAAoa/j7BiXBgfg6LMw5pQVF9fzP/i1Jyc65lpLmSJpWYCGzYWa2yMwWbdq06Z91CgAAAAAocwoTStdKqpPndm1J6/Pu4Jz7yzmXEb75tKSjCirknHvKOdfaOde6WrVq/6ZfAAAAAEAZUphQulBSYzOrb2YxkvpKmpl3BzM7OM/NXpKWBtciAAAAAKCs2uc5pc65bDMbIek9SZGSJjvnfjSzOyQtcs7NlHSZmfWSlC1pi6RBRdgzAAAAAKCMMF8nx7Zu3dotWrTIy9iFVZZOHvaFY7j/OIb7j2O4/ziG+49jGAyO4/7jGO4/juH+4xjuv9JwDM3sa+dc633tV5jluwAAAAAAFAlCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAG0IpAAAAAMAbQikAAAAAwBtCKQAAAADAmyjfDQRl1apVSklJCbzukiVLAqt1yCGHKCKC9wEAAAAAYLcyEUrXr1+vZs1bKLFarUDrRpavps49zwikVtq2vzR+7H0aNmxYIPUAAAAAoCwoE6E0PT1dcRWqKKn/+EDrJgVZbO5k7dy5M8iKAAAAAFDqlYlQCgAAUBYd0eYYrVq1MtCaFhWrilWrB1bv9ltv0eWXjgisHoADD6EUAMqYRYsW6fpbbg+0psUkqGvPUwOrd1D16pr81OOKiYkJrCZQFn3/zULVvHCyFOA1KRIDqySl/vCRvvnu+wArAjgQEUqxxw8//KB169YFXve9994LrFabNm1UuXLlwOoBZdF7772nL1ZuUULzzoHVrNqzpYL8s3PHK49ozN13qmbNmgFWBcqmyMRKsohI320UKCI2QVKm7zYAlHKEUuxx9DHtlFi7aaC/+GIObqIBV94SSK1dWzfqrJO7asrEpwKpB5Rl0VVqK6Hpsb7b2Ktdc5/23QKAA8RFIy7Te+/PCbSmxSSofpMWgdU75+wzNebu0YHVA0obQin2yMrMVPwpN8qigltOFx9YJSn7+/eVnrEtwIoAAKCse/e9D5R66BmKrhTcyowaHaSMgGpl/rlC73/wYYkOpfPmzdOrr78RbNGISF0x8qrAyh16SAsNPf/8wOqheBFKAQBAkVi4cKF27NgReN0PP/wwsFodO3ZUdHR0YPVQMkVXrqWYasm+2yhQbkaqtMl3F3/v3gce0icr/lJM9QaB1azYcYCmfbs1kFouK105k6cQSksxQikAAAjcqlWr1KFjJ1WoF9wSR0mKOaih+l1yXSC1Ujf+rnH33KkLL7wwkHpAWZbQqK3KNe/ou40C5aRt17Yf3vXdBvYDoRQAAAQuMzNTCZWqK+H0YK8EnRBgrdyPn1ZGRlCLMAFg71544QUNGXpBsEUjo5SQmBRYucpVqmrZTz+oXLlygdUsLEIpAAAAABShZcuWKfrQbqrQrrfvVvZq86QLlZqaSigFAAAAgLLIImMUEVv8ga+wzMzb2MF9EjMAAAAAAP8QM6UASpSlS5fq/fffD7zu+PHjA6vVpEkT9ejRI7B6AAAABzJCKRCg7du3a+PGjYHXXb58eWC1atSooaSk4E6KD9r1N9+mD79bodgqdQKrmXh4d9353MeB1HI5WUpfdrNSdwb/MRcAAAAHIkIpEKCuPU7R0mXLFRkTE1jNiISKatPxhEBqZWdm6MjDW+nTD4OfiQyKkxTbvIsSD+kSWM3EwCpJuekp2rVsXoAVAQAADmyEUiBAW7ZuU1KvGxRzUMPAalYMrJKUvnaptvz0YoAVAQAAgP3DhY4AAAAAAN4QSgEAAAAA3rB8FwCAfJxzmjt3rrKzswOtO2fOnMBqxcbGqkOHDl4/Vw4AgCAQSgEAyGf27Nk6u++5SqwZ3PnhMTUa6T8jrg+s3s51v+j9WW/ruOOOC6wmAAA+EEoBAMgnIyNDScmHKL7nqMBqJgRWKeyN25SRkRF0VQAAih3nlAIAAAAAvCGUAgAAAAC8IZQCAAAAALwhlAIAAAAAvCGUAgAAAAC8IZQCAAAAALwhlAIAAAAAvCGUAgAAAAC8IZQCAAAAALwhlAIAAAAAvCGUAgAAAAC8IZQCAAAAALwhlAIAAAAAvCGUAgAAAAC8IZQCAAAAALwhlAIAAAAAvCGUAgAAAAC8IZQCAAAAALwhlAIAAAAAvCGUAgAAAAC8IZQCAAAAALwhlAIAAAAAvCGUAgAAAAC8IZQCAAAAALwpVCg1s+5mtszMVpjZqL/Z72wzc2bWOrgWAQAAAABl1T5DqZlFSnpMUg9JLST1M7MWBeyXJOkySV8G3SQAAAAAoGwqzEzp0ZJWOOdWOucyJb0g6bQC9rtT0hhJ6QH2BwAAAAAowwoTSmtJWpPn9trwfXuY2RGS6jjn3g6wNwAAAABAGVeYUGoF3Of2bDSLkPSgpKv2WchsmJktMrNFmzZtKnyXAAAAAIAyqTChdK2kOnlu15a0Ps/tJEmHSvrEzFZLOkbSzIIuduSce8o519o517patWr/vmsAAAAAQJlQmFC6UFJjM6tvZjGS+kqauXujc267c66qcy7ZOZcsaYGkXs65RUXSMQAAAACgzNhnKHXOZUsaIek9SUslveSc+9HM7jCzXkXdIAAAAACg7IoqzE7OuVmSZuW775a97Nt5/9sCAAAAABwICrN8FwAAAACAIkEoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeFMmQmlCQoIydm7TznnTlblpte92/kdO2nbt/PYdZa1cqPLly/tuBwAAAABKlCjfDQShRo0a+mrBfE2ZNl3Tnxmt9JhEWePjFN+so6LKVyv2fnIz07VrxQLp18+U+tuPOql7dw2dNEHdu3cv9l4AAAAAoCQrE6FUklq2bKkHxz6gsfeP0aeffqpJU6fr9WevVFz1ZLmGHZTQrIMi4xKLbHyXm6P0Vd8qd8U8pSz/Sq2PbqsLr79Ep59+uhITi25cAAAAACjNykwo3S0iIkKdO3dW586dlfHk45o1a5aenjJdHz89VIn1Wyni0B6KTz48sPGyt29U+jdvaNeyz9SgQQNdMPg89e37nKpXrx7YGAAAAABQVpW5UJpXbGysevXqpQoVKigywvTOzDdUvWJNrW42JrhB4iT1kBptq6lrrrhUZ5xxhpKSkoKrDwAAAABlWJkMpc45ffvtt5oyfYaefe55Ka6CrHEH1bxoiqLKV1Vy+sBAx8vNSteuZl/pqvue1IXDL9GJXbvpgsED1L17d8XExAQ6FgAAAACUJWUqlK5cuVIznnlWk6ZO1/aUXYpq0kFxvW5VTNW6RTpuRHScyjXvKDXvqLhdOzT/58/0xciblTFgkM4+6ywNGTRA7du3V0REmbjYMQAAAAAEpkyE0rS0NHU8vquWLl2q+GYdFHXshapYq5nMrNh7iYwvr6QjekpH9FTc9o2a+fNcvdFvsKJz0/Xmqy+rXbt2xd4TAADAAc053x3sXUnuDSgmZSKU/vnnn1q2/FdVHTZFFllynlJUhepKanuO1PYc7ZgzQfPnzyeUAgAAFKNjjz1Gb752q+KbHqfoZp0VU6ORl4mLvJxzylj3s7KXfay0ZfN11uDBXvsBfCs5CW4/RURFlahAmp9Fx/luAQAA4IDz3LQpWr58uaZOm65JUx/U9twIRTbtpPjmnRRVvng/LSFr6x9KX/qxspd9qgqJ8bpk8EANOG+C6tWrV6x9ACVNyU1xAAAAQAAaN26su0bfqdF33qH58+frqUlT9OqzIxV/UAO5RscpoWl7RcQmFMnYObt2Ku3nebIV85S1db369u2rC8a9rqOOOsr7jC1QUnDlHQAAABwQzEzt27fXtMkTtXnDH3rinht0eO5ybXpqiFLfHRfoWC47U6nv3KfNEy9Qu/g/NXX83dq84Q89OeFRtW7dmkAK5MFMKQAAAA44ZrbnS+GAuDruP8EOcpYU/4ARQIF9IJQCAADggOCc0/z58zVx8lS98uqriq2eLDXqqGoXDFREbIKS04Mfs+r5OzV/2Tx9cfkNyto6WP369dXQwQNZvgvkQSgFAABAmbZixQpNnT5dk6ZMV3puhCKadFSF/4wtlgsdRcYnKenwntLhPZW1db1eXvKJnj/5dFVITNAFgwdq4IDzVLdu3SLvAyjJOKcUAAAAZVb/gYPVqnVbPfHeYun4K1T+3IeV1PacYr/yriRFV6qppGP/o4qDnlBmu2Ea/+YXanZoS424fGSx9wKUJMyUAgCAwCUmJmrXts2KnDtFMS26KKZasu+W9shJ3aq0pXOVufIrVap0lu92UMQ+n79AFc68XTHV6/tuZQ8zU1zt5lLt5opq3EGfzX/Td0uAV4RSAChjGjVqpIwl9yslN0exrXooumIN3y1JkpzLVfrKb5Tz47uKiYpQ+fLlfbeEIlSrVi19980iTZoyVVOm3aX02CRZ405KaN5JkeUqFns/LjtTacu/lFbMVepvP+rU007ThS8+o86dOxd7L/CgJJ+7WZJ7A4oJy3cBoIzp06ePfvjuG/VvV087XrhGKTPv0q6VX8u5XC/95KanaOeiN7Rt2iWq+NOruv+aYVr72yolJiZ66QfFp1mzZrr/vnu1Yd0avTTpMXWpkqK/plys1LfuVuqyz+Wys4p0fOecMtYtVcqHj2vTU4PV4K8vNPa6i7Xxz/V68dkZOv744xURwZ9CAOAbP4kBoAxq0KCBxo8bqw3r1+q+q85XhR9e1rZpI5Ty9UzlZqQWSw+Zm1Yr9cMJ2jxpmDpW2ql3X39RP//wvQYPHqz4+Phi6QElQ2RkpE488US98uJz2vDHOt1/zQVK3jBPm54arJSPnlDmxpWBjpeTulU7F7yo7dMuUdRnT+jy09rp5x8W68vP5mrAgAGl6g2RhIREpX42Q9k7Nvpu5X8457Rr1TfK+WWuKlVg1QOA/cPyXQAowxISEnT++edryJAhmj9/vsaMe0jvT7xA5Zodp5iWPQI/z8/lZCtt+QK5H9+V2/6HRlxysS5+a5Jq1CgZS4jhX1JSkgYNGqTjjz9ed44erYlPPy3X89NgB4mT1Fk6eldzPfrIw2rTpk2p/eiN77/7RmMfeljTp49UXN3DFHVYT8XWOdTb88nN3KXUHz9SzpLZqlo+QXfdeKX69+/vpRcAZQehFAAOAGam9u3b68327fXHH3/o4uGX6M3JI1T9jBsU3+TYwMbZ+c59Svn1az3y8MMaMmSIoqOjA6uN0i8lJUWvvfaaHntqkpYs/l4JTTuoxrn3q96uZoGHrKxtfyo99hN1PfVMVUxM0IXnD9aAAeepdu3agY5T1Bo2bKgJj4zXfXeP1rRp0zRm3HjtmCdFHNpT5Vp0UkR0bLH0kbV1vTK+n620nz5Sp06dNer5qerYsWOpDftAcatVq5ayfnpUKcpRbPMuiq5ax3dLkiSXk6Vdvy6SW/6pIi30ZrYPhFIAOEBkZmbqtdde071jH9KvK1erSqeB2tDyUUmPBjdIb6namx115dXX6vMFX+mqKy5Tq1atgquPUic3N1dz587VExMn660331S5uodIjTup2rArZFExRTZudMUaij62r1y7PspY95PGvfaJ7rznXh1x5FG6ZNj5OuOMM7z98fVvJCUlacSIERo+fLjmzJmje+4fp68mzlD8oScqrlWPIvl4E+ec0ld/q5wls5X5xzJdMPR8Xfbig6pXr17gYwFl3dChQ3X00Udr8tRpmj7jVqWXqyJr0lHxzToqMqFCsfbinFPmH78oa9kn2vXzZ2rWvLmGXzFEZ599trfTGwil2MPMtOu37xXf4CiZlazTjXMz0pT15zJFVDvYdytAqbN+/XpNePwJTXjiSUVUqq2IQ7qrcpe2sohIJaefE/h45bpJcR22atbi9/X68SepUcMGGnXVFTrzzDOZOT2A5Obm6tbbbteTT09SdlS8rGlnVRk0QZGJlYq1j9BHbxwi1T5ECZ2GavnyBbpi9MMadtFw9Tr9dI0bc68OPrj0/G6JiIjQSSdHg7LuAAAgAElEQVSdpJNOOknLly/XuPGPaNq0y1W+RUfFdbk4sHGytqxT6tv3qHJinK6/ZqT69+9fqkI8gnVQ9arK/uBtpUVGK77R0bKISN8t7ZHx5wplfTdTlSpX9t3KPrVs2VIPjRurB8bcpw8//FBPTpyidydfpHLJLWWNOym+YRtZVNH9nszesVG7fvpEOcvmqlxMhEYMGaSBL4xX/fr+Py6JUIo9ZkyfphtvvVObP5uqqMN6KOGQ4xUR6/cXUNaWdcr4fpbSfvpYJ3Y9SbfedL3XfoDSwjmnzz//XGPGjdcHH7yvcs07Kq7XrYqpVjwzHJHlKimxXR+VO/osrV2+QJfcdK8uvvRyjRh+kYZfdBHnmB4Ali9frrHjH1XFM29TwkENfLcjSYqIjlW5Fp2kFp0Uu/MvvfPew2r38su67LLLfLf2jznntHbtWq1ctUqyCOXGVdTquP8EN0BNScOkCg/HaOWq1dq6dSuh9AD2xGOPqOurr+quMWP1+2dTFHNYD8Uf1lWRcX5m1VxOttJ+ma/cH2YrMm2Lrr7iMg27YKiXXv6NqKgodevWTd26ddOOHTv2nNbw44cTlNC0veLanK2o8tUCGy/tly/kfpitXRtW6pxzztGwe57XMcccU6KW3xNKsUefPn3Uu3dvzZs3T/fcP1ZzJ12ghBbHK7ZVT0VXKr53kZ3LVfqqb5WzZJayNqzQhcMu0GWvPlrqzgMCfPn222/Vp/8AbdiyQ1GHdVe1oU8rIracl14sMkrlmnWQmnVQ5qbVmvDObN0/tqlOOeVUPT9jqqKi+DVUlsWWS1JMCQmk+UUlVVFMldL3eyUtLU0zZszQfWMf0paUdEUe1lNVh05UREy8ktP7BT5e4llr9PSHs/TQwy10UreTdN1VV6pdu3Yl6o9ZFL2YmBj169dP/fr105dffqkxYx/SrEnDQhfNa3WyYqrWLZY+ctK2a9fi95Sx5F01bdJYN9x/m3r16lWqf5eUL19eZ555ppxzGn3PGK3/eZ42nzwv2EFaShFvpqhPv/66YMigEhdIJUIp8jEzdezYUR07dtRvv/2m8Y88qomTrlNMzWaKPKyn4uq1KrIXcW5GmtJ+/FDZS2areqXyuv7akerbty8fHQH8Q7NmzdKfsbVUaeDdJWopfky1ZMWccLES2p+ntyddqI0bN6pmzZq+2wJKhdWrV+uhhx/RpMlTFFuruSKPPFcVivB38m7RVeoousuFSjj2XH2+ZI66n9FbNatX1Q3XjFSfPn0UG1s8F1pCydG2bVu9+tLz+uOPP/TYhMf16IRblVm1niIO66H4Bq2LZGlv5oZflbl4ltJ++UJnnnmmrhn/fqm/XkF2drbmzJmjJydO0Xvvzla5+q1kh56p6r3aKDk9+CW8NS/cqPeXztXss/srIUq6YMggDRxwnho0KBlvHJacv1ZQ4tSrV0/jHrhff65bo3uuGKRy3z6r7TMuVcqSDwIdJ2fXTqV8MlGbJ12gYxI26a2XntGyHxfzWYbAfohMqFCiAmleEXGJioouugvcAGXJ559/rm4nn6oWLQ/XM/NXqXyfMUo89QbFJx9erDMdEbHllNj6NFUa9Li2NDtdI+9+VNUPrq3rb7hRmzdvLrY+UHIcfPDBGn3nHdqwfo0euvlyJS15RVunDlf2jk2BjpP65h1KfeN2XXFmR/2+6lc9M21KqQ6kixcv1uVXjlS1GjXV/+KR+nRHZVUZ8oTKnTxKCU2OLbJzSqPKV1dS23NU4bxH5Lpcpsdmf6NDj2itI9seq4kTJ2r79u1FMm6h+/M6OkqFhIQE9e/fX7m5ubrmulFyy+ZqdZspwQ0QJ6m7dLQdqmtHXq4OHTqUuCUFAAD40OWEE5TUcfCeJbq+mUUooWEbqWEbZf21RhNenaqt23boiQmP+G4NnqxYsUIffjxXf6xfp3KNjtba6pcHO0AfqdyDkfrwk091dJvW6tq1a6n8O/Grr75S/4FD9OfGzYpt1lkJp9/h5WNhzEyxBzdR7MFNVK7jYK1d+bWuHz9dl15+pU7pdZq3U2sIpfhbv/32W57lQs0U3/1qxdVrpeT0YH8Y5GakKS3pI/Xqc56qVUzS9deMVL9+/ZgpBQAc0LKzspR4eI8SdbXT3aKr1FFU8lHalZHuuxUUs5ycHL3zzju65/5xWvLDj4o7rJuqDHxMkYmVlFwEL4cqQzO0ZOlc9Tn/EpWPMV078nINHDjQ28eX/BuzZ8/Wn3F1VXnoPSVmJZNFRiuh8TFS42MUt2un3po4TFu2bFH16sF/xNS+FOqImFl3M1tmZivMbFQB2y8ysyVm9p2ZfWZmLYJvFcXFOae5c+eq+ym91PywVprx+Uol9RmjcqcU3XKhiNgEJR55iioMeFQ7Duut6x54WgfVqqNrrhultWvXBj4eAAAA/plt27Zp7Lhxqp3cUIMuG6XlFVurytCnlXhsvyL9uKeI6FgltjxJ5f/zoDLaDtFtT7yoGrXq6NLLr9TKlSuLbNygRSZULDGBNL/I+CRFerxg1D5HNrNISY9J6ippraSFZjbTOfdTnt2ec849Ed6/l6RxkroXQb8oQrt27dJzzz2ne+4fp8070hR1WI9iXy5kFqH4BkdJDY5SzJZ1mjJ3th6bcIhOOLGrRl19pY499thSuWQDAACgtPrtt99059336Lnnnle5hq0V1flSla/VrNj7MDPF1T1MqnuYYrdv1PMLZ2nyEa3V7th2uv2mG9S+ffti7wnBKEwcPlrSCufcSkkysxcknSZpTyh1zu3Is385SS7IJlE8mrY4VGnx1RXZqq8qJB/u/Z2c6Mq1FN15qBLa/Ufzf/hQPc7oo8suGqrRd9zmtS8AAIADyZVXX6cPl21W5QGPKCqpiu92JElRFaorseMgJbTrq28Xv69TTz9TWzZt8N0W/qXChNJaktbkub1WUtv8O5nZJZJGSoqRdHwg3aFY/bFurWpe9oAsqmRdFTMiNkFJR52qnVGx+nXVKt/tAAAAHFAyMrMUU/+oEhNI84qIjlNC807atvBl361gPxRmKqygtZL/bybUOfeYc66hpOsk3VRgIbNhZrbIzBZt2hTs5aIBAAAAAKVPYULpWkl5r1dcW9L6v9n/BUmnF7TBOfeUc661c651tWrVCt8lAAAAAKBMKkwoXSipsZnVN7MYSX0lzcy7g5k1znPzZEnLg2sRAAAAAFBW7fOcUudctpmNkPSepEhJk51zP5rZHZIWOedmShphZidKypK0VdLAomwaAAAAAFA2FOrDaJxzsyTNynffLXn+//KA+wIAAAAAHABK5qe3AgAAAAAOCIRSAAAAAIA3hFIAAAAAgDeEUgAAAACAN4RSAAAAAIA3hFIAAAAAgDeEUgAAAACAN4RSAAAAAIA3hFIAAAAAgDeEUgAAAACAN4RSAAAAAIA3hFIAAAAAgDeEUgAAAACAN4RSAAAAAIA3hFIAAAAAgDeEUgAAAACAN4RSACVKVGSEsv5cLpeT5buV/8c5p/R1SxVh/OgEAAAICn9ZAShRRt92i1olpWnr1OFKWfy+XE6275YkSelrflDKqzcrZuEMzZg2xXc7AAAAZQahFECJ0rx5c3360Qea/cbLarTjO22bPkIpP3wkl5vjpZ+M9cu08/XbpE8e0/03XKrVK5bpjDPO8NILildudsl4Q6QgzrkSuZoAAIB/g1AKoERq3769Fnw2V2+8MF11Nn6u7c9crtSl8+RcbrGMn7nhV6XMvEtZ7z+g0Vecr99XrdCQIUMUFRVVLOPDr1atWilq53rtfO1Wpa/9yXc7ezjntGv1d0p55UbF56SoRYsWvlsCAGC/EUoBlGhdunTRN19+oRenPKHqv32gHc+NVNovX8g5VyTjZW5arZR37lP623fp5ov6au3qlbr44osUExNTJOOhZGrQoIF+X7lC9117kfTxI0p5/TZlrFvqtaf03xcr5dWbFL1gssbffo1WLV+mmjVreu0JAIAgEEqBAEVHRyl99bcl5jzIvFxOljJ++1Yx0dG+W/nHzEzdunXTD98u0vQJ41Rx2UzteP5qZe/YFOg4aR88rLQ3bte1A3pp7W+rdPlllykuLi7QMVB6REdHa+jQofp91QrdfdUFyv1ovFLeuF0Z65cVax/pa37QzldvVuTnT+vBW0Zq9YplOvfccxUZGVmsfQAAUFRYhwYEaMbkp3XJ5SP18/QRim7TW+VadJJF+P3D0eVkK/XHj5S58BW1OrSZHnlkgtd+9ld8fLzi4uLktqVobfXLgy1+qlRlQrzKJcTzBz/2iImJ0bBhwzRo0CBNmjRZt9xxp7Iq1lF0m96Krdm0yMZNX/ODsha+pJhdf+mBO27Veeedx/JxAECZxG83IECtW7fWl59/qo8//lhXXXeDfv36VcW06a2E5sfJivljRFxujlJ/+kRZC19W04b1Nfa1F3TccccVaw9B+vTTTzXy2uu1fPUaxRzdRxVOOE7J6cEHx5gey3XrozN0+1336K7bb9WgQYMUXQpnlxG8mJgYXXzxRRoyZLAmTZqsW+8crayKdRRz7HmKqZYc2DgZf65Q5vwZiknbpDG336IBAwbwGgQAlGmEUqAIdOnSRV9/OV8ffPCBrhp1g9Z8/aqi2/RWQtNjizycOpertKXzlLnwJTWsc7DGvjBdXbp0KdIxi9KCBQt01XU3aMnSXxTT5hxV7Hhdkc4+xx7cWLGn3az0tUs16v4ndeudd+ve0berf//+zJ5CkhQbG6thwy5QdHSUhg+/RFl9Aj7XNDn0NaPho+rfv78iIjjTBgBQthFKgSJiZjrppJO0uGtXzZ49W1ddd4P+/PpVRR/dR/GN2srMAh3PuVyl/fKFsr56UXUOqqSx055S165dAx+nuHz99de6etQNWvTtYsW2OVuVBl4uiyy+H1lxtZsrrvbtSv99sa68faxuuu1Ojbn7TvXu3ZuQcADLycnRc889p+tvvk27opJU5ZzblZzeMtAxnHNKX/WNsqbdq1vuuEv33XWHzj77bF53AIAyi99wQBEzM/Xs2VM/Lf5WUx+9X+V/fE07XrhauVkZgY2Rm5Gm7c9dpcq/vK1nn3pYi79eqJNOOqlUBtIVK1ao+ymnqVPXHloS0VCVBz2uxFbdizWQ5hVXt6USz75LGW0GaPioO9S4+aGaPXu2l17gT05Ojp599lklN2qqK257QFntzlfiWaMVVzfYQCqFfmbENzhKSb3vU9rh/XTRtbepUfND9corryg3t3g+EgkAgOLETClQTDIyMrRy1Spt+WuzIms01e9Jg4MrHifpfKnK4wlavXq1MjMzFRsbG1z9YjTqplv1xfpsVR78hCKiS8ZzMDPF1z9ScclHaMeyz9W773+0c/tW322hGOTk5Oill17SqJtu1U4Xq+i2g5VYr1WxvOFjZopv2FpxDY5SyspFGnb1zRp1820ac9cdOv3005k5BQCUGYRSoIhlZGRo4sSJuuWO0bKqDRR3yo2KOaihktODHyum+38v0jP6tls0ePDgUneBlKzsbEXXal5iAmleZqb45MOVmpvjuxUUsdzcXL388su67qZbtSM7StFHn6ek5CO8rD4wMyU0bCPXoLV2/PqVzr/iBl130y26/+7ROu2000rliggAAPLibVagiGRlZempp55SnfoNdeujMxTT/VolnnqDYg5qWGRjxh7cWImn3azIE67Q9WOfUt36jTRlyhRlZ5e8z00FSrI5c+Zo6CVXaNcR/1FS73sVX/9I7+HPzJTQqK3K93tA25ufqXMHX6AFCxZ47QkAgCAQSoGAZWdna8qUKapbv5GuH/uUIk+4Qomn3azYg5sUWw+xtZor6Yzb5ToN19V3Pax6DZvomWeeUU4OM3xAYaSmpiqhZiPFNzjKexjNz8yU0LitEg5KVmpqqu92AADYb4RSIECzZs1ScqOmuvquh+U6DVfSGbcrtlZzb/3E1TlUiWfdqexjh+qyW+5TgybN9cEHH3jrBwAAAMiPc0qBAF121XVKP/R0JR7SpcTMrpiZ4uu1Ulzdlkr5/j2NvO5GLena1XdbAAAAgCRmSoHAxVSrV2ICaV5mpuiq9eTkfLcCAAAA7EEoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeEMoBQAAAAB4QygFAAAAAHhDKAUAAAAAeFOoUGpm3c1smZmtMLNRBWwfaWY/mdliM/vQzOoF3yoAAAAAoKzZZyg1s0hJj0nqIamFpH5m1iLfbt9Kau2caynpFUljgm4UAAAAAFD2FGam9GhJK5xzK51zmZJekHRa3h2ccx8759LCNxdIqh1smwAAAACAsqgwobSWpDV5bq8N37c350uavT9NAQAAAAAODFGF2McKuM8VuKPZuZJaS+q0l+3DJA2TpLp16xayRQAAAABAWVWYmdK1kurkuV1b0vr8O5nZiZJulNTLOZdRUCHn3FPOudbOudbVqlX7N/0CAAAAAMqQwoTShZIam1l9M4uR1FfSzLw7mNkRkp5UKJBuDL5NAAAAAEBZtM9Q6pzLljRC0nuSlkp6yTn3o5ndYWa9wrvdLylR0stm9p2ZzdxLOQAAAAAA9ijMOaVyzs2SNCvffbfk+f8TA+4LAAAAAHAAKMzyXQAAAAAAigShFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4E2U7waCkJSUpOy0HVr34Fm+W9mr3Oxs1br8HN9tAAAAAECJUiZCabVq1bT1r83KyckJtG5CQoLS0tICqxcfHx9YLQAAAAAoC8pEKJWkmJiYIqlLkAQAAACAosM5pQAAAAAAbwilAAAAAABvCKUAAAAAAG8IpQAAAAAAbwilAAAAAABvCKUAAAAAAG8IpQAAAAAAbwilAAAAAABvCKUAAAAAAG+ifDcAAAAAFJX4uFjt+GiCsmPjfbdSoKy0nYpPrum7DcArQikAAADKrHffflMrVqwItOaJJ56oOXPmBFavWbNmgdUCSiNCKQAAAMqsevXqqV69eoHXPeGEEwKvCRyoOKcUAAAAAOANoRQAAAAA4A3LdwEAAACgiLnsDOXs2uG7jb1yznkbm1AKAACKRE52lrK2rPPdxl7lpKf4bgHAAaJ58+bKGfeQtv78UWA1U1NTVa5cucDqVatWNdB6/wShFAAABK569eqqXaOatr17T6B1/9ywQTUOOiiQWglmatWqVSC1AODv9O7dW7179w60ppkpZfvWQGv6QigFAACBq1Spkn5a/F3gdc1M639fFXhdAIA/XOgIAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4A2hFAAAAADgDaEUAAAAAOANoRQAAAAA4E2U7wYAAMGKjIxU1ooF2pW60Xcre5WeukMREbwvCgAACKUAUOYMHz5c9evXD7Rm37599cILLwRWr0qVG1WjRo3A6gEAgNKLUAoAZUz58uXVp0+fQGv27ds38JoAAAASoRQAAKBEy/zjF8lK5nL37O0bpNqVfLcBoJQjlAIAAJRQJ/c6Q8u/fzbQmr/88ouaNGkSSK0Kkroef00gtQAcuAilAAAAJdRbr78SeE0z08+Lvw28LgD8WyVzLQgAAAAA4IDATCn2MDOlffK0FBHpu5UCZW76XRHHHem7DQAAAAABIpRij7dmvqmVK1cGWnP48OGaMGFCYPW6desWWC0AAAAA/hFKsUdRBL7hw4fr4ov/r717j7K1LusA/n04aIAXSsMoRSAkDRQVjlinixW6ktWFKBJoVejSXJmkrZamaSW5upklamBEZSJrGXTTwLiUGZqIwkEFAU3wlicyj9fSUAR+/fG+E9vNwAzsPec3e+bzWWvW2e9lv/PMc2a/7/6+7+/d86y5bxcAgF3jfve7T7588dlpV53Xu5Rl3XbbLdlrr716l8EMhFIAAOBOnfma03L99dfPdZtbt27N9u3b57a9fffdd27bYtcTSgEAgDt1//vfP0ccccTct7sW22Qx+fRdAAAAullVKK2qJ1fVv1XVDVX1wmWWf29Vvaeqbqmq4+ZfJgAAABvRiqG0qrYkOT3J0UkOSXJiVR0ytdq/J3lqkjfMu0AAAAA2rtXcU3pkkhtaax9Jkqo6J8kxSa5bWqG19rFx2W1rUCMAAAAb1GqG7z44yScmpneM8wAAAGAmqwmltcy8dk++WVU9s6q2V9X2nTt33pNNAAAAsIGsJpTuSLLfxPRDktx4T75Za+3M1trW1trWffbZ555sAgAAgA1kNaH0iiQHV9WBVXXvJCckOW9tywIAAGAzWDGUttZuSXJykouTfCDJX7XWrq2ql1bVjyZJVT2uqnYk+ckkf1JV165l0QAAAGwMq/n03bTWLkhywdS835h4fEWGYb0AAACwaqsZvgsAAABrQigFAACgG6EUAACAboRSAAAAuhFKAQAA6EYoBQAAoBuhFAAAgG6EUgAAALoRSgEAAOhGKAUAAKAboRQAAIBuhFIAAAC6EUoBAADoRigFAACgG6EUAACAboRSAAAAuhFKAQAA6EYoBQAAoBuhFAAAgG6EUgAAALoRSgEAAOhGKAUAAKAboRQAAIBuhFIAAAC6EUoBAADoRigFAACgG6EUAACAboRSAAAAuhFKAQAA6EYoBQAAoBuhFAAAgG6EUgAAALoRSgEAAOhGKAUAAKAboRQAAIBuhFIAAAC6EUoBAADoRigFAACgG6EUAACAboRSAAAAuhFKAQAA6EYoBQAAoBuhFAAAgG6EUgAAALoRSgEAAOhGKAUAAKAboRQAAIBuhFIAAAC6EUoBAADoRigFAACgG6EUAACAbnbvXQBsJAcesH/e8rrnpqp6l7Ks1lqOPPa43mXcpS277ZavXnNx/vfGq3uXsqx2y1ez227O520GN//P55KPXNm7jDt185e+0LsEAJgLoRTm6B8vOD+ttbluc8uWLbn11lvntr31GpiXnPry38vb3va2uW7zpJNOyllnnTW37R100G/NbVusT4cffngOO+jBufmTl8xtm+++/PI8/sgj57a9PQ95WA499NC5bQ8Aeql5v4Fera1bt7bt27d3+d6rVVVzDxibjR7OTg9np4ez08PZ6eF86OPs9HB2ejg7PZzdIvSwqq5srW1daT1j0AAAAOhGKAUAAKAboRQAAIBuhFIAAAC6EUoBAADoRigFAACgG6EUAACAboRSAAAAuhFKAQAA6EYoBQAAoBuhFAAAgG6EUgAAALoRSgEAAOhGKAUAAKAboRQAAIBuhFIAAAC6EUoBAADoRigFAACgG6EUAACAboRSAAAAuhFKAQAA6EYoBQAAoBuhFAAAgG6EUgAAALoRSgEAAOhmVaG0qp5cVf9WVTdU1QuXWf51VXXuuPzdVXXAvAsFAABg41kxlFbVliSnJzk6ySFJTqyqQ6ZWe3qSz7XWHpbk1CQvm3ehAAAAbDyruVJ6ZJIbWmsfaa3dnOScJMdMrXNMkrPGx3+T5KiqqvmVCQAAwEa0mlD64CSfmJjeMc5bdp3W2i1JvpDkgfMoEAAAgI1r91Wss9wVz3YP1klVPTPJM5PkoQ996Cq+9dq4OxdxV7tua3f4cTc0PZydHs5OD2enh7Nbix4mm6uPd3dwld/F5Xk9z04PZ6eHs9uMPVxNKN2RZL+J6YckufFO1tlRVbsn2TvJZ6c31Fo7M8mZSbJ169ZunVnv/ymLQA9np4ez08PZ6eHs9HB2ejgf+jg7PZydHs5uM/ZwNcN3r0hycFUdWFX3TnJCkvOm1jkvyUnj4+OSvLVtxm4CAABwt6x4pbS1dktVnZzk4iRbkry2tXZtVb00yfbW2nlJ/jzJ2VV1Q4YrpCesZdEAAABsDKsZvpvW2gVJLpia9xsTj7+c5CfnWxoAAAAb3WqG7wIAAMCaEEoBAADoRigFAACgG6EUAACAboRSAAAAuhFKAQAA6EYoBQAAoBuhFAAAgG6EUgAAALoRSgEAAOhGKAUAAKAboRQAAIBuhFIAAAC6EUoBAADoRigFAACgG6EUAACAboRSAAAAuqnWWp9vXLUzyce7fPPV+8Ykn+5dxILTw9np4ez0cHZ6ODs9nA99nJ0ezk4PZ6eHs1uEHu7fWttnpZW6hdJFUFXbW2tbe9exyPRwdno4Oz2cnR7OTg/nQx9np4ez08PZ6eHsNlIPDd8FAACgG6EUAACAboTSu3Zm7wI2AD2cnR7OTg9np4ez08P50MfZ6eHs9HB2eji7DUCTOJAAAAkMSURBVNND95QCAADQjSulAAAAdLMpQ2lVHVtVraoeMU7vVlWvrqprqur9VXVFVR1YVc+tqldOPO9PquotE9O/WFWv7vEzrAdVtW9VnVNVH66q66rqgqr6tqq6qareN/H1s+P6X5x6/lOr6rQ+1a8Pd9HDa6bWO6Wqnjc+fl1VfXTs7Qer6iV9ql8fZuzhVVX1oap6fVU9uM9PsH5U1TdV1Ruq6iNVdWVVXTbuL7+vqt48rnOH121VXVJVG+LT/+alqr5YVQdM7A+vq6ozqurRE/vGz068lt+y8lY3tmWOEV9fVZ+pqhqnv3M8dj9knN577OFu42v6uHH+A6rqvVX1tF3/U/Q19ufsiendq2rn1Ot359Qx+pCq2jq+B7r3uN5B437g/r1+ll5W2cPTxse7VdVZVfXaGnysqr5xXHbE+Pp+bJ+fpL+qeuDE79knq+o/JqZfUlXXVtXV4/Tjx+c4nkwYf6/eUVVHT8x7SlVdVFW3jr27pqrOr6qvH5cfMP0eaBFsylCa5MQk70hywjh9fJJvSXJYa+1RSY5N8vkk70yybeJ5j0myd1VtGae3Jbl0l1S8zoxvEt6Y5JLW2kGttUOSvCjJNyX5cGvtMRNfr+9a7Dq1Qg9X8vzW2mMy/E6eVFUHrmGp69YcevjoJA9P8t4k/7L0hmwzGnv5piRvb619a2vtiAz7yIf0rWzhfXh8rR6W5JAkBy3tG5Ocl/G13Fp7Ytcq16HW2ueTfDLJt4+ztmV4rS4dl78jybtba7ctPaeq9k5ycZIzW2t/sQvLXS++lOSRVbXnOP2kJP8xtc65U8fo61pr25O8PcnzxnVOT/Li1tp/75qy15XV9HBpn3lGknsleUabuB+uqg5L8jdJjm+tvXftS16fWmufmdjfnZHk1PHxs5I8OcnhrbXDkjwxySc6lrpujb9XP5/kFVW1R1XdJ8lvJ3l2kpvG/j4yyWfHeQtr04XSqrpvku9K8vTcHkq/Ocl/Lh3YWms7Wmufy3Dw+7aq2nM80P1vkvcledT4vG0Zgutm9P1JvtpaO2NpRmvtfbFTuTvm0cM9xn+/NM/CFsjMPWyDUzO8+T16pfU3sB9IcvNULz/eWvujjjVtGK21WzIcLx7Wu5YFc2luD6Hbkpw6NT15DL5vkguTvKG19se7rML158IkPzQ+PjHJX67yeS9K8oyq+pUk92qtrfZ5G9FqeviqJA9M8rOTJ0YynER5U5Kfaa1dvqZVLq5vTvLp1tpXkqS19unW2o2da1q3WmvXJDk/yQuSvCTJ61trH55a7bIkCz3ia9OF0iQ/luSi1tqHkny2qg5P8ldJfmS8BP6HS0MtxjcR70vyuIxnZJO8K8m2qvqWDB8UtVlD2COTXHknyw6aGhr0PbuysAWy6h5mOEs26eXj/B1JzmmtfWotC13HZunhtPckecRcq1ssh2bowWocP9VbQ61WUFV7JTkqyft717JgJkcsfWuSv87tv2/To5VekeQd40mmzeycJCdU1R4ZrtC/e2r58VPH6D2T/78y/bIkv5vkF3ZpxevPSj38qSRHJDlhfK846e+TnNxae8fal7mw/jHJfjXcPvOaqnpC74IWwG9m+L07OsnvTy4YR3AelWH0zcLajKH0xAw7m4z/ntha25FhCN+vJrktyT9X1VHjOktnabdlOAtx2fj4u7J5r5KuZHr47r/exbo+/nl5X9PDDMNeJi0N3903yVFVte2Om9j0VurhtNoVRS2Kqjq9hntur1hm8blTvd2+q+tbIAeNwf3SJP/QWruwd0EL5tIMJ4IPTPKx1tqXM4ycvG+GUDB5JeqtSY6pqgd1qHPdaK1dneSADO93LlhmlenhuzdNLDs6yX9lGGq+aa2ih+9Jsn+SI5dZ9pYMV5y3LLOMJK21L2Z4/T4zyc4k51bVU7sWtc611r6U5NwkZy9dYU6y53h8+UySByT5p171zcOmCqVV9cAMQ9T+rKo+luT5Gc4YVmvtK621C1trz0/yOxmuqCa3n6X9zgyB9AMZdtab9n7S0bUZdih3x01T9+w9IMmn51fSwrknPfwa4479kiTfPY+CFtDMPZzw2Ayv783q2iSHL0201p6d4czrPt0q2hiWTo48trV2Su9iFk1r7fok35DkRzIcg5NhdMTTknx03AcuOSfJHye5oKrut0sLXX/OS/IHWf3Q3VTVDyfZO8kPZhiNs9ca1bYo7qqHH0zylAxh6tCpZSeP/75mDWtbeK21W1trl7TWXpKhZz/Ru6YFcNv4teSm8cTw/knuHfeULpTjMozD3r+1dkBrbb8kH03yveNw3FTVbhmGanx8fM47Mwzd3ae19qnxhuOdSY7J5r5S+tYkX1dVP7c0o6oel+GFcWfeluSnx3X3zLBD/5e1LHKduyc9/BpVtXuSxyeZvrdgs5hHD6uqnpPhHpeL5l/iwnhrkj2q6lkT8zb7m1LWh8uSPDe3h9LLkvxSljkGt9ZemeSfk7xxM39wWZLXJnlpa21Vw8XHY/IfJnn2+Jy/T/LiNaxvEdxlD1tr78xwW8g/VNVDJxbdluEK68Or6qVrX+biqaqHV9XBE7Mek9vfd3M3tda+kOQ5SZ5XVffqXc89tdlC6YkZPqlz0t8meV2S88ePT746yS1JTkuS8QOPdma4irDksiQPSnLVGte7bo3h/NgkT6rhT3Fcm+SUJDfmjveUPmd82nOT/Pg41OBdSf66tfb2HvWvByv0cCVL95ReneEetb9bs0LXsTn08KokH8pw3/j3t9ZuXrNi17mxlz+W5Ak1/BmDy5OcleGDFbgbxpNFX1lxRabtVVU7Jr5+eZx/aZL9cvsw8csy3F+67Inh1toLMnzY2dnjieZNZ/zAxlfdyeLpe0q3Jfn1JG9qrV03rnNKhnsqD76TbWx4K/RwaZ03Z7jX76JxNN7S/K9kuHjxo1W10Fev1sh9k5xVw5/KujrDCMRT+pa02MZPeb4qt3+I68KpiU+wBgBmVFWPTvKnrbXl7jcDAKZsyjOIALAWqurnM9yD9mu9awGAReFKKQAAAN24UgoAAEA3QikAAADdCKUAAAB0I5QCAADQjVAKAABAN0IpAAAA3fwfAUQd37DVNpYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1152x648 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#XXX correct this\n",
    "fig = plt.figure(figsize=(16, 9))\n",
    "ax = fig.add_subplot(111)\n",
    "bplot = ax.boxplot(exp_hes, notch=True, patch_artist=True)\n",
    "#print(bplot.keys())\n",
    "#bplot['boxes'][0].set_facecolor('red')\n",
    "#bplot['boxes'][0].set_edgecolor('red')\n",
    "ax.set_title('Distribution of expected Heterozygosity')\n",
    "ax.set_xticks(1 + np.arange(num_pops))\n",
    "ax.set_xticklabels(my_pops)\n",
    "None"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
